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Die magnetischen CGS-Einheiten 
und das Internationale Einheitensystem 


Von 
EUGEN FLEGLER, Aachen 


(Eingegangen am 12. Dezember 1959) 


Inhaltstibersicht: Die Einheiten der Dreier-Systeme. — Der Ubergang zum Vierer-System. — Die 
Einheit Maxwell und das Internationale Einheitensystem. — Die Einheiten Gau®B und Oersted. — Weitere 
Entwicklung. — Literatur. 


Von den Einheiten der CGS-Systeme haben die elektromagnetischen Einheiten eine 
besondere Bedeutung gewonnen, einmal, weil die Einheiten der magnetischen Gré8en heute 
noch haufig verwendet werden (wie z. B. die Einheiten Maxwell, Gau8 und Oersted), und zum 
anderen, weil die elektrischen Einheiten des internationalen Einheitensystems (wie z.B. 
Volt, Ampere usw.) ihr Dasein dem elektromagnetischen CGS-System verdanken. Viele Jahr- 
zehnte sind die genannten Einheiten wie Maxwell, GauB, Oersted, Volt und Ampere in Gleichun- 
gen nebeneinander benutzt und durch Gleichheitszeichen miteinander verbunden worden. 
Erst neuerdings werden Bedenken gegen derartige Gleichsetzungen laut und es wird die 
Forderung erhoben, in Beziehungen zwischen den Einheiten Maxwell, GauB und Oersted einer- 
seits und den Einheiten Volt und Ampere andererseits nur das Entsprichtzeichen zu ver- 
wenden [1]. Ware diese Forderung berechtigt und ihre Erfillung unausweichlich, so kénnten 
z.B. die Einheiten Gau8 und Volt nicht mehr gleichzeitig in dieselbe Gleichung eingesetzt 
werden, sondern sie diirften nur in einer Entsprichtbeziehung verwendet werden. Es ist zu 
bezweifeln, ob der praktisch rechnende Ingenieur mit derartigen Entspricht-Beziehungen, etwa 
bei der Berechnung elektrischer Maschinen und Transformatoren, etwas anfangen kann. 


Die Einheiten der Dreier-Systeme 


Wir betrachten die Beziehungen zwischen den Einheiten der elektrischen Spannung und 
des magnetischen Flusses im elektromagnetischen CGS-System. Dieses System beruht 
auf drei mechanischen Grundeinheiten (cm, g, s) und gehért damit auch zu einem auf drei 
mechanischen Grundgr6Ben (Lange, Masse, Zeit) beruhenden Dreier-GroBensystem. Ent- 
sprechend der physikalischen Gleichung U = d@®/dt besteht im elektromagnetischen CGS- 
System die Einheitenbeziehung 


@) 
Te eles 
1 [ ive 1 alle (1) 
‘und mit 

eer ees eee CUS oes a (2) 

1 [Jem = 18 (3) 
ergibt sich 

IRC ae ticen == a CINIose BS 2 (4) 


Das aus dem elektromagnetischen CGS-System abgeleitete Quadrantsystem ist ebenfalls 
ein Dreier-System. Die Spannungseinheit dieses Systems ist die sogenannte absolute prak- 
tische Einheit, fiir die sich schon 1873 der Name Volt findet. Er wurde 1881 international 
angenommen [2]. Hiernach ist 

PV =108(0) 
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und mit Gl. (4) 
1 M1079 Vs. (6) 


Alle drei in Gl. (6) verwendeten Einheiten (M, V und s) sind hier noch echte Dreiersystem- 
Einheiten. 


Der Ubergang zum Vierer-System 


Der Ubergang von den Einheiten der Dreier- zu denen der Vierer-Systeme ist (fiir unsere 
Betrachtung) gekennzeichnet durch zwei Schritte: 


a) Durch die Einfithrung einer vierten (elektrischen oder magnetischen) GrundgréBe und 
damit einer vierten Grundeinheit. Verbunden war dieser Ubergang zum Vierersystem mit der 
gleichzeitigen Einfiihrung der Influenzkonstanten e, und der Induktionskonstanten fu, deren 
zahlenmaBige Festlegung den gleichzeitigen Ubergang zur rationalen Schreibweise der Gleichun- 
gen ermoglichte. (Die Einfiihrung der beiden GréBen ¢) und wy kann man auch als einen 
Schritt zur Einfiihrung eines Fiinfer-Systems ansehen, das aus drei mechanischen, einer elek- 
trischen und einer magnetischen Grundgr6Be besteht. In den folgenden Ausfiihrungen wollen 
wir jedoch die Frage der Fiinfer-Systeme ausklammern und uns lediglich mit den oben ge- 
nannten Beziehungen zwischen den Dreier- und Vierer-Systemen befassen. Ebenso soll die 
bereits von anderer Seite [3] angeschnittene Frage ausgeklammert werden, ob nicht schon - 
in den CGS-Systemen die elektrischen und magnetischen GréBen als VierergréBen angesehen 
wurden.) . 

b) Durch den BeschluB, als elektrische Einheiten der Vierer-Systeme die ,,praktischen 
Einheiten“ zu verwenden, also z. B. als Spannungseinheit das Volt, als Stromstarkeeinheit das 
Ampere [4]. Dies gilt inzwischen fiir alle bisher bekannt gewordenen Vierer-Systeme (GIORGI- 
System, MiE-System und Internationales Einheitensystem). 


Die Einheit Maxwell und das Internationale Einheitensystem 


Wir betrachten nochmals die Einheitengleichung (6). Auf der rechten Seite stehen zwei 
Einheiten, die beide urspriinglich Einheiten eines Dreier-Systems waren und heute als Vierer- 
system-Einheiten aufgefaBt werden kénnen. Die eine dieser beiden Einheiten (s) ist als solche 
in die Vierer-Systeme tibernommen, die andere (V) durch BeschluB zur Vierersystem-Einheit 
ernannt worden. 

Natiirlich kann man sich heute auf den Standpunkt stellen, daB Gl. (6) dann nicht mehr 
giiltig sein konne, wenn man die Einheiten V und s als Vierersystem-Einheiten auffasse, da ja 
die Einheit Maxwell (wie die Einheiten Gau8 und Oersted) Dreier-Einheit geblieben sei. 
Aber Gl. (6) ist heute nochimmer richtig, wenn man V und s als Dreier-Einheiten auffaBt. Es 
ist also nicht mdglich, die Entspricht-Beziehung anstelle (Gl. 6) als einzig mégliche Beziehung 
zu bezeichnen, sondern es mu8 heute daneben auch Gl. (6) als giiltig angesehen werden. DaB 
dieses Nebeneinander nicht nur unschén ist, sondern auch bedenklich verwirrt, liegt auf der 
Hand. Die Unklarkeiten kommen aber nicht daher, daB heute ganz tiberwiegend Gl. (6) und 
nicht die Entspricht-Beziehung benutzt wird, sondern daher, da8 man seinerzeit die elektri- 
schen Einheiten des Quadrant-Systems gleichzeitig als Dreier- und Vierereinheiten hat gelten 
lassen, aber diese Freiheit den magnetischen Einheiten des CGS-Systems nicht zugebilligt 
hat. Will man heute diese Unklarheiten beseitigen, so sind grundsatzlich zwei Wege méglich: 
Entweder man hebt die Giiltigkeit der Namen Volt, Ampere usw. fiir die praktischen Einheiten 
als Dreier-Einheiten auf und erklart damit Gl. (6) als nunmehr ungiiltig, oder aber man erkennt 
auch die Einheiten Maxwell, Gau8 und Oersted als Vierer-Einheiten an, holt also auf diese 
Weise das seinerzeit Versiumte nach. Es bleibt dann Gl. (6) giiltig, wahrend die Entspricht- 
Beziehung tiberfltissig wird. DaB sich die praktisch rechnende Elektrotechnik fiir den zweiten 
Weg entscheiden wird, der ja dem bisher durchweg geiibten Brauch entspricht und der durch 
den AEF schon vor Jahren festgelegt wurde [5], dariiber diirfte wohl kaum ein Zweifel be- 
stehen. Damit ware wahrscheinlich auch das wirksamste Mittel gefunden, die Dreier-Systeme 
als solche fiir die Elektrotechnik endgiiltig in der Versenkung verschwinden zu lassen. 
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Die Einheiten Gau8B und Oersted 


Fiir die Einheit GauB des Dreier-Systems gilt sinngemaB das zur Einheit Maxwell 
Gesagte. Im Dreier-System ist 


i Ggee1 Om) SCM, (7) 


Auf der rechten Seite stehen nur Einheiten, die inzwischen im Internationalen Einheiten- 
System unmittelbar (s), mittelbar (cm = 10-% m) oder durch Beschlu8 (V) Heimat gefunden 
haben. Auch hier sollte man nicht mehr die Einheit GauB als Dreiersystem-Einheit auffassen 
und das Entspricht-Zeichen setzen, sondern es bei der heute durchweg benutzten Gl. (7) 
belassen und sowohl V als auch G als Vierer-Einheit betrachten. 
Bei der Einheit Oersted scheinen die Verhialtnisse nicht ganz so einfach zu liegen, wenn 
man von der Einheitenbeziehung 
G 
Oe = i | (8) 
ausgeht. Hier erscheint in der Einheitenbeziehung die Induktionskonstante fm, deren Ein- 
fiihrung, wie bereits erwahnt, bei dem Ubergang vom Dreier- zum Vierer-System eine Rolle 
gespielt hatte. In der Tat beeinfluBt auch die Induktionskonstante die Beziehung zwischen 
den Feldstarkeeinheiten der Vierer- und Dreier-Systeme durch das Auftauchen des Faktors 4 a 
(Verhaltnis der rationalen zur entsprechenden nichtrationalen Einheit). Da die Stromstarke- 
einheit des elektromagnetischen CGS-Systems zehnmal gréBer als die Einheit Ampere des 
Quadrant-Systems ist, hat man nach Ubernahme der Einheit Ampere in die Vierer-Systeme 
A A 
We gO) Na SOy ee 10° (9) 


m 
gesetzt. Auch hier stehen, auBer auf der linken Seite von Gl. (g), Einheiten, die in die Vierer- 
Systeme iibernommen worden sind. Holt man dies fiir die Einheit Oe nach, so kann Gl. (9) 
unbedenklich weiter benutzt werden, und die schwer verwendbare Entspricht-Beziehung wird 
hier ebenfalls iiberfliissig. 
Weitere Entwicklung 


Wenn die Einheiten Maxwell, GauB und Oersted als Vierer-Einheiten iibernommen sind 
und damit die weitere Giiltigkeit der vorstehenden Gl. (6), (7) und (9) gewahrleistet ist, dann 
kénnte man schlieBlich — wenn man Wert darauf legt, vdllig klare Verhdaltnisse zu schaffen 
und jede Doppeldeutigkeit zu vermeiden — in Aussicht nehmen, die Namen sowohl der prak- 
tischen Einheiten (Volt, Ampere usw.) als auch der Einheiten Maxwell, GauB und Oersted von 
den Dreier-Einheiten zu trennen und sie nur noch fiir die Vierer-Einheiten zu verwenden. 
Bis zur Verwirklichung dieser Lésung kann man folgendermafen verfahren: Will man zum 
Ausdruck bringen, da8 in den Gln. (6), (7) und (8) Dreier-Einheiten gemeint sind, mége man 
die magnetischen und elektrischen Einheiten mit dem Index III versehen, im anderen Fall mit 
dem Index IV. Nach den bisherigen Erfahrungen kann man wohl sagen, daf sich in den mei- 
sten Fallen eine derartige Index-Kennzeichnung ertibrigt. MiBverstandnisse sind nicht méglich. 
Eine Beziehung zwischen Dreier-Einheiten und Vierer-Einheiten sollte man nicht tiber ein Ent- 
spricht-Zeichen herzustellen versuchen, sondern eine Gleichung mit einem System-Verbin- 
dungsfaktor (k) benutzen, z. B. fiir die Verbindung zwischen dem Dreier-GauB und der Vierer- 
Einheit Tesla (1 T = 1 Vsm_2) die Gleichung 


Gy = &-10%T (10) 
mit k= Giy/Gyy =V10’Am/Vs = 4%] - 
Literatur 
[1] Srixve, U.: Arch. Elektrotechnik Bd. 406 (1952) S. 254, insbesondere auch FuBn. 5.— [2] STILLE, U.: 
Messen und Rechnen in der Physik, S. 213 u. 214. Braunschweig 1955. —— [3] Wattor, J.: Arch. Elektro- 
technik Bd. 40 (1952) S. 256; J. FiscHer, ebenda, S. 257. — [4] wie [2], S. 216. — [5] DIN 1325, Ausgabe 


1939; DIN 1301, Ausgabe 1955; DIN 1339, Ausgabe 1958. 
Professor Dr.-Ing. E. FLEGLER, Aachen, Schinkelstr. 


Technische Hochschule, Rogowski-Institut fiir Elektrotechnik 


1* 


: : ; Archiv fur 
4 J. Fiscuer: Die Einheiten Maxwell, GauB, Oersted Hicktrotechaile 


Die Einheiten Maxwell, GauB, Oersted 
Von 


J. FiscHER, Karlsruhe 


(Eingegangen am 12. Dezember 1959) 


Dem praktischen Elektrotechniker und Physiker sind die folgenden Einheiten fur den 
magnetischen InduktionsfluB @, fiir die magnetische Induktion (Flu8dichte) B und fir die 
magnetische Feldstarke H gelaufig: 


[O) =a Maxwell ="10Volts, 


[B] = 1 GauB = 1 Maxwell/cm* = 10°°——_, (1) 


Ampe 
(Fl) =a \Ocerstedi 2s ae 


AN = GM 


Wir wollen fiir sie im Folgenden die Kurzzeichen M, G, Oe benutzen und fiir Volt V, fiir Am- 
pere A schreiben. 

Man findet diese Einheitenbeziehungen von jeher in Lehr- und Nachschlagebtichern der 
Elektrotechnik, in Tabellenwerken, in Normblattern! und in Lehrbiichern der Physik®. Der 
Praktiker benutzt diese Einheiten manchmal weniger deshalb, weil sich mit ihnen handliche 
Zahlenwerte ergeben, was ja nicht in jedem denkbaren Fall zutreffen kann, sondern mehr aus 
dem Grunde, weil er dann Zahlenwertgleichungen und Zahlenangaben, die er in Nachschlage- 
werken, in Tabellenwerken und in der wissenschaftlichen Literatur von Jahrzehnten findet, 
in seinen Zahlenrechnungen unmittelbar und ohne Umformungen und Einheiten-Umrechnun- 
gen einsetzen kann. Selbst wenn der Praktiker in dem Wunsch, einer Empfehlung zu folgen 
oder einer Vorschrift zu gehorchen, kiinftig diese Einheiten nicht mehr selbst benutzen wiirde, 
so diirfte er trotzdem die Kenntnis der Beziehungen (1) nicht verlieren, denn damit wiirde er 
sich selbst die Literatur verschlieBen. 

Um die Herkunft und die Berechtigung der Beziehungen (1) wird sich der praktische 
Elektrotechniker und Physiker im allgemeinen wenig kiimmern, er wird jedoch mit einem 
Blick auf die rechten Seiten der Gleichungen feststellen, da8 ihm auch von der praktischen 
meBtechnischen Seite her keine Schwierigkeiten erwachsen, denn er findet seine Strommesser 
immer in Ampere, seine Spannungsmesser immer in Volt, seine FluBmesser in Maxwell aus- 
geteilt. Den Zahlungsfaktor 10/4 a = 0,796... wird der Praktiker nicht anders hinnehmen, 
als die Umrechnung von Zoll in Zentimeter, die ebenso zu seinen Berufsfertigkeiten gehdrt, 
wie die Kenntnis der englischen Sprache. 

Gegen die Beziehungen (1) sind in den letzten Jahren Einwendungen laut geworden. Da 
diese geeignet sind, Unsicherheit hervorzurufen, stellen wir hier zwei der wichtigsten Gesichts- 
punkte noch einmal zusammen. 


I, 


Der erste Einwand beruft sich auf Beschliisse zwischenstaatlicher Kérperschaften, die 
1900 und 1930 gefaBt worden sind‘; nach diesen gilt 


M/cm? = G = Oe = — \a- eas (2) 


s cm cm 
1 Zum Beispiel DIN 1339, Einheiten magnetischer GréBen, Ausgabe 1958; DIN 1301, Einheiten, Kurz- 
zeichen, Ausgabe 1955. 
* Becker, R.: Theorie der Elektrizitat 1. Band, 16. Aufl., neu bearbeitet von F. SAUTER. Stuttgart 1957. 
* Zum Beispiel ist in der Elektrotechnik recht haufig die Einheit Kilogau8 bequem, wenn es sich um 
magnetisch weiches oder magnetisch hartes Eisen handelt. 
* Zusammenfassende Darstellung in: U. StiLLz, Messen und Rechnen in der Physik. Braunschweig 1955. 


XLV. Band 


Heft 1 — 1060 J. FiscHer: Die Einheiten Maxwell, GauB, Oersted 5 


Hiernach sind die Einheiten G und Oe gleich, und alle drei Einheiten sind ersichtlich mechani- 
sche GroBen deswegen, weil cm, g, s, drei Grundeinheiten der Mechanik sind. (Die Einheiten 
(2) wurden und werden absolute elektromagnetische Einheiten genannt.) Nun sind aber V 
und A anerkannt als Einheiten eines Einheitensystemes, das nicht drei, sondern vier unab- 
hangige Grundeinheiten hat, namlich heute m, kg, s, A®. Die Gleichungen (1) kénnten, so 
wird gesagt, deswegen nicht bestehen, weil M, G, Oe nach (2) einem Einheitensystem mit drei, 
dagegen V und A einem Einheitensystem mit vier unabhangigen Grundeinheiten angehoren. 
G und 10-* V s/cm? zum Beispiel entsprachen zwar einander, kénnten aber deswegen nicht 
- einander gleich sein. 


Der praktisch messende Physiker und Elektrotechniker (und ihm soll ja durch Festsetzung 
von Einheiten geholfen werden!) wird allerdings vielleicht gegen diese Auffassung Bedenken 
haben. Er wird geltend machen, daB er mit demselben MeBinstrument jederzeit dieselbe 
GréBe gemessen habe, ganz unabhangig davon, ob zu zwei verschiedenen Zeiten zwei ver- 
schiedene Definitionen giiltig waren, oder ob in den beiden Fallen dieselbe Definition gegolten 
hat. In der Tat war es auch stets, wenn eine Anderung einer Definition beschlossen wurde, 
eines der wichtigsten Anliegen der BeschluBfassenden, den Betrag der Einheit nicht zu andern, 
mit anderen Worten, den MeBtechniker die Anderung der Definition nicht merken zu lassen*. 


Man kann gegen diesen ersten Einwand heute das Folgende sagen: 


1. Man hat heute nicht mehr den Wunsch, die elektrischen und magnetischen Erscheinun- 
gen und Vorgange mittels mechanischer GroBen zu beschreiben, also mittels solcher GréBen, 
deren Grunddimensionen nur die Masse, die Zeit und die Lange sind. Nur zu diesen mechani- 
schen GroéBen passen die in (2) genannten mechanischen Einheiten. 


2. Die in (2) genannten Einheiten entstehen, wenn man annimmt, da8 man die Induk- 
tionskonstante fy (die Permeabilitat des leeren Raumes) gleich der reinen (unbenannten) Zahl 
eins setzen kénne. Die Meinung, da8 diese willkiirliche Verfiigung niitzlich ware, wird heute 
iiberwiegend nicht mehr aufrecht erhalten, es lassen sich im Gegenteil Griinde dafiir angeben, 
daB sie nicht einmal zulassig sei. Ein wesentlicher Fortschritt in der Theorie, namlich eine 
auBerordentliche Vereinfachung und zugleich Vervollstandigung, ist von J. WALLOT erzielt 
worden®. Nach dieser Theorie sind cm, g,s und my die Grundeinheiten der absoluten elektro- 
magnetischen Einheiten, und an Stelle von (2) gilt 


C= 1 Lo ae —| Mo - dyn We 
s | 4ncm cm 4% cm?’ 


z pe (3) 
1 g 1 dyn 
Oe ) ay 
s | 4up~)>.cm cm |/ 47M 


daher G = py: Oe. Diese Beziehungen sind aber mit den Beziehungen (1) identisch’. 


m? kg 
. s8 A 

54 Den Namen Ampere haben hintereinander drei Stromstarkeeinheiten erhalten: die Quadranteinheit 
(ein Potenzenprodukt aus den drei mechanischen Einheiten cm, g, s), das selbstandige Silbervoltameter- 
Ampere und das Ampere des heutigen ,,Systemes der Internationalen Einheiten“, und innerhalb einer gewissen 
Me8genauigkeit sind die drei verschieden definierten Ampere gleich groB. 


5 Daher zum Beispiel V = 


8 Wattot, J.: Handbuch der Physik, herausgegeben von H. GercER und K. ScHEEL, Band II, Kapitel I, 
GréBengleichungen, Einheiten und Dimensionen, Berlin 1926; 2. Auflage, Leipzig 1957. 


? Man rechnet die Gleichungen (3) in die Gleichungen (1) leicht um, indem man beachtet 
cme Se —= dyn 105’ Ve Avs/eni-, 


4% VS 4% dyn 


ly = = 
tae 10? Acm 100. «= A2 


Beides sind exakte Beziehungen, keine Naherungen. 
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Der zweite Einwand (der mit dem ersten zusammenhangt), stellt sich fiir den Praktiker 
so dar: Die in (1) genannte Beziehung 


A 
1 Oe =— — (4) 


ist unzutreffend; ist sie in der Praxis benutzt worden, so beruht das auf einer irrtiimlichen Auf- 
fassung; zutreffend ist die Beziehung 


A 
1 Oe = 10—. (5) 
cm 


Dieser Einwand wiegt wesentlich schwerer, als der erste, denn man kann ihm nicht, wie 
bei I., mit dem Hinweis begegnen, da8 die Theorie der Einheiten seit dem Zeitpunkt Alterer 
BeschluBfassungen weiter fortgeschritten ist, vielmehr handelt es sich hier um einen Unter- 
schied in der absoluten GroéBe der Einheit Oersted um den Zahlenfaktor 42. Es lohnt sich, 
iiber Zahlenfaktoren von der GréB8enordnung einer Zehnerpotenz nachzudenken. 

Als im Jahr 1930 beschlossen wurde, der absoluten elektromagnetischen Einheit der 
magnetischen Feldstarke den Namen ,,Oersted“ zu geben, war mit gr6Bter Wahrscheinlichkeit 
die Einheit des (iiblichen und verbreiteten) nicht rationalen elektromagnetischen CGS-Syste- 
mes gemeint. Wir wollen dies annehmen. Es waren zwar damals auch schon die Begriffe: 
GréBengleichungen und Zahlenwertgleichungen bekannt, jedoch ist das Oersted nicht durch 
eine Gleichung festgelegt worden, weder durch (4), noch durch (5), und auf GréBengleichungen 
irgendwelcher Form ist nicht Bezug genommen worden. Es ist daher hier sicher zweckmaBig, 
zuerst (a) die magnetische Feldstarke, die 1 Oersted genannt werden soll, durch eine erdachte 
Versuchsanordnung mit bestimmten Versuchsbedingungen festzulegen, hierauf (b) die 
Beziehungen zu den GréBendefinitionen und GréBengleichungen herzustellen. Dabei wird es 
sich zeigen, daB dieser zweite Vorgang unerlaBlich ist, denn die Einheitengleichungen eines 
Systemes hangen nicht allein von den zugehorigen Zahlenwertgleichungen ab, sondern auch 
von den GréBengleichungen. 

a) Experimentelle Definition des Oersted durch eine Normalanordnung 

Besonders einfach ist der gerade, unendlich lange fadenformige Stromleiter von vernach- 
lassigbar kleinem Durchmesser; er selbst und der ganze Raum, in den er eingebettet ist, 
soll dieselbe Permeabilitat haben. Wird er von einem Strom J durchflossen, so besteht in 
einem radialen Abstand 7 eine magnetische Feldstarke H. Da es sich um die Feldstarkeeinheit 
des nichtrationalen elektromagnetischen CGS-Systemes handeln soll, miissen die nichtrational 
geschriebenen Zahlenwertgleichungen dieses Systemes herangezogen werden, hier also die 
Gleichung 


Peay 
ae (6) 
Der Zahlenwert {H} der magnetischen Feldstarke ist auf die Einheit 1 Oersted bezogen, der 
Zahlenwert {J} der Stromstarke auf 1 elektromagnetische CGS-Einheit gleich 10 Ampere, 
der Zahlenwert {7} des senkrechten Abstandes von der Achse auf 1 cm. (Zahlenwertgleichun- 
gen sind sinnlos, wenn nicht angegeben wird, auf welche Einheiten sie bezogen sind.) Wird 


daher der Stromleiter von 10 A durchflossen, so ist im Abstand von 2 cm der Zahlenwert der 
magnetischen Feldstarke 


WER its (7) 
und dieses ist also der Zahlenwert der magnetischen Feldstarke, bezogen auf die Einheit 
1 Oersted. — Natiirlich kann man auch andere Normalanordnungen wahlen, etwa eine lange 


diinne Zylinderspule, oder einen kreisférmigen Stromleiter; es miissen dann in entsprechender 
Weise gewisse Zahlenwerte (Versuchsbedingungen) gegeben werden’. Alle diese Méglichkeiten 


* Zum Beispiel: Die sehr lange diinne zylindrische Spule hat auf 1 m axialer Lange w = 1000 Windungen, 
die vom Strom 10 A durchflossen werden. Dann besteht im Inneren der Spule die magnetische Feldstarke 
4% Oersted. — Hier gilt namlich die Zahlenwertgleichung {H} = 4 a w {I}/{1} des nichtrationalen elektro- 
magnetischen CGS-Systemes; die Lange ist auf die Einheit 1cm, die Stromstirke ist auf die nichtrationale 
elektromagnetische CSG-Einheit gleich 10 A bezogen, daher ist {H} = 4a. 
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wollen wir ,,experimentelle Definitionen‘‘ der magnetischen Feldstarke nennen. Das allen 
Gemeinsame ist die Form 


p {H} {ds} = 40 {1} (8) 
des Durchflutungsgesetzes, dieses als Zahlenwertgleichung geschrieben, bezogen auf nicht- 
rationale elektromagnetische CGS-Einheiten. : 


Fiir diese experimentelle Definition ist also keine Voraussetzung dariiber gemacht, wie 
etwa ihrerseits die elektromagnetische CGS-Einheit der Stromstarke definiert ist. 


b) Definition des Oersted mit Hilfe von GréBengleichungen 


Hierfiir ist die WALLotTsche ,, Verkniipfungsbeziehung‘ von entscheidendem Nutzen. Sie 
sagt Folgendes aus®: Besteht zwischen GréBen A, G,...G, eine allgemeine GréBengleichung 
von der Form . 


A ZG Go ci GO, (9) 
in welcher Z der (exakte) Zahlenfaktor ist (zum Beispiel Z = 1/2 in W = mv®/2), sind 
CAV NG A ee Gea (10) 
die zugehérige Zahlenwertgleichung und 
[A] = ¢- [GJ*- [Ge]*-.-.+ [Go] ey 
die zugehérige Einheitengleichung, so gilt stets die Verkniipfungsbeziehung 
Lee (12) 


wie man leicht durch Einsetzen allgemein bestatigt. Man sieht, da in der Tat die Einheiten- 
gleichungen (¢) nicht allein von den Zahlenwertgleichungen (z) abhangen und umgekehrt, 
sondern daB die GréBengleichungen (Z) an dem Zusammenhang maBgebend beteiligt sind*®. — 
Einheiten heiBen koharent, wenn in der sie verkniipfenden Einheitengleichung ¢ = 1 ist. Sind 
die Einheiten koharent, so haben die GréBengleichung und die ihr entsprechende Zahlenwert- 
gleichung dieselbe Form, denn es ist dann € = 1 = Z/z oder z = Z. 
In Gl. (6) 

2 {1} 

(6) 
ist also der Zahlenfaktor z = 2. Alles Weitere folgt mit Riicksicht auf die Verknipfungs- 
gleichung (12) zwanglaufig: 

«) Definiert man eine GréBe H’ durch die GréBengleichung 


{H} = 


Ox AL 
Speen (13) 


if 
so folgt aus dieser und der Zahlenwertgleichung (6) die Einheitengleichung 
[J] 10 A 
Yo om 


[H’] 1 Oe’. (14) 
[In der GréBengleichung (13) ist Z = 2, in der Zahlenwertgleichung (6) ist z = 2, in der 
Einheitengleichung (14) ist daher € = Z/z =1.] Nach der GréBengleichung hat in der ge- 
wahlten Anordnung die magnetische Feldstaérke im Aufpunkt die GréBe 

2eLORE A 


H' = = to0-—== 1-0e 
2 Cin cm 


9 Wattot, J.: Arch. f. Elektrotechn. Bd. 40 (1952) S. 325—331; GroBengleichungen, Einheiten und 
Dimensionen, Paragraph 35. Leipzig 1957. 

10 Beispiel: Gehért zu der GrdBengleichung W = mv?/2 die Zahlenwertgleichung {W} = {m} {v}?/2, 
so ist notwendig die Einheitengleichung [W] = 1 [m] [/]?/[£}?; aus Z = 2 etOleG notwendig € = 1. Hatte 
man zu derselben Zahlenwertgleichung die GréBengleichung W’ = m v? gesetzt, so ware aus Z = 1, 2 = 2 ge- 
kommen ¢ = 2, die Einheitengleichung hieBe [W’] = 2 [m] [7]?/[#]?. Das triviale Beispiel zeigt den Nutzen 
der Verkniipfungsgleichung (12) fiir die Darstellung. Hierzu vergleiche auch FuBn. 11. 
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8) Definiert man eine andere GréBe H durch die GréBengleichung 


I 
H=—_, (15) 
270 7 
so folgt aus dieser und der Zahlenwertgleichung (6) die Einheitengleichung 
1- [1] 1 10 A 
[7] = — — = — — = 16 
eu 4x [7] Wege (enon Or ( ) 


[In der GréBengleichung (15) ist Z = 1/2 z, in der Zahlenwertgleichung (6) ist z = 2, in der 
Einheitengleichung (16) ist daher € = Z/z = 1/4.] Nach der GréBengleichung (15) hat in der 
gewahlten Normalanordnung die magnetische Feldstaérke im Aufpunkt die GroBe 

10_A 1SeLORA 


[ig a = = AOE 


2U7-2CM 4m7 Cm 


Das magnetische Feld als physikalischer Zustand ist natiirlich in beiden Fallen a), £) 
dasselbe. Stromstarke und Abstand sind ja in beiden Fallen dieselben. Dasselbe magnetische 
Feld wird nur durch zwei verschieden definierte GréBen H’ und H gekennzeichnet. Wir 
vergleichen die beiden Gr6Ben und die beiden Einheiten: 


H =anH, Oe =42 0e, (17) 
daher , S 


Wie wir mit (6) vorausgesetzt haben, ist der Zahlenwert in beiden Fallen gleich, und in dem 
besonderen Falle der Normalanordnung ist er also {H}= 1. Die GréBe H’ ist um genau 
denselben Zahlenfaktor gréBer, als die GréBe H, um welchen die Einheit Oe’ gréBer ist, als 
die Einheit Oe. 

Wir nennen H’ und Oe’ die nichtrational definierte Gré8e und Einheit, H und Oe die 
rational definierte GréBe und Einheit. Die Beziehungen, die zwischen rational und nicht- 
rational definierten Grd8en und Einheiten bestehen, lassen sich hiernach ersichtlich mit 
Hilfe der WatLotschen Verkniipfungsgleichung besonders einfach darstellen ©. 

Fiir die Praxis ist es sicherlich unzweckmaBig, zwei verschiedene Gr6B8en und zwei ver- 
schiedene Einheiten nebeneinander zu fithren. Da der erwahnte Beschlu8 iiber die Namen- 
gebung ,,Oersted“ auf GrdBendefinitionen und GroBengleichungen, wie schon erwahnt, 
keinen Bezug genommen hat, kann hier eine zusatzliche Entscheidung auf gar keine 
Weise umgangen werden. Entscheidet man sich fiir die nichtrationale Definition H’ der 
Gr6Be, so hat man die Einheitengleichung (14), entscheidet man sich fiir die rationale De- 
finition H der GroBe, so hat man die Einheitengleichung (16). Mit der nichtrational definierten 
GroBe A’ gilt das Durchflutungsgesetz in der Form | 


PH ds = Aa, (19) 
mit der rational definierten GroBe H gilt es in der Form 
f Joh Cky — oe (20) 


11 Ein anderes Beispiel fiir eine experimentelle Definition und fiir die Anwendung der WaLLotschen Ver- 
knipfungsgleichung ist das CouLtompsche Kraftgesetz der Elektrostatik und die Definition der elektro- 
statischen CGS-Einheit [Q], der elektrischen Ladung: Fiir zwei gleich groBe elektrische Ladungen Q schreiben 
wir das CouLomssche Gesetz als GréSengleichung 


ON Gn Ge le, MISC) iss thang. (B 1) 
Ausgeschrieben : 
{QO}? [OQ = 4% {e} {7}? {F}. le] W? [F) . (B 2) 
Die elektrostatische CGS-Einheit [Q], ist definiert durch die Bedingung, daB {Q} = 1 seifiir {e} = 1, {v7} = 1, 
{Ff} = 1; [€] = &, [yr] = cm, [F] = dyn. Setzt man diese Bedingungen in (B 2) ein, so ergibt sich die Ein- 
heitengleichung 
[013 =4me,cm@%dyn, also C=4qn. (B 3) 
Dividiert man (B 2) durch (B 3), so erhalt man * 
Omer AT 8 ty aleO. a sae, (B 4) 


in der Tat ist {Q} = Z/& = 1, wie fiir die experimentelle Definition verlangt war. 
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J. FiscHER: Die Einheiten Maxwell, GauB, Oersted 9g 
Die Entscheidung fallt klar zugunsten der rational definierten magnetischen Feldstarke H und 
daher der Einheit Oersted nach (16) und (4) aus, wenn man folgende Gesichtspunkte bedenkt: 
1. GroBengleichungen sind insofern die ideale Darstellungsform fiir physikalische Sach- 
verhalte, als sie von der Willkiir der Einheiten frei sind. Sie haben die einfachste Form, wenn 
sie rational geschrieben werden [man vergleiche zum Beispiel Gl. (19) mit (20).] 

2. Hinsichtlich der Einheiten geht die Entwicklung immer deutlicher auf die Benutzung 
der Einheiten hin, die aus den Grundeinheiten kg, m, s, A koharent abgeleitet werden. Dieses 
ist ein rationales System, die Zahlenwertgleichungen werden rational geschrieben. Wegen 
der Koharenz (¢ = 1) ist Z= 2, mit anderen Worten: Die Zahlenwertgleichungen dieses 
Systemes haben dieselbe Form, wie die rational geschriebenen GréBengleichungen. Verwendet 
man dagegen nicht rational definierte GréBen, so geht dieser groBe Vorteil verloren. 

Beide Griinde sprechen dafiir, die rational definierten GroBen und daher die rational 
geschriebenen GréBengleichungen zu benutzen!®. Ihnen gegeniiber ist der Nachteil gering, 
daB in der Einheitengleichung (16) der Zahlenfaktor 4 7 auftritt. 

Wir fiihren noch einen dritten Grund an: Oben war festgestellt worden, daB iiber die Be- 
schluBfassung der Namengebung ,,Oersted“ hinaus eine zusatzliche Entscheidung unumgang- 
lich ist, wenn man Eindeutigkeit herstellen will. Diese zusatzliche Entscheidung darf man 
darin sehen, daB in Technik und Physik bisher praktisch ausnahmslos die Gleichung 


4% cm 
anerkannt und benutzt worden ist. Die Praxis hat also die notwendige zusatzliche Entschei- 
dung schon gefallt, sie hat sich fiir diese Einheit und daher fiir die rational definierte magne- 
tische Feldstarke H schon entschieden. Sie benutzt diese und keine andere Einheitenglei- 
chung, und zwar ganzlich unabhangig davon, wie die Stromstarkeeinheit Ampere definiert ist. 

Liegt die Einheit Oe fest, so sind auch die Einheiten G = wy: Oe und M = cm?- G ge- 
geben. Es ist fp = 4%: 10-° V-s/A cm. 

III. 

Indem wir verallgemeinern, k6nnen wir zu der Frage der rationalen und nichtrationalen 
Definitionen von GréBen und Einheiten Folgendes sagen: 

a) Die Zahlenwertgleichungen der CGS-Systeme sind von jeher ganz tiberwiegend nicht- 
rational geschrieben worden (z = 1 in einigen Gleichungen). Dies ist eine historische Tat- 
sache, die man, weil sie historisch ist, nicht mehr andern kann. 

b) Die beste Schreibweise der GroéBengleichungen der Elektrizitatslehre ist die rationale, 
aus den beiden oben angefiihrten Griinden. 

c) Verwendet man zusammen miteinander rational geschriebene GroBengleichungen und 
nichtrational geschriebene Zahlenwertgleichungen, so ist notwendig in einigen Einheiten- 
gleichungen ¢ =: 1, einige Einheiten koharieren nicht. 

d) Erhebt man jedoch die Kohiarenz aller Einheiten zum Postulat (¢ = 1 ausnahmslos), so 
mu8 man einerseits bei den nichtrationalen CGS-Systemen nichtrational definierte GroBen, 
also nichtrational geschriebene GroBengleichungen einfiihren, andererseits bei dem heute 
eingefithrten praktischen System rational definierte GroBen, also rational geschriebene GroBen- 
gleichungen. 

e) Dies ist eine in physikalischer Hinsicht ganz unnotige Doppelspur; man darf urteilen, 
daB sie dem Sinn der GréBengleichungen widerspricht. Dem gegeniiber ist die Inkoharenz 
einiger Einheiten bei dem Vorgehen nach c) offenbar der kleinere Nachteil. Sie wird ja auch 
sonst in der Praxis hingenommen. Es ist wahrscheinlich einfacher, in einigen Einheiten- 
gleichungen einen Zahlenfaktor ¢ + 1 zu haben, als fiir eine Anzahl von GréBen zwei vonein- 
ander verschiedene Definitionen kennen und auch wissen zu miissen, wann diese, und wann 
jene Definition benutzt werden muB. 

12 So auch in dem Normblatt DIN 1313, Schreibweise physikalischer Gleichungen, Ausgabe 1931. 


Professor Dr.-Ing. Jou. FiscHER, Technische Hochschule Karlsruhe, 
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Beitrag zur Kenntnis der Ionenbewegung in Halbleitern 
Von 


FRANZ OLLENDORFF 


(Mitteilung aus dem Forschungsinstitut fiir Elektronik am Technion, 
Israel Institute of Technology, Haifa) 


Mit 7 Textabbildungen 


(Eingegangen am 25. August 1959) 


Inhaltsiibersicht: 1. Ziel der Arbeit. — 2. Die Grundgleichungen. — 3. Lésung des Potentialproblemes. _ 
4. Entwicklung des Halbleiter-Potentiales in der Umgebung des Kanales.— 5. Das quasistatische Potential. — 
6. Die Bremskraft. — 7. Zahlenbeispiel und Kritik. — 8. Zusammenfassung. 


1. Ziel der Arbeit 


Die moderne Halbleiterphysik beruht auf den Ausbreitungsgesetzen der Wahrscheinlich- 
keitswellen, welche die Elektrizitatstrager bei deren Bewegung durch das sozusagen ein- 
gepragte, dreifach-periodische Potentialfeld des jeweils vorliegenden Kristallgitters begleiten. 
Auf Grund dieses einfachen Mechanismus kénnen — im Grenzfalle verschwindend niedriger 
Absoluttemperatur — immer genau jene Elektrizitaétstrager ihren Wohnraum hemmungslos 
durchwandern, deren individuelle Gesamtenergie 7 einem der ,,erlaubten‘‘ Energiebander 
An angehort. Um daher den bei hinreichend hoher Temperatur tatsachlich auftretenden 
Bewegungswiderstand genetisch zu verstehen, muB man die thermischen Vibrationen des 
Kristallgitters in Rechnung stellen, welche gema8 ihrer wesentlich EINSTEIN-DEBYEschen 
Quantenstatistik dem vorher statischen Potentialfeld 6rtlich und zeitlich stochastisch schwan- 
kende Stdérungen aufzwingen; die urspriinglich ebenen Wahrscheinlichkeitswellen werden 
nunmehr teilweise gestreut und reflektiert. Ungeachtet der erfolgreichen Anwendung dieser 
Theorie auf das Leitfahigkeitsproblem verschweigt sie doch sowohl die dynamische Ein- 
wirkung je eines wandernden Ladungstragers auf seine beweglichen Nachbarn wie auch die 
von ihm geweckte Polarisation der Gitteratome. Da nun diese Erscheinungen in das physi- 
kalische Verhalten der Halbleiter entscheidend eingreifen, sollte man versuchen, ihnen rechne- 
risch nachzugehen; die vorliegende Arbeit will einen Beitrag zur Lésung dieser Frage liefern. 


2. Die Grundgleichungen 


Es hegt gewiB nahe, die gestellte Aufgabe als Vielkérper-Problem der Wellenmechanik 
anzusetzen; doch reichen meines Wissens die bisher bekannten mathematischen Hilfsmittel 
zur Loésung der entsprechenden SCHROEDINGER-Gleichung durchaus nicht hin. Unter dem 
Zwang dieser Sachlage werden wir daher die phanomenologischen, MAxwELLschen Feldglei- 
chungen als adaquate Beschreibung des elektromagnetischen Geschehens im Halbleiter anzu- 
sehen haben. Der Einfachheit halber ersetzen wir den tatsdchlich vorgegebenen. Kristall 
durch einen homogenen und isotropen Kérper der elektrischen Leitfahigkeit x und der Di- 
elektrizitatskonstanten ¢, den wir tiberdies mit der allerdings nur fiktiven Eigenschaft ver- 
schwindender Permeabilitat ~ ausstatten; in ihm orientieren wir uns an Hand des Bezugs- 
systemes der Zylinderkoordinaten z (Achse), 7 (Radialdistanz) und « (Azimut). Nun richten 
wir unsere Aufmerksamkeit auf ein individuelles Ion der invarianten Ladung g, welches 
wir vorerst als materiellen Punkt im Sinne der klassischen Mechanik auffassen. Indem wir 
sowohl seinen elektrischen wie seinen magnetischen Spin geflissentlich auBer Betracht lassen, 
begniigen wir uns mit der Angabe seiner tragen Masse m, welche, im Einklang mit der voraus- 
gesetzten Isotropie des halbleitenden Korpers, als Skalar aufzufassen ist; doch mag ihre 
GréBe von der jeweiligen lonengeschwindigkeit v abhangen, welche ihrerseits wahrend 
der Dauer der beabsichtigten Kontrolle als gleichférmig gelte. Die dann gewiB geradlinige 


XLV. Band 


Heft 1 — 1960 F. OLLENDORFF: Ionenbewegung in Halbleitern 11 


Ionenbahn mit der z-Achse des Bezugssystemes identifizierend, diirfen wir den Augenblick 
der Ionenpassage durch den Ursprung als Nullpunkt der laufenden Zeit ¢ wahlen, so daB die 
Gleichung 
z= ut; 7, =0[—o <t<o] (1) 

den Ort (z,;7,;) des kontrollierten Ions oder, genauer gesagt, seines Schwerpunktes fiir alle 
Zeiten schildert. 

Um uns von allen kinematischen Schwierigkeiten bei der Bewegung des Ions durch den 
Halbleiter zu befreien, statten wir diesen mit dem virtuellen ,,Kanal‘‘ 


Ta Vy et OOnee 2-22 OO (2) 

aus, welcher von Materie frei sei und also phanomenologisch durch die Angaben 
Xp = 0 ; & = 1 (3) 

beschrieben wird; seine Permeabilitat werde wiederum gleich Null gesetzt. Uberdies ersetzen 
wir die bisher in (z,; 7;) konzentriert gedachte Punktladung voriibergehend durch eine Linien- 
ladung, welche langs der Strecke 
2 pare Ole 2 Om OS 0 (4) 
der Achse in der gleichférmigen Dichte 


stetig verteilt sei, auBerhalb dieser Strecke jedoch identisch verschwinde; sie wird somit durch 
das FourtErsche Integral 


lo) a6 lo) 
i ilz’ , in 1 —-il(z— 
A(z; t) == (ieee dl | aoe dz =- faite Lie 99 di] 
—oo 46 —oo 
ex . 
ee Up eae i(z — vt) dl (6) 
na 16 


in ein kontinuierliches Spektrum harmonischer Wellen je der Wellenzahl / zerlegt, 
welche samtlich mit der einheitlichen Phasengeschwindigkeit v in Richtung der positiven 
z-Achse fortschreiten. Eine in der Ebene z = const im Kanal aufgespannte Kreisflache wird 
daher von dem , eingepragten“ Konvektionsstrom 
oO 
Hest) = ves) = 42 | cos (2 — v4) al 3 (7) 
rt 16 

durchflossen; gesucht wird das von ihm im Raume 7 > 7% des sonst unbegrenzt gedachten 
Halbleiters geweckte elektrische Feld. 

Nachdem die Permeabilitat sowohl in dem hier benutzten, fiktiven Substitut des wahren 
Halbleiters wie innerhalb des Kanales y < % definitionsgema8 verschwindet, muB sich iiberall 
dort auch die magnetische Induktion annullieren. Daher erweist sich der Vektor E der 
elektrischen Feldstarke als wirbelfrei, sodaB er als negativer Gradient des skalaren Po- 
tentials  dargestellt werden kann 

E = — grad (8) 

Durch A, die sogenannte Dielektrizitatskonstante des leeren Raumes bezeichnend, 
tritt nun innerhalb des Kanales gemaB (3) neben dem eingepragten Konvektionsstrom (7) 
ein Maxwellscher Verschiebungsstrom der Dichte 


jy = dy ey = — Ay grad S ; 41<iN (9) 


auf, wahrend sich ebendort die Leitungsstromdichte annulliert 


1, = My L = — x grad yp = ; TeV, (10) 
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a 


Im Halbleiter dagegen gesellt sich zu der Verschiebungsstromdichte 


pee eee te prade” yt (11) 
ot ot 
die Leitungsstromdichte 
j, =“xE=—kegrady; 41>" (12) 
Hier wie dort geniigt die Gesamtstromdichte 
f= hy ti» (13) 
der Kontinuitatsgleichung 
div 7x==07 (14) 
Die Annahme einer je gleichférmigen Leitfahigkeit enthalt implizite innerhalb des Kanales 
einerseits, im Halbleiter andererseits — nicht jedoch an deren gemeinsamer Grenze! — die 


Voraussetzung der Quasineutralitat, an welcher wir weiterhin strikte festhalten. Innerhalb 
der genannten Gebiete befriedigt daher das Potential m die LAPLAcEsche Gleichung 
ey ey 1 Op 1 Op 
a ne Les < 1 
Nis emir ers Fae ers ome (15) 
so daB die kinematische Forderung (14) sogar fiir den Verschiebungsstrom und den Leitungs- 


strom einzeln erfiillt ist. 
Wir erganzen diese allgemeinen Angaben durch die Randbedingungen, welche der 


Funktion  aufzuerlegen sind: 
1. Die Radialkomponente D, der elektrischen Induktion 
DIAG ot age (16) 
im Kanal ist der dort linienhaft in der Dichte 2 gemaB (6) langs der Achse verteilten Ladung 
genetisch verkniipft: 


2 20 


lim [| D,rdx=—A,lim 1 da Qian (17) 
(Late) 10 if 
a=0 G= 


2. An der Grenze r = 7% bleibt das Potential  stetig 
lim @ (% — Ar) = lim g (% + Ar) ; Ar>o. (18) 
Aro Ar—o 
3. Die Radialkomponente 7 der Gesamtstromdichte durchflieBt den Grenzzylinder 7 = 1% 
stetig 
lm 


Ar—o 


ay : ein a 
A a= lity) | Ae se merge a ; 1 
Cova yen aa O° & at Or |r, 4 Ar (19) 


4. Mit wachsendem Radialabstand 7 des Aufpunktes von der Systemachse soll das Potential 
gegen Null konvergieren 


lhimg =o. (20) 
3. Lésung des Potentialproblemes 
Aus Symmetriegriinden hangt das Potential » nicht vom Azimut « ab, so daB sich die 
LapLacesche Gleichung (15) auf 
prmer a e a (21) 


reduziert. 
Unter dem Primarpotential m® verstehen wir jenes Potential, welches aus gm durch 
den virtuellen ProzeB 
E=> & = 1; %—> X= 0 (24) 


hervorgeht. Da hierbei der gesamte Raum 7 > 0 frei von Materie wird, verwandeln sich (18) 
und (19) in Identitaten. Durch die Symbole H® und H® beziehentlich die HANKELschen 
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Zylinderfunktionen m-ter Ordnung der Arten 1 und 2 bezeichnend, geniigen wir nun der 
LapraceEschen Gleichung (21) unter der Randbedingung (20) nach Wahl einer vorerst willkiir- 
lichen Amplitudenfunktion a = a(l) durch das Integral 
gf?) — 4 | a(l) i H(i 17) cosl (2 — v1) dl, (23) 
14 


dessen Konvergenz ebenso wie jene der spater zu entwickelnden Integrale weiterhin voraus- 
gesetzt wird. Die Radialkomponente E) der elektrischen Primarfeldstarke berechnet sich 
also zu 


ap) q re ‘ 
Hp = =e 4 few [— H(i 1r)] cost (eg — vt) 1 dl. (24) 
Mit Hilfe der Formel 
lim2aar[— H®(ilr)] 1 = 4 (25) 


hiefert demnach (6) im Verein mit (17) die Relation 


sin 61 1 sin dl 
44ya=—- 3 al) = Ta a (26) 
so daB das Primarpotential durch 
(p) q e sin Ol. Vay, 
Pe eA, morrmeaees (cl 7) cost (zg —v#) dl (27) 


dargestellt wird. 


Zu unserer urspriinglichen Aufgabe zuriickkehrend, erganzen wir zunachst die HANKELschen 
Zylinderfunktionen jeweils der m-ten Ordnung durch die BEssELschen Funktionen 


= — [HD + HP) (28) 
und die NEUMANNSchen Funktionen 
N, = 7H? — He). (29) 


Nun rufen wir vier vorerst unbekannte Amplitudenfunktionen c(/); s(l); C(l); S(2) zu 
Hilfe, mittels derer wir die gesuchte Lésung in folgender Form ansetzen: 


1. Wir begeben uns in das Innere des Kanales 7 < 7. Da wir dort der Randbedingung (17) 
bereits durch das Primarpotential geniigt haben, kann sich das gesuchte Gesamtpotential » 
von gy) nur um ein Sekunddrpotential  unterscheiden, welches in 0 <7 < 7% aus- 
nahmslos stetig ist; da es iiberdies die LApLAcEsche Gleichung (21) befriedigen muB, werden 
wir somit auf die Integraldarstellung 


og) = rik i Jo(t 2 7) [c(l) cos 1 (zg — vt) 4+ s(d) sind (z—v4)] dl ; r<n (30) 


gefiihrt. Aus (27) und (30) resultiert die Radialkomponente 


E,=— a SA: i [— Hy (il r)] cosl (zg —vt)ldl+ 


oO 


co 


= sop eRe) [c(l) cos 2 (zg — vt) + s(l) sind (z — vt)] ld] (31) 


oO 
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der elektrischen Feldstarke; sie erregt gem4B (g), (10) und (13) die radiale Gesamtstromdichte 


oE in 
jp = gt Ef PEE OG I] +6 HIN el] sind @—wt) Pal + 


[oe] 


oe i J,(il7) s(l) cosd(z —vt) dl. (32) 
2. Wir gehen in den Halbleiter 7 > 7% tiber, léschen dort das Primarpotential gp!) und 
wahlen im Hinblick auf (20) die Lésung 


[o,e) 


e- 3 | t HYG 17) (C(l) cosl (g —vt) + Sl) sind (z —v2)] dl ; 1; (33) 


der LapLAceschen Gleichung (21) als Ansatz des Gesamtpotentiales; es erzeugt die Radial- 


komponente 
CO 


ee — ns Pa i [— HY (il 7)] (C(l) cost (z —vt) + S(l) sind (z—vt)]2 dl, (34) 


welche ihrerseits gemaB (11), (12) und (13) die Radialkomponente 


[o,0) 


; : 
i,=(« + AcS)E =e | | HP (il) {xS(l) + 1 Aye v C(l)} sind (2 — vt) Ldl+ 
ei 
ate al [— AP(ilr)] {x Cl) —1 Aye v S(l)} cos 1 (zg —v dt) idl, (35) 


der Gesamtstromdichte zur Folge hat. 

Um nun die vier unbekannten Amplitudenfunktionen zu bestimmen, bedienen wir uns der 
Stetigkeitsbedingungen (18) und (19). Das gemeinsame Argument (7 / 7) der unterschiedlichen 
Zylinderfunktionen der Kiirze halber unterdriickend, finden wir dann die vier Gleichungen 

sin oo 


i HY + Joc) =i HP CQ), (36) 
Jy s(t) = 7 H® Si), (37) 

se ae H® +i J,c()} =— Ao wre S(l) + ec], (38) 
—i J, s() = — H® as) e S(l)}. (39) 


In ihnen diirfen wir den Grenziibergang 6-> 0 zur Punktladung des bewegten Ions aus- 
fiihren und erhalten mit 


DexJo HY AG v aay ' ne {¢ Jo Hy? — J, Hy} (40) 
die Angaben 
si) _ Jo HY) Sra Fs 7 en) 
Cay so | tae (42) 
s() = 5) = (43) 
a(t) = "2 (cay. (44) 


Es wird sich zeigen, daB es in den spater zu berechnenden Integralen hauptsachlich auf 
den Bereich 


i 
oS (45) 
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der Wellenzahl / ankommt. Mit Riicksicht auf das dort jeweils maBgebliche Verhalten der 
unterschiedlichen Zylinderfunktionen vereinfachen sich daher die vorstehenden Ausdriicke 
in ausreichender Genauigkeit zu 


w 7 Ao é v)? + x? 2 (lL Age v)? + x? 
D= Hy Ne Sl, Ao Ke (46) 
: 1 Agu 
S(/) ( Ae v)2 + 2’ (47) 
(1 Agv)? € 
CW) Ta area (48) 
No 1 Ayu x HD (7d 1 
s( ) —_ (1 Aye ve” Ay (2 %) , (49) 


€ (e— 1) (Av)? +22 . ; 
el) = are 1 HEH). (50) 
In der gleichen Genauigkeit liefert daher (33) fiir das Potential im Halbleiter die Integral- 
darstellung 


co 
qv [ lAjevcosl (g—vt) + xsin/ (z— vi?) 
2G: (1 Ay ev)? + x? 


te) 


i Hil) Lal, (51) 


mit welcher wir uns fortan zu beschaftigen haben. Insbesondere ergibt sich im Grenzfalle 
#%—> 0 [Isolator) das relativ zur wandernden Ebene z =z; elektrostatische Feld der 
Punktladung g 


(oe) 


+ ye ee a /; 
Peele Age cosd (2 —v1)7 Hy’(lr) dl = 


q t 
4m Aye y'(g— vi? + 7 : (52) 


“x—->O0O 


wahrend der allerdings nur fiktive Grenzfall e ~ o des unpolarisierbaren Leiters in 
limg = rae t)i Hil) lal 
imp = 2° | sin (g—vt)t H(ilr) lc 


&—0 
oO 


qu & Mn aU lee al gee ( 

== 4mx O2V (zy iP +P 4%%[(z —vt)? + PB 53 

die Strémung schildert, welche einer mit der gleichfoérmigen Geschwindigkeit v langs der 
z-Achse fortschreitenden Doppelquelle vom Betrage (gv) ihres Momentes entspringt. 


4. Entwicklung des Halbleiter-Potentiales in der Umgebung des Kanales 


Um die Struktur der Strémung im Halbleiter dem physikalischen Verstandnis zu er- 
schlieBen, kehren wir zu Gl. (51) zuriick und versuchen, aus ihr eine explizite Beschreibung 


des Potentiales y in der Umgebung des Kanales zu gewinnen. Zu diesem Zwecke begeben 
wir uns zunachst in die komplexe 


orc a (54) 
Ebene und schreiben 
qu 
eae [®; + Dy] (55) 
mit 
1 o etl(z—vt) ; 1 al eb (s— v2) 
DP, =e | Tee ae ame BOG l 7) Lol) = a | Aven H (pr) pdp , (56) 
sowle | 
1 rs e—tl(z—vt) ; 1 estan 
Py = = | jAgsumzin Ho (t!7) ldl =| Pare H®(p vr) p dp. (57) 


re} oO 
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Von nun an haben wir zwei Falle analytisch wesentlich verschiedenen Charakters zu 
behandeln: 


1. Es werde 


(zg —vt) >0 (58) 
vorausgesetzt, so daB also der Aufpunkt vor der wandernden Ebene z — 2; liegt. Die Relation 
i HOG lr) = —1 HP (—tl7) (59) 
ausnutzend, erhalten wir dann mittels der konformen Abbildung 
pak +i =—p, (60) 
in der ~’-Ebene fiir ®; die Darstellung 
0,=—-+ [ay Pn bab = 3 | Fee PO nv ae 6) 


> 
= 
= 
@ 
= 
aS 
a 
neuer Weg 2 
a 
= 
Sa 
= 
x. 
s 
= neuer Weg 
7 
= 
a 
5 
p-Ebene 
Bild 1. Zur Auswertung des Integrales (61). Bild 2. Zur Auswertung des Integrales (63). 


Da zufolge (58) der letztangegebene Integrand mit |f’|—> oo in der Halbebene x > 0 
verschwindet, diirfen wir den urspriinglich langs der negativ-imaginaren Achse der #’-Ebene 
verlaufenden Integrationsweg entsprechend Bild 1 in deren positiv-reelle Achse verlegen 
und finden 


fo} 


1 ek’ (z ma th t) (2) ; , . 
© ie rig reir (eet) Rade >. (62) 


Die namliche Uberlegung liefert in der p-Ebene entsprechend Bild 2 fiir ®,; die Dar- 
stellung 


co 


1 e—k(z—v!) a) 
ones RAS (kr) kdk. (63) 

Die in (62) und (63) beziehentlich benutzten Zeichen k’ und k der Integrationsvariabeln 
mit dem einheitlichen Symbol y vertauschend, erhalten wir somit aus (55) im Hinblick auf (28) 
fiir das Potential m des Halbleiters die Darstellung 


qv e—y (e—v?) 
eS an pAsraie di ae Ute Oe rites (64) 


oO 
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In ihr ersetzen wir die Strecken (z — v #) und ¢ beziehentlich durch ihre, in der ,,natiirlichen 
Langeneinheit“ 


Ayev 
M=—— (65) 
gemessenen, dimensionsfreien Werte 
z—vt ¥ 
eines oe. Oe aay Oe 
schreiben 
EVE nT (67) 
und finden 
q fete 
> 4 ise NE Poi Joel) Par (68) 
Mittels der Potenzreihe 
(sat 
fore) Bee 
y\ 2 (6) 
Nel) = yey! ve 


der BEssELschen Funktion nullter Ordnung geht daher (68) in die Summe 
se) ps a\” 2¥ q2” i —{T 


P40 AyeM deol (v2 \2} atey} P44 


Oo 


iiber. Um das hier verbleibende, bestimmte Integral 


eons 
fe) = | Fzrar, (72) 


18) 


auf tabulierte Funktionen zuriickzufiihren, substituieren wir anstelle von I" die Veranderliche 


u=C(f +1). (72) 
Mittels der in 


CO 


Ba) =| ~~ du (73) 


x 


gegebenen Definition des Exponential-Integrals von negativ-reellem Argumente wird 


dann 
= FN Nar oe 7 ao ee 1+¢feEi(—) 
w=] fe au =e | u dy] = SHES EHCD (74) 
g 


entsprechend Bild 3a. Aus (73) und (74) entnimmt man die Relationen 


df 1 Clie it 
a OP a eae (75) 


durch deren Restitution in (70) die fiir hinreichend kleine numerische Radialabstande @ und 
fiir gentigend groBe numerische Abstande ¢ > 0 des Aufpunktes von der wandernden Ebene 
z = 2; brauchbare Entwicklung 


a ee 2 S | | eee 
= see Ola — pti} (76) 


entsteht. 
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2. Wir gehen zum Falle 
(g —vt) <0 (77) 


iiber, so daB nunmehr der Aufpunkt hinter der wandernden Ebene z = 2; liegt. GemaB (56) 
verschwindet daher jetzt der Integrand von ®, fir |p|> oo in der Halbebene k > 0. Nun 
konstruieren wir entsprechend Bild 4 um den Pol 


x 1 


bo =F + th ; No tg Re =o, , (78) 
als Zentrum den in/ > 0 gelegenen Halbkreis 
p=fyt Ape”; TZOLZO; Apy < Ro (79) 


urspringlicher Weg 


Bild 3. a Die Funktion /(£) nach (74); b Die Funktion g(¢’) = g(—€) nach (92). Bild 4. Zur Berechnung des Integrales (80). 


und erhalten, unter Berufung auf den CAucHyschen Hauptsatz der Funktionentheorie, 


Ry — Apo 
1 ek (z — vt?) 


@, = ; HV(kr)k dk + 


2 k Ayev— 


1 @ * 
ere, i eto? HD (Lk + Ap €'°\7) (ty + Apy &*) i a8 + 
1 e ek (z— vt) 
+t f GSP amen sar. (80) 
Ro + APo 


Zur Berechnung von ®;; gema8 (57) ziehen wir die konforme Abbildung (60) heran 


1 qs e—P (21) 1 ts eb’ (2 —vt) 
= (1) , , , 
On= =| Saag hOon sap =+[ —" noe npvay. x) 


In der p’-Ebene zeichnen wir entsprechend Bild 5 den Halbkreis 


bp =ppt+4he*; —axS0' So (82) 
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und erhalten, durch abermalige Anwendung des CaucHyschen Satzes, 


kyo— APy 


1 ek’ (2 — vt) 2) / 27 D , 
os eee (k r) Rk dk’ + 


k' Ay ev—x 


ie) 
oO 


es HOCK ae Abo oy 7) (Ro ae Ap, li dwy ae 


1 ek’ (2 — 012) ve Da s¥ 


urspringlicher Weg 


Bild 5. Zur Berechnung des Integrales (83). 


Fihren wir jetzt den Grenziibergang 
Ap) > 0 (84) 
aus, vertauschen, wie oben, k und k’ mit dem einheitlichen Symbol y und bedienen uns der 
Relation (29), so folgt aus (55), (80) und (83) 


kyo—Afo seed 
qu 3 f ey (7— vt) ev 05) 
eae alee A See ee Ay an ee 
4 Ap —+-O } 0 A ? 0 
itis ky + Abo 
x 
- wipe (¢— vt?) 


be 
: xi( 2) ; (g—vt)<0; rN - (85) 


Ayev es Ayev 
Wir setzen hier, mit Benutzung der natitirlichen Langeneinheit M nach (65), in Erginzung 
von (66) 


aaron ee 0 (86) 
und erhalten mit Riicksicht auf (67) und (85) die Darstellung 
pA 
: 0 prrar eT ery rar = 
~ 4% AyeM ae | r—1 Jo(b L’) ai ; Tea Jo(e L) +e A) (87) 
: o Apy 
Saas 


2* 


> ; Archiv fiir 
D() F. OLLENDORFF: Ionenbewegung in Halbleitern Blalirotechnik 
es 


des Halbleiterpotentiales; sie verwandelt sich mit Hilfe der Potenzreihe (69) in 


q Seles ee ae Saf rciae | oF a —' N,(p) 
7 yadyeM| < l)? \2] a" yy : Lae eee 


Apy->0 i Ady 
aad Py 
(88) 
Zur Berechnung der Funktion 
Apo 
Miele cae 
ss ta 
Boies: ite | SP ae | ian (89) 
Apo Sie Up wa 
14+ a 
substituieren wir anstelle von J’ die Veranderliche 
Wir en) ! (90) 
und erhalten auf Grund der Definition 
cae pte n etw 
11 (90) it | ~ du’ 4 | ; aa : SuceiO: 4 x 0 (91) 
Eo mes oe é Us 
des Exponential-Integrals von positiv-reellem Argument x aus (89) 
1 ble iC" 
gc) ==, (92) 
entsprechend Bild 3b. Aus (91) und (92) resultieren die Beziehungen 
dg 1 st Areas ree rn. 
ae ea 80) 5 ace = ga t pat 80); bah? AQ (93) 


durch deren Einsetzen in (88) die fiir hinreichend kleine numerische Radialabstande @ und 
fiir geniigend groBe numerische Abstande ¢’ > 0 des Aufpunktes von der wandernden Ebene 
z = 2, brauchbare Entwicklung 
q “ oe? [2 1 : Pe 
0= chau 6) “lat me PeBlC ci mae rhea om NGO (94) 
des Halbleiterpotentiales entsteht. 


5. Das quasistatische Potential 


Im Ausdruck (94) des fiir den Bereich (z — v#) < 0 zustandigen Halbleiterpotentiales 
erscheint neben einem nach ganzen, positiven Potenzen des radialen, numerischen Aufpunkt- 
abstandes @ fortschreitenden Anteil die singulare Komponente 


q —l’ / 
Pt =Tahem =e NolQ);  U >o, (95) 


die schon allein fiir sich der LapLacEschen Gleichung geniigt. Welche physikalische Bedeutung 
kommt ihr zu? 

Wir behaupten, daB das Potential y,, ungeachtet seiner gemaB (86) impliziten Abhangig- 
keit von der laufenden Zeit ¢ ein wesentlich elektrostatisches Feld definiert und deshalb als 
quasistatisches Potential bezeichnet werden darf. 

Zum Beweise dieses Satzes bilden wir zunachst die achsial weisende Komponente der 
quasistatischen elektrischen Feldstérke und erhalten mit Riicksicht auf (86) 


OP st nn OPst q 52 
ee dae Ter'aicmnmne ne UC (96) 


Sie erregt die ebenso gerichteten Komponenten 


OE s z Ayev OE ¢ Zz g z 
2.5 2 hing ’ = 
ot Mat pga pee Ne) (97) 


Ist, 2(V) Aye 
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der Verschiebungsstromdichte und 


; q Po 
Ist,2(L) — * Dae sz Sn Aye?” e~ N,(@) (98) 


der Leitungsstromdichte, welche einander zu verschwindender Achsialkomponente der Ge- 
samtstromdichte erganzen 


ey rai Ist, tv) ee Ist, (2) == 0% (99) 


Ahnlich finden wir aus (95) fiir die Radialkomponente der quasistatischen elektrischen 
Feldstarke unter Beachtung der Definition (66) den Ausdruck 


TG LOS, Sa REIN te q 


= =~ aos 
SEY Or M 20 47 Aye M2 @ N,(@) (100) 
so daB 
A OE s+ Ave v OE t, q op 
Ist,0v) = Ay € at ae M SE ee tae Ae ae : N,(@) (101) 


die von ihr hervorgerufene Komponente der Verschiebungsstromdichte und 


; q ; 
1st, (L) aaa Dee) aa oS press Mt = N,(e) , (102) 


jene der Leitungsstromdichte in Richtung wachsenden numerischen Radialabstandes o miBt; 
wiederum kompensieren einander beide Stromanteile zu verschwindender Radialkomponente 
der Gesamtstromdichte 


lay xy Ist,r(V) aie Ist, r(L) cere i (103) 


Im Lichte der Aussagen (99) und (103) offenbart sich also das Feld des Potentiales p, als 
stromfrei; damit ist der gewiinschte Beweis seines quasistatischen Charakters erbracht. 
Wir erganzen und vertiefen diese Uberlegungen an Hand der Analyse der flachenhaft 
langs der Kanalwandung 7 = ~% ausgebreiteten elektrischen Ladung. 
Zunachst richten wir unsere Aufmerksamkeit auf jenen Teil der Kanalwandung, welcher 
vor dem bewegten Ion liegt. Die dort herrschende radiale Induktion D, ergibt sich aus (76) 
zu 


2 


op Aye op q o |2 1 * y 
D,=— Aye eee oF 20 xs al an a +10)| — Sean bls 0 = =aF° (104) 
Nun sei die Ionengeschwindigkeit v als so groB vorausgesetzt, daB der numerische Kanal- 
halbmesser @, der Ungleichung 


Y% % 


Lo FM Awe (105) 


gehorcht. Bei AusschluB der nachsten Umgebung des Ions kommt dann gemaB (104) an 
dem betrachteten Teil der Kanalwandung auch nur eine sehr schwache Radialinduktion 
zustande, so daB wir die von ihr gebundene Flachenladung vernachlassigen diirfen. 

Wenden wir uns jetzt zu dem hinter dem fliegenden Ion gelegenen Teil der Kanalwandung, 
so haben wir uns an die Potentialdarstellung (94) zu halten und finden fiir die Radialinduktion 
die Summe 


Aye ay Aye d(p—Hs) Ane Ope E Cae ’ | 
Smee NT a de ee e+ pet ett) — + zi 
q i Y 
eno Nil0)s © 0 = O= 5, (106) 


Nun gilt fiir hinreichend kleine @ die Entwicklung 


Nie) = te (107) 
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Bei abermaligem Ausschlu8 der unmittelbaren Umgebung des Ions fallt daher — unter 
der Voraussetzung (105) — der absolute Betrag des zweiten Postens der Summe (106) sehr 


groB gegen jenen des ersten aus, so daB wir (106) hinreichend genau durch 
q py Ear WS 
Die re are : N (0): 5 at 0 = 0 (108) 


approximieren diirfen, Demnach enthalt jede Langeneinheit der Kanalwandung ¢’ > 0 die 


Ladung 


(as 
= 2D, =20%M OD, =— ze » (109) 

so daB auf der gesamten Kanalwandung die Ladung 
Go No dedea AE Ge (110) 


gebunden ist. 

Man kénnte versucht sein, die Gegenladung (+ q) in der Ebene € = ¢’ = 0 zu suchen, die 
ja den Existenzbereich des quasistatischen Potentiales y, begrenzt und daher eine Un- 
stetigkeitsflache eben dieses Potentiales definiert. Doch ist diese Meinung durchaus urig: 
Ungeachtet der genannten Unstetigkeit des Teilpotentiales @,, bleibt das Gesamt- 
potential m des Halbleiters auch an jener Ebene stetig! Im Lichte dieser Dialektik mit 
Einschlu8 der vorausgesetzten Quasineutralitat des gesamten Halbleiters erweist sich also die 
Ladung (— q) der Kanalwandung als eine Art raéumlich verteilter [,,verschmierter‘] Senke 
fiir den von der Ionenladung (+ q) ausgehenden Induktionsflu8; dabei sind folgende Falle 
zu unterscheiden: 

1. Solange sich das Ion in einem Isolator befindet [x —>0; M— oo], verschwindet ’ 
- gem&B (109); die Feldlinien verlaufen merklich geradlinig langs der im Ion zentrierten, gegen 
eine unendlich groBe Hiillkugel zielenden Strahlen. ' 

2. Bei der Bewegung des Ions im Halbleiter [M endlich] kriimmen sich die Feldlinien 
riickwarts und folgen ihrer Quelle wie der Schweif eines Kometen seinem Kopfe, wobei ihre 
lineare Senkendichte mit wachsendem Abstand vom Ion exponentiell abfallt. Diese phano- 
menologische Aussage also beschreibt in iiberaus anschaulicher Weise die Wirkung des fliegen- 
den Ions auf die seiner Bahn benachbarten, antipolaren Elektrizitatstrager, welche ihm 
schwarmartig nachlaufen! 

3. Falls der Halbleiter zu einem ,,vollstandigen“ Leiter entartet [x oo, M -> 0], konzen- 
triert sich die Gegenladung (— q) auf die nachste Umgebung des Ions, wahrend der Leiter 
feldfrei bleibt ; dasselbe tritt, bei endlichem Werte der Leitfaihigkeit x, im Falle verschwinden- 
der Ionengeschwindigkeit v ein. 


6. Die Bremskraft 


Das elektrische Feld des Halbleiters greift am fliegenden Ion mit einer antiparallel der 
z-Achse wirkenden Bremskraft F an. Welches ist ihre GréBe ? 

Zur Beantwortung dieser fundamentalen Frage bedienen wir uns des Energieprinzipes: 
Um ungeachtet der Wirkung der Bremskraft die nach Voraussetzung gleichférmige Bewegung 
des Ions aufrecht zu erhalten, muB man dieses mit der Antriebskraft 


f= (111) 


in Richtung der z-Achse vorwarts ziehen. Die hierbei aufgewandte Leistung 
Pes (112) 


wird durch Vermittelung der elektrischen Strémung innerhalb des Halbleiters in JouLEsche 
Warme verwandelt. Durch @ und 7, beziehentlich die Werte des Halbleiterpotentiales und 
der Radialkomponente der Gesamtstromdichte am Kanalumfang bezeichnend, finden wir 


nun 
foe} 


P= J Pot, 2%% dz. (113) 
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Zu Gl. (51) zuriickkehrend, berechnen wir zunachst die Radialkomponente der elektri- 
schen Feldstarke 


co 


ye ap qu [ lA,evcosl (z—vt) + xsin] (g—v?) 4). ; 
B= —P=—2 | (1 Aye v)? + x? Hy (il) P dl (114) 
und aus ihr die Radialkomponente der Gesamtstromdichte 
Pee Oe pie oe ect ot) OG Ty) ital (115) 
Vy 0” oat Te Am 1 r 


Auf Grund ihrer Antimetrie in bezug auf das Argument (z—v#) geht in das Produkt- 
integral (113) nur der ebenso antimetrische Anteil 
qu x sin 1 (z— vt?) 
4k (1 Ay ev)? + x? 


ce) 


gp tH alr) dl (116) 
des Potentiales @ leistungsbildend ein. 

Zur Auswertung jenes Integrales machen wir von folgendem Spezialfall des PARSEVALschen 
Theoremes Gebrauch: 

Es seien a, = a,(%) und a, = a,(x) zwei je reelle Funktionen von x, welche gleichzeitig 
die Antimetrie-Eigenschaften 


a, (—%*) =—a,(%) 5 ay (— x) = — (2) (117) 
aufweisen; daher sind sie der FourIER-Darstellungen 
a,(x) = f o,(l)-sinl x dl ; a,(x) = f o,(l) sin 1 x di (118) 


fahig, deren Spektraldichten o,(/) und o,(/) sich aus den Ausgangsfunktionen a,(x) und a(x) 


nach den Vorschriften 
[oe) CO 


g(t) = ={ Gale | sini aes o(l) = i as(% ysindineds, (119) 


oO ce) 


berechnen. Nun bilden wir — unter der Annahme seiner Existenz — das Integral 


if reales An(h) Ota 82 fala) Gs(X) da. (120) 


Ersetzen wir in ihm etwa a, durch die aus (118) zu entnehmende Fourt&rR-Darstellung, so 
entsteht zunachst 
T= 2 [|  f-4,@) 0,0) sin 1 4-di dx (121) 

Onl =o 
und wenn wir hier die Integration nach x ausfiihren und beachten, daB die Integrations- 
variable x’ der Spektraldarstellungen (119) ohne Anderung des Ergebnisses in x umbenannt 
werden darf, 

T=x f o,(l) o,(2) dl. (122) 

l=o 


Bei der Anwendung dieses Satzes auf (113) haben wir gemaB (115) und (116) 


0,0) =— 22 AMG InP: oy = Lage HOI) (123) 


gat (LAg ev)? + x? 
zu setzen und erhalten fiir die Leistung P den Ausdruck 


= DU erate St 870, Oy 3 : 
Wet ann( 22] x | (dentate ’ Ho (il %) [— HyP(tly)) Bd. (124) 


te) 
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Mittels der Substitution 


verwandelt er sich in 


foe) 


1 q \2 ° 4 HYG A) (— A BG Ad) 
P= slam) WN) ne et oy 


oO 


Wir setzen abkiirzend 


NyedA. (127) 


if HDG A) [— A HDG A 
vio) = | a 
Bild 6 zeigt den mittels numerischer Integration ermittelten Verlauf dieser Funktion; 
fiir 0) > o konvergiert sie gegen den Festwert 


aoe ae | i HOG A)[—A HG A] dA=0,5, (128) 


09—>0 


oO 


4 2 3 4 5 
lie Lok ee 
% 7% Y 
Bild 6. Die Funktion y(g9) nach (127). Bild 7. Die {numerische] Bremskraft als Funktion 


der [numerischen] Ionengeschwindigkeit. 


wahrend ihr asymptotisches Verhalten fiir @) > 1 durch die Naherung 


v(eo) & a3 | 6 HOG A) [A HME A] Arad =" (126) 


Qo” 0 
Oo 


beschrieben wird. 
Bringen wir jetzt die Definition des numerischen Kanalhalbmessers fy in die Gestalt 


1 Ay € 
soo ee (130) 


so liefert (112) im Verein mit (111) fiir die gesuchte Bremskraft die Darstellung 


pert 
[F | =F’ = 3H 573% PlO0) (131) 


Bild 7 veranschaulicht ihren Verlauf in Abhangigkeit vom Kehrwert des numerischen 

Kanalhalbmessers 

ii ee Ay € 

00 cae Py % v ’ (132) 
welcher im Lichte dieser Gleichung ein dimensionsfreies MaB der Ionengeschwindig- 
keit liefert. Bei hinreichend kleinen Werten dieser numerischen Geschwindigkeit besteht 
somit Proportionalitat zwischen ihr und der Bremskraft; das Verhaltnis 
. |qv 2 Fae 1 8 x“ 
= lim |—~| = 8 lin =- Ta 
p po | i q| ae Qo” Y(@o) 0,17 lq\ (133) 


% 
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definiert die Beweglichkeit des Ions, welche hiernach nur mit der Leitfaihigkeit x seines 
Wirtsk6rpers, nicht jedoch mit dessen Dielektrizitatskonstante ¢ verkniipft ist. Dieses Ergebnis 
stimmt mit jenem der klassischen Korpuskulartheorie der elektrischen Leitung wohl iiberein: 
Sei m die Dichte der beweglichen Elektrizitatstrager je Raumeinheit, so errechnet sich die 
Leitfahigkeit x aus 


x= |q|nB. (134) 
Der Vergleich von (133) mit (134) fiihrt also auf die Relation 


ee ae 5 (Ot7 4/2 Sf 
on7 m= Aas A oe \ 8 V3 = 0,277 V5 , (135) 
welche anzeigt, daB der Kanalhalbmesser 7 bis auf einen Zahlenfaktor der GréBenordnung 1 
dem mittleren Abstand der genannten Elektrizitatstrager gleicht. 
Die so einfache Beschreibung der Ionenbewegung mittels einer konstanten Beweglichkeit 


versagt jedoch, sobald die numerische Geschwindigkeit des fliegenden Ions sich der GréBe 1 
nahert: Die Bremskraft nimmt dann langsamer als proportional zu, um bei der dem Werte 


id page FL 6 
Qo /kr 1G uN Ua ess (13 ) 

entsprechenden kritischen Geschwindigkeit 

x 1% 
ee YW (137) 
den Hochstwert 

F e668 e : 
Lea cornne poe kta (138) 


zu erreichen. Falls also die Triebkraft F durch ein am Ion angreifendes elektrisches Trieb- 
feld der Starke # geliefert wird 


Fog, (130) 
korrespondiert der Maximalkraft (138) die kritische Feldstarke 
Zz a 
|Z, | = 0,01 Wee (140) 


welche sonach wesentlich unabhangig von der Leitfahigkeit x ist. Wird die kritische Feldstarke 
iiberschritten, so werden die von ihr ergriffenen Ionen hemmungslos beschleunigt: Die kriti- 
sche Feldstarke leitet den inneren Durchbruch des Halbleiters ein! Dieses wohlbekannte 
Phanomen der Halbleiterphysik bedarf demnach zu seinem Verstandnis nicht der wellen- 
mechanischen Theorie des erzwungenen Uberganges der Ionen zwischen zwei durch eine 
,verbotene Zone getrennten Energiebandern, sondern ergibt sich bereits zwanglos aus den 
folgerichtig angewandten Vorstellungen der klassischen Mechanik und Elektrodynamik; doch 
beachte man die kritischen Bemerkungen des folgenden Abschnittes. 


7. Zahlenbeispiel und Kritik 


Wir erlautern die vorstehenden Uberlegungen am Beispiel eines Halbleiters der Daten 
a 1085 ; fetes n= 10% = 
wahrend wir dem bewegten Ion die Ladung 
|¢| = 1,60- 10°! Cb 
gleich dem absoluten Betrage der Elektronenladung zuschreiben. Aus (135) ergibt sich dann 
der Kanalhalbmesser iiber Erwarten groB 


1 10,27 Gee NO. TL, 
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wahrend fiir die kritische Geschwindigkeit gemaB (137) der auffallend niedrige Wert 


10? - 0,277 - 1078 ee 3 mM 
Ferree 16 AT. 1:9 * 10 25 eee 


resultiert. Fiir die Durchbruchsfeldstarke nach (140) finden wir schlieBlich die Angabe 


1,60+ 10719. 4 7-9 - 109 


=. Vv 
|E,,| = 0,04 = A750 gee 


16 - 0,2772 + 10728 
die sowohl mit unseren Vorstellungen wie mit den Beobachtungen einigermaBen vertraglich 
ist. 

Der schwache Punkt der mitgeteilten Zahlenrechnung — und damit der gesamten, vor- 
stehend entwickelten Theorie — liegt ersichtlich in dem jeweils einzusetzenden Wert 7% des 
Kanalhalbmessers oder, tiefergreifend, in der Konzeption des materiefreien Kanales im sonst 
den Raum kontinuierlich erfiillenden Halbleiters; in der Tat wiirde der Grenziibergang zu 
einem ionenleeren Stoff [7 —> 0] gemaB (115) zu dem gewi8 absurden Schlusse eines tiber alle 
Ma8e groBen Kanalhalbmessers fiihren, und hieraus sollte nach (140) sogar eine nur infinite- 
simal schwache Durchschlagsfeldstarke eines solchen Stoffes resultieren! 

Die genannten gedanklichen Schwierigkeiten riihren wesentlich von der Wahl unseres 
kinematischen Modelles her, die indes fiir die mathematische Durchfiihrung der Aufgabe 
notwendig ist und nur wenige Varianten gestattet; zu ihr treten die im Verlaufe der Rechnung 
zusatzlich eingefiihrten, einschrinkenden Annahmen vorwiegend numerischer Natur. Hier 
also liegen die Ansatzpunkte fiir eine erneute, verscharfte Behandlung unserer Aufgabe, um 
die angezeigten Schwachen der phanomenologischen Theorie zu tiberwinden. 


8. Zusammenfassung 


Es wird versucht, die ,,eingepragte“ Bewegung eines Ions in einem homogenen und iso- 
tropen elektrischen Halbleiter gegebener Werte seiner elektrischen Leitfahigkeit x und seiner 
Dielektrizitatskonstanten ¢ phanomenologisch zu beschreiben. Durch systematische Vernach- 
lassigung der magnetischen Induktion gelingt es, die Kenntnis des Feldes auf ein elektrisches 
Skalarpotential y zuriickzufithren, welches iiberall der LApPLAcEschen Gleichung unter ein- 
fachen Stetigkeits- und Randbedingungen gehorcht. Die Lésung wird zunachst in Gestalt 
eines FourrERschen Integrales angegeben. Die hieraus gewonnene, explizite Darstellung des 
Potentiales in der Umgebung der Ionenbahn 1a8t sich anschaulich interpretieren: Das Ion 
schleppt einen Schwarm antipolarer Elektrizitatstrager nach sich, welche dem ,,Kopf‘‘ nach 
Art eines Kometenschweifes folgen. Die hieraus resultierende Bremskraft auf das Ion wird 
energetisch formuliert und unter Benutzung des Parsevalschen Theorems auf das Produkt- 
integral der spektralen Dichten zuriickgefiihrt, welche die Stromdichte und das Potential des 
Halbleiters in FourtERscher Weise schildern. Als Funktion der Ionengeschwindigkeit zeigt 
die Bremskraft nach anfanglich linearem Anstieg ein Maximum, aus welchem die Existenz 
einer kritischen Feldstarke erschlossen wird, die den inneren Durchschlag des Halbleiters 
einleitet. Die Ergebnisse der theoretischen Untersuchung werden an einem Zahlenbeispiel 
erlautert, dessen Aussagen die Kritik herausfordern; die Ansatzpunkte fiir eine verscharfte 
Behandlung des Problemes werden angegeben. 

Es mag schlieBlich gesagt werden, daB die vorliegende, in gewissem Sinne sozusagen ,,alt- 
modische Arbeit die modernen, wellenmechanischen Methoden der Halbleiterphysik nicht 
ersetzen, sondern erganzen will. In der Tat gelten die hier gefundenen Ergebnisse nicht allein 
fiir feste Halbleiter mikrokristallinischer Struktur, sondern schildern in gleicher Weise die 
Ionenbewegung in quasineutralen Elektrolyten oder im quasirieutralen Plasma von elektrischen 
Gasentladungen. 

Fir wertvolle Anregungen bei der endgiiltigen Fassung dieser Arbeit bin ich Herrn 
Prof. J. FISCHER zu Dank verpflichtet. 


Prof. Dr. F. R. OLLENDorFrF, Haifa (Israel) 
Israel Institute of Technology. 
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Elektromagnetische Felder mit axialer Symmetrie 
bei Begrenzungen durch Kegel und Kugel 


Herrn Staatssekretaér Prof. Dr.-Ing. e. h. Dr. phil. FRIEDRICH GLADENBECK 
in Dankbarkeit zu seinem 60. Geburtstag gewidmet 


Von 
H. BucuHoiz, Darmstadt 


Mit 6 Textabbildungen 
(Eingegangen am 16. August 1958) 


Ubersicht. In der vorliegenden Arbeit wird zu Anfang das axialsymmetrische Feld eines hypotheti- 
schen, freistehenden, magnetischen Stromringes bestimmt, der im Inneren eines von zwei vollkommen leitenden 
Kkegeln und zwei ebensolchen Kugeln begrenzten Raums steht, so da ein in sich geschlossener Schwingtopf 
vorliegt. Die geschilderte Art der Anregung des Feldes gewinnt sofort eine physikalische Realitat, wenn der 
magnetische Stromring in eine der Begrenzungsflachen des Raumes hineinriickt. Wird die Intensitat des 
magnetischen Stroms im Ring in Analogie zu dem elektrischen Fall mit Uy in Volt bezeichnet, so entspricht 
diese GréBe der Spannung zwischen den Kanten des sehr schmalen, ringférmigen Spaltes, der dasselbe 4uBere 
Feld erzeugen wiirde. 

Die obige allgemeine Lésung wird dann nach verschiedenen Richtungen hin schrittweise vereinfacht. 
Unter den dabei vorkommenden, offenen Strahlungsfeldern interessiert insbesondere der schlanke, vollkommen 
leitende Kegel mit einer an der Spitze aufsitzenden, vollkommen leitenden Kugel, wobei dann entweder die 
Kugel oder der Kegel den Spalt enthalt. Auch fiir diesen Falle werden die beiden maSgebenden Reihenent- 
wicklungen angegeben. 

In einer anschlieBenden Arbeit wird im letzten Fall dieses Strahlungsfeld auch numerisch berechnet, um 
so zu einem quantitativen Urteil itber den EinfluB der kapazitiven Beschwerung durch die Kugel zu kommen. 


1. Einleitung 


In zahlreichen Aufsatzen, von denen hier nur-einige [1], [2], [3], [4], [5], [6] und [7] er- 
wahnt seien, wird das Strahlungsfeld einer Antenne durch einen ringformigen Spalt von sehr 
kleiner Offnungsbreite erzeugt gedacht. In der Oberflache des Antennenleiters besteht dann am 
Ort des Spaltes ein sehr kraftiges elektrisches Feld in Richtung der Mantelliniedes Antennenleiters 
und senkrecht zu den Konturen des Spaltes. In einer heute mehr als fiinfzehn Jahre zuriick- 
liegenden Arbeit, die in Form einer Broschiire erschien, hat E. METZLER [5] die Eigenschaft von 
rotationssymmetrischen, elektromagnetischen Strémungsfeldern untersucht, die in der Ober- 
flache des Antennenleiters durch eine ,,zonale EMK“ hervorgerufen werden, um eine von dem 
Autor selbst benutzte Bezeichnungsweise zu gebrauchen. SchlieBlich hat der Verfasser der vor- 
liegenden Arbeit in zwei teils noch alteren Aufsatzen [6], [7] tiber axialsymmetrische Strah- 
lungsfelder, die ihre Grenzen an Kegeln, Hohlkegeln oder an Drehparabolen haben, u. a. eine 
Erregung des Feldes durch koaxiale magnetische Stromringe angenommen. La8t man einen 
zunachst freistehenden, zur Achse eines Kegels symmetrischen elektrischen Stromring in die 
vollkommen leitende Oberflache eines Kegels oder einer Kugel mit dem Mittelpunkt in der 
Spitze des Kegels hineinriicken, so wird der Ring kurzgeschlossen, und es hért jede merkliche 
Strahlung nach auBen auf. Handelt es sich aber in dem sonst gleichen Fall um einen magneti- 
schen Stromring, so entsteht ein Feld, das im wesentlichen mit dem elektromagnetischen Feld 
des schmalen Spaltes iibereinstimmt, so daB dieses Feld zum mindesten in seiner Grund- 
struktur auf diese Weise vielleicht am einfachsten hergeleitet werden kann. Auch die An- 
nahme einer zonalen EMK lauft auf die Anregung durch einen Spalt hinaus. 


2.1. Das primiére Strahlungsfeld eines freistehenden magnetischen Stromringes im 
unbegrenzten Raum mittels der skalaren Hilfsfunktion 


Geht die Erregung des Feldes von einem magnetischen Strom mit der Verteilungsdichte 
C™ (7, 9, g’) aus, so berechnet sich bekanntlich (s. z. B. [6]) das von ihm abgestrahlte elektro- 
magnetische Feld iiber den FITZGERALD-Vektor 
et ikyR (oO) = Am 


—1W My SOU As, pea ea [ome ep LY? 


R (om) = V/m? (1) 


vu’ 
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mittels der beiden Formeln: 


E(r, 9, y) = +iawpy: rot (7, 9) , (4 a) 
H(r, 0, ~) = RB- D™(r, 0) = rot rot O™(r, 8) (1b) 
und es werden dann in der Tat die beiden Feldgleichungen 
rotGE=—=+iam:§, (14) 
rot) = —iwes-€ (1B) 


erfullt, wenn Q” die partielle und auBerhalb des Stromringes homogene Differentialgleichung 
rot rot Or, 3) — B.D (r, 3) = 0 (1 A) 
befriedigt. Aus der Gl. (1A) resultiert auch die doppelte Schreibweise der Gl. (1b). In den 


Gl. (1) ist bei einem rein sinusférmigen zeitlichen Verlauf gem4B dem Zeitgesetz exp (— 10 2) 
+14 Mo 
der induktive Widerstand der freien Raumwelle in Ohm/m, 
—10& +x 
die iiberwiegend kapazitive Leitfahigkeit der freien Raumwelle in S/m, 
hy = 0? fy & = (2 0/Ag)° 
die Wellenzahl in m~ im freien Raum, 


fe 7h, Cig ee fee 


die Wellenzahl fiir x 40 mit Im(k) > o0.und & = @- (3 + a) , 
Ao 
die freie Raumwellenlange in m ane hy == 2405 Aes 
dV’ =dr «7 - dd «7 sin} - dg’ 
ein Volumenelement des von dem magnetischen Stromring erfiillten Raumes in m’, 
R?=r+7%—aryr-cosy 

das Quadrat des Abstandes zwischen dem Aufpunkt (7, #, y) und dem Ort (r’, ®, y’) von dV’. 

Dabei ist in dieser letzten Gleichung wegen 
cosy = cos }- cos # + sind. sin # - cos (p’ —q) (2) 


y die dem Winkel w — q’ gegeniiberliegende Seite eines spharischen Dreiecks mit den beiden 
anderen Seiten #, 2%’. 


Da im vorliegenden Fall eine magnetische Strémung nur langs der Kreislinie (7’, 0) in 
immer derselben GréBe und Phase vorhanden sein soll, so kann in Gl. (1) geschrieben werden: 


Cr, Wo’). dV! = U,:7'-sin & -dg’- £3(0, o’) , (3) 
rtA w+ 

lim ff |e, 0 @)|- dr’ or aw =U, [UJ =V. — (3a) 
(4=0,6=0) \’—A #—6 


Liegt gema8 Bild 1 der Stromring zu dem Kugelkoordinatensystem der 7, 3, p (Aufpunkt) 
und 7’, &, yp’ (Quellpunkt) in der Mantelflache des Kegels #’ = const., so daB er gleichzeitig 
senkrecht steht zur z-Achse, so hangt das Feld allein vom Winkel % ab und nicht von gy, und 
es ist demnach mit # = 1, 2, 3 
. (m) E Uy , : , fs / , , eb kk , 
= 40 py DIM, Bf, = Fr’ sin | (B08 9)» dg" (4) 


oO 
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Hierin haben die drei Einheitsvektoren f,mit p = 1, 2, 3, die jeweils in die Richtung zunehmen- 
der Werte von 7, 9 und @ zeigen, aber hinsichtlich ihrer Richtung nur von @ und  abhangen, 
fiir die ganze Integrationsbahn ein und dieselbe Richtung, da sie sich auf den Aufpunkt 
beziehen. Aus Bild 1 ist dann sofort zu ersehen, daB 


i, - 1; = cos (p — 9’) (4a) 
ist. Im Falle # == 1 oder 2 kann geschrieben werden: 
+7 m 


5 et ike Ro) , : rime SaaS ett hye RQ . 
J GB Sea dy’ = [ Ge) <6 @—2a)) dg’ =o, 


—n oO 


Ally—Achse 


Stromring 


Bild 1. Das Bild erlautert die Lage des zunadchst noch freistehenden Stromringes (r’, #’) und die Lage eines beliebigen 

Ringes von Aufpunkten (vy, #) zu dem Kugelkoordinatensystem mit den drei Einheitsvektoren {,, f,,f;, die jeweils in die 

Richtungen zunehmender Werte von 7, # und @ weisen. Die Richtungen dieser drei Einheitsvektoren sind Funktionen 

von #@ und gq, jedoch nicht von ry. Der Raum, in den hinein die zu berechnende Ausstrahlung erfolgt, ist bestimmt durch 

die Angaben: « 29S 8 und y2a. In den Abschnitten 2.1...2.3 kommt noch eine duBere Begrenzung durch die 
Kugel r = b >a hinzu. 


da R(y’) = Rw’ —22) ist. Mithin ist f,-O” =f,-O%) =o. Im Hinblick auf die 
Gl. (4a) ist also 


22 


oy oy : 
Uy» { 00s (p— 9!) dp (5) 


19) 


‘ ' Ry: sin & 
SS 4) [Mo A OM (r, 0) = A iad 


Da nach Gl. (5) 0/€ = 0 ist, wenn die Differentiation an der Funktion Q\” wirkt, so redu- 
ziert sich die Vektordifferentialgl. (1 A) auf die einzige, skalare Differentialgleichung: 

ike Ge 1 Oh al Ort 

ae ~O™) Foy a Seat ny tty ara) 2) 7) , 

(7 QW) +S + (sin 0-5 (Gin B-O0")) + &§- QL" = 0 (1B) 
Sie ist mit der Wellengleichung nicht identisch. 

Das elektromagnetische Feld des magnetischen Stromrings hat nach den Gln. (1a, b) die 

folgenden drei Komponenten: 


| — 

E,y, 8) = 4 2 sind 0) , (6a) 
} a) 

E,(r, 0) =—* =F. ie tee (6b) 

lela 0) = +k. Qi 3 (6c) 


An dieser Darstellung desStromringfeldes durch die Vektorkomponente Q\” des FITZGERALD- 
Vektors 0), die an sich einfach genug ist, ist jedoch auszusetzen, daB die Funktion Q”, wie 
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schon oben bemerkt wurde, nicht der separierbaren Wellengleichung geniigt. Diese Forde- 
rungen befriedigen aber bekanntlich, wenn es sich um Kugelkoordinaten handelt, die beiden 
skalaren, von G. MIE eingefiihrten Hilfsfunktionen, die von ihm selbst mit u(v, 3, gm) und 
u(r, &, p) bezeichnet worden sind. Da es bei Verwendung dieser Funktionen aus dem an- 
gegebenen Grunde leichter méglich ist, andere Partikularlésungen der fiir w und v verbind- 
lichen Wellengleichung anzugeben, als sie in der Darstellung der primaren Erregung bereits 
vorkommen, so mag auch in unserem Falle der Ubergang von Q zu den Funktionen 
oder v vollzogen werden. 

Im axialsymmetrischen Falle stellt sich bei einem magnetischen Stromring der Zu- 
sammenhang zwischen den Komponenten der Felder ©, § und der in diesem Fall zustandigen 
Funktion w in folgender Weise dar: 


Er, 8) =7-|k- u(r, 3) 4 aoe een) ; (7 a4) 
if or” 
1 @(r- u(y, 8)) 
E4(7, 0) mit; Be a oo Bro (7b) 
3 uly, 8 
es 0) = —10 €&° SS . (7 c) 
Da aber hierin (7, 0) der Wellengleichung in der besonderen Form 

1 @(r- uly, 9)) 1 Of: 9 & 2 ee 
y or a .sin® ad (sin is 5) TEs Le ane “) 


geniigt, so laBt sich die Feldkomponente E£, (7, @) statt durch die Gl. (7a,) auch durch die 
Gleichung 
1 ) 


EY, ae y-sind® 0d 


(7 a2) 


(sin 9. duly, 8) 


oo 
berechnen. Der Vergleich von Gl. (7a,) mit Gl. (6a) fithrt iiberdies sofort zu der folgenden 
Beziehung zwischen den Hilfsfunktionen Q”(7, #) und u(r, 8): 


. m G oe 
+ 4 @ My? OM (7, 8) = } ie 


Nun gilt im Hinblick auf Gl. (5) fiir die Funktion exp (+ 7 hy R)/hky R: 
ay Ry fe ap ts | j 1 (ho7) H, 1 (fo7’) 
Cre Sih = + — 24+ 1): $s, a a COs VEG NEN LO 
nO Tee Paes) <r), C00 
falls r’ > vist. Im Falle 7’ <r bekommt die HANKELsche Funktion in Gl. (10) das dann gréBere 
Argument ky 7. 
Die Integration nach g’, auf die es gemaB Gl. (5) nur noch ankommt, ist nun leicht aus- 
fiihrbar, denn es gilt ja bekanntlich fiir die Kugelfunktion erster Art P,(x) die bei beliebigen 
reellen oder komplexen Werten von y unendliche Reihe: 


P, (cos 0+ cos # + sin #- sin 9 - cos (p — ¢’)) 


es! = cee) is —m -m / / 
==.9 Pa ae P,” (cos 0) - P*™ (cos &) - cosm (py —q’). (11a) 


Fiir ein positiv ganzzahliges »y = n bricht die rechtsstehende Reihe hinter dem Glied »y = » 
ab. 

Fihrt man unter Benutzung von Gl. (11a) die in Gl. (5) geforderte Integration aus, so 
erhalt man die Beziehung: : 


(m) _ «Uy, 7’-sin YY Pipi, © rae a 
Om (r, 9) ae 2 Crane Vig i | P=*(cos 8) P¥ (cos #') . 
4, =1 0 0 


(12a) 
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Wegen der Gleichung 


aft ( se 2 + m) a™P_ 1(#) 
atl a (x) 7 ( Ne j= up (x) = (1 xeyeel? _ _ 
—+ r(' nt 7) —= ax 
2 
(ee A 10s 2 es) (11b) 


die bei ganzzahligem m auch fiir die Kugelfunktionzweiter Art Om 1(%) + (—)” gilt, berechnet 


; 2 
sich dann aus Gl. (12a) fiir die Hilfsfunktion u(r, 0, p) selbst, die das primare Strah- 
lungsfeld beschreibt, die Darstellung: 


ler) (7, 0) —— or 1 3 Ro y’ E sin y Ff U,|4 
= Git iho) HG er) 
2h +1 A+ — py, ae 
x = LSS x SS aes Be 1 x VAY A so 
Pa AU ad lige ae (cos §’)  P,(cos #) 


(Ogee = OO)e (12b) 


Geht man in dieser Gleichung in bekannter Weise von der Reihe zum Integral tiber, so erhalt 
man 


yer) (r, 8) oe tag Ug: sin a t 


4aynr' cos at 
S 
(o2r<r’ <oo) (020 <d<a2) 
T (ho 1) + A (hg 7’) ae (+ cos 0”) Lae (— cos #) 
: i < - dt. (12) 
Teg’) (For) { \() + P, (—c0s 6") P__ 4 (cos 6) 
(031 <7r<oo) ; (o26<# sz) ; 


In Gl. (12) ist © ein auf dem Hinweg dicht unterhalb und auf dem Riickweg dicht oberhalb 
der reellen Achse der ¢-Ebene verlaufender Integrationsweg, der von + co kommend, den 
Punkt ¢ = + 3/2 auf der linken Seite in einem Abstand |f—3/2|< 1 umlauft und auf dem oberen 
Ufer zum Punkt + oo zuriickkehrt. Soll die Gl. (12) fiir ein 7 aus dem Bereich oS r<7’ 
gelten, so muB der von der Variablen 7 abhangende Faktor des Integranden die Funktion 
J,(o 7) enthalten, weil nur sie fiir 7 = 0 endlich bleibt, d. h. es gilt in diesem Falle die obere 
Zeile der ersten geschweiften Klammer von Gl. (12). Bei der geschweiften Klammer mit den 
Kugelfunktionen im Inneren ist zu beachten, daf fiir beliebiges komplexes ¢ die Funktion 
P,_1)2 (cos 8) fiir cos @ ~ — 1 unendlich wird, wahrend sie fiir eincos # > + 1 endlich bleibt. 
Dieses verschiedene Verhalten ist die Ursache dafiir, daB auch in der zweiten Klammer des 
Integranden in Gl. (12) nach den angegebenen beiden Fallen unterschieden werden muB. 
Der Vorzeichenwechsel in den beiden Zeilen dieser zweiten Klammer ist eine Folge der Be- 
ziehung: 

P® (— cos 8) = (—)**+™. P% (cos 9) (11¢) 


falls 4, m ganzzahlig sind. 
Da P, *(#) mit —1= % S +1 verschwindet, so darf im tibrigen der Weg auch bis an 
den Punkt ¢ = 0 als dem 4uBersten linken Punkt der Schleife herangefiihrt werden. 


2.2. Das vollstandige Strahlungsfeld eines freischwebenden, magnetischen Stromringes 
im Raum zwischen den Kugeln 7 = a und 7 = b > a und den Kegeln 0? = « und 0 = B>« 


Obwohl dieser allgemeinste Fall raumlicher Begrenzung durch Flachen +7 = const, 
# = const gemaB Bild 2 spater nicht ausfiihrlicher behandelt werden wird, mége er jedoch bei 
der Angabe der Lésung mit beriicksichtigt werden, denn er kann ohne zusatzliche Schwierig- 
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keiten auf dieselbe Weise erledigt werden wie der hauptsichlich interessierende einfachere 
Fall. Da die Wandungen von Kegel und Kugel als vollkommen leitend angesehen werden 
sollen, mu8 die vollstandige Lésung stets 
die Randbedingung erfiillen, daB in der 
Mantelflache der Kegel iiberall E, ver- 
schwindet und in der Oberflache der 
Kugeln dies durchweg bei der Kompo- 
nente EF, der Fall ist. 

Aus Griinden, die erst spater ersicht- 
lich werden, ist es jedoch in Riicksicht 
—— + auf die angestrebte, méglichst allgemeine 
Lésung zweckmaBig, voriibergehend in 
dem Raum zwischen den Kugeln 
a=2r<b und den KegelnazZzdsf 
neben dem Verschiebungsstrom auch 
noch den reinen Leitungsstrom zu _ be- 
riicksichtigen. In den bisherigen Glei- 
chungen braucht dann, sofern es noch 


Stromring 


Bild 2. Dieses Bild stellt ein die z-Achse als Symmetrieachse enthalten- nicht geschehen Site lediglich Ro durch k 

den Schnitt durch einen Schwingtopf dar, der nach Abb. 2 von den oe 

Kugelflachen 7 = a und y = b > a und durch die Kegelflachen & = « oe ; 

und = 6 >« begrenzt wird. Das elektromagnetische Feld, das der und €9 durch & = €9 (3 ai ba | ersetzt 

eingezeichnete, magnetische Stromring im Inneren dieses Schwing- 2 x 0 

topfes erzeugt, 1aBt sich durch die Hilfsfunktionen u(r, #) von Gl. 2 (18) zu werden. Ein solcher Leitungsstrom 
oder 2 (18’’) beschreiben. = 

Im Falle « =o fehlt der rechte Kegel, im Falle 8 = 2 der linke. k6nnte etwa dadurch zustandekommen, - 


daB das Dielektrikum in den Schwing- 
tépfen nicht verlustfrei ist. 
Um den oben angegebenen Forderungen gerecht zu werden, mogen in Gl. (12) zunachst 
die erste der beiden geschweiften Klammern durch den Klammerausdruck 
(r<r’) 
Jk) HENRY’) 
| + A,(t)-J,(kr) + A,(t)- H(z Hh == Kt; Ro 15.0 3G; 0) (13) 
Jer’) -H(b7) 
(r>1) : 
ersetzt werden. Er soll im ganzen Bereich a =r < 6 giiltig sein. Der Lésungsansatz (15a) 
darf also sowohl die Funktion J,(% 7) wie die Funktion H,(k 7) enthalten. Im Hinblick auf die 
Forderung E,(7, ?) = o fiir alle Winkel # aus dem Bereich a = % <P und fir 7 = a oder 
vy = 6 miissen dann wegen Gl. (7a,) und Gl. (13) die beiden Gleichungen gelten: 


ree fe i" 34) Lees eer ! Ria E Ci nt) | 
x=ak Ca Pi 


Ox Ox Ox 
(14a) 
el ae a 
A(t) (2b*- se) + A,(t) -{° Gi np) = — tyr» (ale? 2%) 
a On 4—OF : Ox x=bk ee eee 3 
(14b) 


Wird also zur Abkiirzung 


§ ale? - 74x) ale® «19> ) 
f(a) — | ae za) (15, b) 
*=ak SG=ak 


Ox 


gesetzt, so ist 
jak) WAR) jak) Aa) 


A,(t) = — hOB B) - 7 , 
v OO Tam HOUR jfoR) WB Cae) 
15,02) WOR) eR) HOB) 
A,(t) = + 7,(@ R) - 3} 
v ae Thr hk) HOW by GbR) WD R) 05°) 
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Die geschweifte Klammer (13) kann dann jetzt in der Form geschrieben werden: 
j(ak) Hp (rk) jfak)  hp(ak) 
Kp(ts bs 7, 1/;.0,0) =| Flr k) WGA) — f(D) HDR) 
Tle b) HMB) Tr B) HO" B) 
Gee) 


‘|g(bk) WB R) 


Fir 7 >~/ vertauschen in der letzten Zeile im Element a,; der Determinante 7 und 7’ ihre 
Platze. 

Die Zahlerdeterminante von Gl. (13a) 1aBt sich in den beiden Fallen eines ry => 7 aber 
auch noch in der Form eines Produktes aus zwei zweireihigen Determinanten darstellen. 
Fira <sr</7’ =} b muB8 man zu diesem Zweck die Elemente der dritten Spalte in der drei- 
reihigen Determinante mit J,(7 k) multiplizieren und sodann von der ersten Spalte abziehen. 
Man erhalt dann 


H(bk) ——hy(b R) 
Si eer eum ue of A) HPC) a, 2 y 
LNG OM eg A Del 0) Wee IW B) | AGB) Na ay (Ger Oat sa) 
) 


|7,(b R) (6 A) 


Fir a<r <r = b hat man nach der oben angegebenen Anderung des Elementes a, in 
der Determinante (13a) die Elemente der zweiten Spalte mit H,(y k) zu multiplizieren und 
von der ersten Spalte abzuziehen. Das fiihrt jetzt zu der Beziehung 


Tr kh) HY (v’ k) 
(ak) (ak) 
(1) k) 

k) 


TGs, 7 4-a,.0) ads hae) ee 
1(0 ®) (bd 


Augenscheinlich erfiillen sowohl die Gl. (13a,) als auch die Gl. (13a,) die Randbedingungen 
fiir 7 = a oder y = b. Ein anderer Vorteil, der fiir unsere Zwecke durch den Ubergang von 
(13a) zu (13a, .) erzielt worden ist, besteht darin, daB dadurch die Funktion Kg hinsichtlich ¢ 
zu einer geraden Funktion gemacht werden konnte, denn es ist bekanntlich [11]: 


Fire) Fyfe 2) 
H(r 2) H®(r’ 2) 
An der oberen Grenze ist hierin das Argument der rechtsstehenden BEssELschen Funktion 


gleich Null, an der unteren Grenze gleich (7 — 7) z. Mit den modifizierten BEssELschen Funk- 
tionen nimmt die Gl. (13a) die Form an: 


Dray We) 
THA? ALIS (G22) 
Die zweite Klammer des Integranden von (12) werde erweitert zu dem Lésungsansatz: 


(A < S 


(eis a a (cos #’) +P _ , (—cos ¥) 


2 


| 
| 


In(r’/r) 


J Gof ve)-JoevF+r2—2r7 -Colg)-de. (13a) 


te) 


In(r’/7) = ze _ 
= f Gof (r)-1) (zy? +72 arr’ .Cofy)-dy. (138) 


2 


| 
ge Pie)el seg (COS 0): +: Bt): Ps (—cos 0) | 
(cos #) | 


1 
t—+ t—+ 


|. P—*, (—cos o)«P 


(8 < 0’) 


= = K,(t; 0, 8's 0, B). (16) 
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Auch dieser Ansatz soll im ganzen Bereich a = <f gelten, so daB es mithin gerechtfertigt 
ist, die beiden fiir beliebiges ¢ voneinander unabhangigen Partikularlésungen ee 1 (+ cos 8) 


gleichzeitig in den Ansatz aufzunehmen. Da im Kegelmantel fir 0= «, = die Feld- 
komponente E, fiir alle in Frage kommenden Werte von 7 > a verschwinden muB, so wird 
gemaB Gl. (7a,) diesmal u(r, «) oder u(r, B) selbst gleich Null. Mithin ist jetzt zu setzen: 

B,(t)+ P_, (cos 0) =-" Ba) ages (a COS arma (— cos &) - ee (cos a) , (17a) 


t 


B,(t)- P_ 1 (cos) + B,t)- P,_ 1 (— cos peat (+ cos #’)-P__, (—cos#). (17b) 


Daraus folgt durch Auflésen dieses Gleichungssystems 


B,(i) =P, (—cos6) 


2 


P_1 (cos a) . P. (— cos a) EA Caso) eleven (— cos a) 


ae eee 
2 


x i. 2 , 2 : 2 2 F (174) 

| ere (cos #’) P_*, (—cos 0) P_ 1 (cos f) Fy palace B) 

pat Ee att (— cos f) 
2 

P_1 (cos) ae — cos B) Be aa (cosa) P_, (—cosa) 
x a ene ey : 17 B) 

zohan COCO Nine, 6088 | a(cosB) P_, (—cos 8) 

~ e ar | | 2) 


Die obige geschweifte Klammer (16) bekommt damit zunachst das folgende Aussehen: 


Lb, Uae ep) 


| Bae ate Bees nue x) pes x) Pan cece) 
2) = Pa cos) 1(—cosp) P,_i(cosB) — P,_ (—¢08) 
ir: een (— cos #) Ge cos?) P ‘3 (cos #) wie (— cos 0”) 
Py (05a) Pea (cosa) 
p en P eer @Zo<aw =f). (16a) 


Fira <& <¥ Sf lautet das Element a,, in der Zahlerdeterminante von Gl. (16a) 


+ ery (cos &) - Pre (— cos #) . 
2 


2 


Auch hier la8t sich wieder die dreireihige Determinante als das Produkt zweier zweireihiger 
Determinanten darstellen. In der Tat bekommt man im Falle 


aZo<w <p: 


Pi (cosa) P_ ,(—¢osa) 

= : | Bea(cos B) a= cos B) P ee 0) Fae ie cos #) 
Kp(t; 9, #5 a, B)= | ee (cos 8’) — Po, (— cos 8’) |" | Ps (00 ee Le cos a) | 

eae, 1(cosB)  P__4 (cos) 
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und im Falle 
Ca aes p 


Ie 1 (COS a) tee (— cos a) 
2 2 
[Pia (605.8) Py (—0088)| | pa, (cos) — P_, (0086) 
Bette 2) 10s) 2 2 Seta. — | 
P._ 1 (cos f) P_ 1 (—cos B) | P__ +{cos:a) P*1 (— cos a) 
2 es oN 
| P_ 1 (cos f) P_1 (— os B) | 


(165) 


In Gl. (16a,) verschwindet die Funktion Kp, wie es sein mu8, fiir & = « und in Gl. (16a,) 
fiir @ = B. Auch die Funktion K, ist wegen P_1=P_,_ 1 hinsichtlich ¢ eine int gerade 
Funktion. E eee 


1 


Mit den neuen, standigen Abkiirzungen Z, (¢;k; 7, 7; a, 6), N, (¢; a,b) und Z, (‘+ ; 
2 
OF +6, B), Np (+; On p) , deren Bedeutung durch unmittelbaren Vergleich der nach- 


stehenden Beziehung (18) mit den Gln. (13a,, 13a.) und Gln. (16a,, 16a,) unschwer hervorgeht, 
lautet nunmehr die formale Lésung der eingangs formulierten Aufgabe: 
x Ug: sin & 


u(r, 3) = + aa 
4 Yor’ 


cos mt 


[ i. Kats kit, 5a, b) Kplt; 8, 0’; 0, B) » dt 


© 


1 
atiaaet | ean ion ae) ai. (28) 
aa 7 ee ct Ge a we a ee foe Oe oe 
4 Vrlr Bl ) Np (> me 6) 


& 


Hierin ist © immer noch der im AnschluB an Gl. (12) beschriebene, nach rechts offene U- 
f6rmige Weg von Bild 1, der parallel zur reellen Achse dicht oberhalb und unterhalb von ihr 
verlauft und dessen linker zuriickfithrender Bogen bis an den Nullpunkt der ¢-Ebene heran- 
gefiihrt werden kann. Der Integrand des Integrals von Gl. (18) besitzt in seinem ersten Glied 
t/cos(a t) eine ungerade Funktion. Die anderen beiden Faktoren K, (¢) und K,(¢) sind hin- 
gegen, wie wir gesehen haben, gerade Funktionen von ¢. 

Alle Pole, die der Integrand in Gl. (18) zu besitzen scheint, sind jedenfalls einfache 
Pole. Sie kommen zustande durch die Nullstellen der drei Nenner des Integranden. Diese 
Nullstellen liegen, falls Im(k)) = 0 ist, bei . 

p= 1/2 +A mit A —.0,1, 2... im*Falle des Nenners cos(x ¢); 

2-t= +1, mit reellwertigem 1,, so da8 Np(+1,— 1/2; 0, 8) = oist ftir g = 1, 2,3..«; 

Zaebe= 57 Gy mit reellwertigem o,, so daB N;, (+12 Oy; Ro; a, 0) SONS tmit po oe oe 

3b. #= +e, mit den endlich vielen Werten Qj, Q2, 03, Q4 - - - Qm, SO daB Np (0p; Ro; 4, 6) 

= oist, falls die Wellenzahl k = ky und also reellwertig ist. 


Auch ein Teil der Polstellen der Polkette 3 scheint also wie die Pole von 1 und 2 in das 
von der Wegschleife © umrandete Gebiet zu fallen. Aber gerade der Umstand, da8 nur ein 
Teil der Nullstellen von 3 davon betroffen wiirde, macht diese Annahme von vornherein 
sehr verdachtig, und es fragt sich, auf welche Weise entschieden werden kann, ob tatsachlich 
im Hinblick auf die Anwendung des Residuensatzes, wo diese Frage von entscheidender 
Bedeutung ist, eine solche verschiedenartige Lage der Polketten io, und 9, zur Wegschleife © 
zu Recht angenommen werden muB. 

Gerade um diese Frage zu entscheiden, ist es zweckmaBig, neben dem besonderen Fall 
eines reellwertigen k) auch den allgemeinen Fall eines komplexwertigen k zu betrachten. In 
der Tat besteht zwischen den Polketten 1 und 2 einerseits und der unter 3 genannten Polkette 
der Unterschied, daB die Lage der Nullstellen 1/2 + 4 und 1, ganzlich unabhangig von der 


a 
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den EinfluB des Materials wiedergebenden Konstanten k ist, wahrend sie im dritten Fall 
durchaus davon abhangt. Sollte sich namlich zeigen, daB bei komplexwertigem k die Null- 
stellen der Funktion Nz (¢; k; a, b) in bezug auf t alle als auBerhalb der reellen Achse und damit 
auch als auBerhalb der Wegschleife liegend angenommen werden miissen, so werden diese 
Nullstellen wegen ihrer stetigen Abhangigkeit vom Argument der beteiligten Funktionen 
auch bei reellwertigem k = ky als auBerhalb des von der Wegschleife © umschlungenen Bereichs 
liegend anzusehen sein. Eine dementsprechende Wegfiihrung laBt sich nun aber auch bei 
reellen Werten von k jederzeit erreichen, wenn beide Halften des Schleifenweges die auf der 
reellen Achse liegenden Nullstellen 3b oberhalb oder unterhalb von ihnen umlaufen, wahrend 
die beiden Teile von © bei den unter 1 und 2 genannten Nullstellen bei jedem Umlauf um eine 
Polstelle den betreffenden Pol zwischen sich haben, wie es in Bild 3 angedeutet worden ist. 

Die eben beschriebenen Schwierigkeiten in der Wahl des richtigen Verlaufs des Inte- 
grationsweges treten nicht auf, wenn, wie in den einfacheren Fallen z. B. der Potentialtheorie, 
das die Anregung beschreibende Ausgangsintegral, in Gestalt eines gewdhnlichen Linien- 
integrals darstellbar ist. Wenn es mdglich ist, sollte man immer zunachst diese neben der 
Reihe einfachste Darstellungsform anstreben und mit ihr auch die Anpassung an die Rand- 
bedingungen vornehmen. Es ist aber hier nicht méglich, so zu verfahren; es sei denn, man 
umgeht diese Schwierigkeit durch eine nachtragliche Reihentransformation [10]. 

DaB die Nullstellen 3, wie bereits erwahnt wurde, bei komplexwertigem # nicht langer 
auf der reellen oder auf der rein imaginaéren Achse liegen, 14Bt sich unmittelbar aus der Glei- 
chung 


Ng (0p + 5; hy + 4A; a,b) Nolopis ha) +. (Np (63 Bo) |), + 4 A- (ON g(t; 2)/Ox)e—2, 


erkennen, die bei geniigend kleinen Werten von 6 und J wegen der stetigen Abhangigkeit der 
Nullstellen von Argument und Parameter ohne weiteres aus dem Taylorschen Lehrsatz folgt. 


Bild 3. Dieses Bild stellt den erforderlichen Verlauf des Integrationsweges © in Gl. 2 (18) dar, wenn die Wellenzahl k in 
der Funktion Ap reell, d.h. k = & ist. Er muB diese besondere Form haben, damit er nur die Nullstellen 7 umschlingt. 


Wegen weiterer Erlauterungen wird auf den Text im AnschluB an die Gl. 2 (18) verwiesen. Wird der Weg © an der imaginaren 
Achse gespiegelt, so 1aBt er sich zu einem Doppelweg oberhalb und unterhalb der reellen Achse von — coe«e -++ oo erweitern, 


Im iibrigen wird wegen der teils reellen, teils rein imagindren Nullstellen, die die Funktion 
Np (t, Xo; a, 6) besitzt, vor allem auf die Arbeit [10] yon G. N. Watson verwiesen. Wegen 
anderer hierzu interessierender Formeln vergleiche man die Arbeiten [11], [12] des Verfassers. 
Aus der schon dort angegebenen Integraldarstellung (13«) geht z. B. sofort hervor, daB fiir 
| (r’ —1r)- z| < Jo, Mit Jo(Jo,) = 0 bei reellem z und y noch keine der Nullstellen 0, besteht. 
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Aus den Gln. (18) und (16a,, a,) ist zu ersehen, daB die Polkette 1 in Wirklichkeit nur 
vorhanden ist, solange im Integranden fiir ¢ = 1/2 + A nicht auch die Funktion K, im Zahler 
verschwindet. Das ist aber gewiB nicht der Fall, wenn « = 0 und f = a ist, weil fiir cos 3 
> —1die Funktion P_, (cos 8) > oo strebt. Es sind dann alle Kegelflachen zwischen den 

2 


beiden Kugeln 7 = a und = 6b verschwunden, und die Funktion 


EGO Ona) P*, (cos w)- Sea (—cosd?) ("<< 9), (19) 


2 2 


reduziert sich wieder nach Gl. (16) auf den allein die Anregung beschreibenden, ganz links 
stehenden Term. 


Wenn auch in dieser Arbeit nicht die Absicht besteht, etwa in den Ergebnissen bis zu 
Zahlenwerten vorzudringen, so hat es aber doch im Hinblick auf die Beziehungen fiir die 
Resonanzfrequenzen des Schwingtopfes einiges Interesse, im Fall der allgemeinen Gl. (18) 
auch die aus der Integraldarstellung der Lésung entflieB8enden Reihenentwicklungen kennen- 
zulernen. Es empfiehlt sich dann aber, zunachst noch den Integrationsweg von Gl. (18) 
umzugestalten. Da der Integrand in Gl. (18) auf Grund der fritheren Angaben im ganzen ohne 
das Glied dé eine ungerade Funktion von ¢ ist, so fiihrt der Ersatz von ¢ durch —?’ im Inte- 
granden zu dem Integrationsweg —(’ in Bild 3. Auch diesen Weg hat man sich dann als 
einen Doppelweg zu denken, bei dem zwischen den einzelnen Polen Hin- und Riickweg zu- 
sammenfallen, der die auf der reellen Achse liegenden Pole = + T, mit dem Hinweg oberhalb, 
mit dem Riickweg unterhalb umfahrt, wahrend die beiden Hialften dieses Weges die Pole 
t = +0, entweder beide oberhalb oder beide unterhalb der Pole umlaufen. Die beiden zu- 
nachst noch getrennt liegenden Wege © und —€ diirfen dann sogar, wenn es zweckmabig 
erscheint, zu den Wegen ©’ + ©” zusammengefaBt werden, wie es Bild 4 erkennen 1aBt. 


t —Fbene 


Bild 4. Bei der Abstrahlung in den unbegrenzten Raum hinein reduziert sich die Nennerfunktion Kg auf die Funktion h() (ak) von Gl. (23a). 


Die Nullstellen liegen nunmehr wegen Gl. 2 (27b) im ersten Quadranten der t-Ebene, und der Weg © kann jetzt ersetzt werden durch den Weg ©’ 
oder ©” oder ©’ + ©’, wenn notwendig durch Aufnahme des Faktors 1/2. 


Der Wert des Integrals mu8 dann allerdings gehalftet werden. Die Verschiedenartigkeit 
der Umlaufe um die Polstellen +7, und +70, wurde auf diesem Bild absichtlich nicht 
hervorgehoben. 


Lést man nun in Gl. (18) etwa das Integral der zweiten Zeile zunachst nach den Null- 
stellen ¢ = +1, des Nenners Np auf, indem man den Doppelweg © auf diese zwischen ihm 
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liegenden Pole der ¢-Ebene zusammenzieht, so erhalt man ohne weiteres die folgende erste 
Reihendarstellung: 


1 
Zp\tz— — ; B, 5 4, 
wae s Tq Z pltqi kos ¥, 1’; 4, ©) >| PEE B ; 


u(r, ?) — —271- = Sill y. : -ko- ab A ae eee 
4 Vor" ean COS (7 Tq) Na(tq; Ro; a, ) (exp (x ia, p) /en), 
ed 


Zz 


(18’) 


Eine ganz anders geartete Entwicklung fiir u(r, #) entsteht dadurch, daB man sich in 
Bild 4 die obere und untere Halfte der beiden Wege © und —€ rechts und links von der 
imaginadren Achse durch Schwenken um die Winkel + z/2 an die imaginare Achse der t-Ebene 
herangebogen denkt und danach einen der beiden nunmehr senkrechten und mit entgegen- 
gesetzter Richtung durchlaufenen Wege iiber die imaginare Achse hinwegzieht und mit dem 
anderen vereinigt. Das erwahnte Herumschwenken um go° und das Verschieben der Vertikal- 
wege ist mit der Konvergenz der Integrale vereinbar und ohne Wertanderung méglich. Der 
Nachweis hierzu gelingt teilweise schon mittels der Formeln 3 (2a,b). Beim Hintiberziehen 
des linken oder rechten vertikalen Weges iiber die imaginaére Achse hinweg werden die dort 
liegenden einfachen Pole +7-o, iiberquert und beim Hochziehen der Weghalften von der 
reellen zur imaginaéren Achse rechts die Pole + @, und links die Pole —@,, und nach dem 
Residuensatz ergibt sich auf diese Weise als zweite Entwicklung der Ausdruck: 


co 


a U,-sind Op Z plop; Ros 7, 7"; a, 0) 
5 ? = Lf . Se ee . Sat pa eee Soe = 
u(r, P) +271 i Vrir os Cos(%q) (ON p(; ho; 4, b)/Ox)x = 0, 

1 

20 (oy %, 9: 6) oo ad Zetgnee 0,8) 

Sere = | net) OS Dees 2A 

1 Coj ao 1 

Np(er— a 8) p=o0 b Np(iap— i 8" x, B) 

SS ENN) (18) 


(ON p(x; hy; 4, b)/Ox)x=+ icy 


denn es sind natiirlich die Residuen in den Polen x = +74, fiir gleiche Werte von # entgegen- 
gesetzt gleich. Wegen der Eigenschaft, gerade zu sein, sind die Funktionen Z, und N, fiir 


t =ic,—— und fiir Z, fiir ¢ = ia, reellwertig. Dagegen ist die Funktion 0N (Ox [Ur == ae, 
p F B p § §es B p 


rein imaginar. 


2.3. Das Strahlungsfeld eines magnetischen Stromringes in den Kugelflachen a =r 
oder b oder in den Kegelflachen # — « oder £ 


Fiir das Hineinriicken in die Kugeloberflache 7’ = a fiir ein # aus dem Bereicha < 0 < B 
ist die Gl. (13a,) maBgebend. Im Grenzfall ist in dieser Gleichung 7’ = a zu setzen. Es be- 
kommen dann alle Zylinderfunktionen der maBgebenden Determinante dasselbe Argument. 
Macht man daher von den Definitionsgleichungen (15a, b) Gebrauch und beachtet eine be- 
kannte Beziehung aus der Theorie der Zylinderfunktionen, so erhalt man schlieBlich den 
Ausdruck: 


ERCP R700, 0) 


ni [Tr k) HM i : (130) 


aYak |j(bk)  W(b R) 


jak) h(a k) 
H(b&)  hy(b R) 
Erfolgte der Grenziibergang von’ gegen die 4uBere Kugel mit dem Radius b, so ware in Gl. (13.4) 


auf der rechten Seite (a k)1/? durch (b k)!/2 zu ersetzen und in der zweiten Zeile der Zahler- 
determinante im Argument der Funktionen j, und h der Buchstabe 6 durch a. 


XLV. Band 


Hietta=— 1960 H. Bucunoiz: Elektromagnetische Felder mit axialer Symmetrie 39 


LaBt man andererseits den Stromring # = & in die Kegeloberflache 9 = a iibergehen, 
so hat man in Gl. (16a,) in der letzten Zahlerdeterminante # = a zu setzen. Es empfiehlt 
sich dann, zundchst mittels der ersten der beiden folgenden Beziehungen 


Pe (— x) = cosa (s + pw) - P#(x) —2/n-sina (s+ m)-Qe(x) (0o<%< 1), (20a) 


P(x) 0) __ 0) PO aja) (-1<xe+1), (20d) 


innerhalb der in Rede stehenden Determinante durchweg zu dem Argument + cos a iiber- 
zugehen und gleichzeitig den oberen Parameter —- 1 vermége Gl. (11b) in den Parameter + 1 
zu verwandeln. Es kann dann auch die Ableitung der Funktion P_1(cos «) selbst ein- 


gefiihrt werden. Im Hinblick auf die Gl. (20b) ergibt sich so die Relation: 
| Peet (cos a) P 


1 1 
t 3 ra 


os an 2) | cos x t 2/% (21) 


= + 


P~*, (cos a) = P54 (—€0s 9) #—i/4 sina 
a 


eek 
2 


Im Falle des in der Kegeloberflache # = « liegenden Stromringes ist mithin: 


2/m% + cos 2 t 


#— 1/4)-sina 


PoE PO Op) = | 


P,_1(cos 0) eee cos #) P_ 1 (cos x) een ae COs a) 
x Bed (cos B) 1 ee cos f) P 1 (cos B) er (— cos B) 
(a<9< 68). (16.4) 


Liegt der Stromring in der Kegelflache @ = B, so muB in der letzten Gleichung statt sin a im 
Nenner des Bruches vor der Determinante sin B treten und unter Umkehr des Vorzeichens in 
der ersten Determinante rechts der Winkel § durch a ersetzt werden. Fiir # = f verschwindet 
Kp nach Gl. (16a) identisch. Fir #—=a tritt, wenn die Funktion Kp, von Gl. (16a) 
unter dem Integralzeichen steht, gleichfalls Verschwinden ein, da dann auf Grund der friheren 
Ausfiihrungen innerhalb der Wegschleife © iiberhaupt keine Singularitaten liegen. 

Ist die kleinere Kugel y = a mit der gréBeren Kugel y = b nur durch den Kegelmantel 
@ = a verbunden, so vereinfacht sich die Gl. (16a) wegen P;_,/,(x) > oo fiir x > — 1 weiter 


zu der Formel: 
1 (-- cos ?) 


Be 
2/m+-cosmt Se 
EGR Ose nase) =: ; 

(t#2— 1/4) -sina& pip den Ol COS &) 
2 


(16.4) 


Bei dem besonderen Aufbau von Ky; in den Gln. (16a) und (164) fallt in diesen beiden Fallen 
die polerzeugende Funktion cos 7 ¢ im Integranden von Gl. (18) heraus. 


2.4. Das Strahlungsfeld des auf Kugel oder Kegel aufsitzenden magnetischen Stromrings 
in den freien Raum hinaus 


Wir wenden uns nun dem uns hauptsachlich interessierenden Fall zu, in dem die Aus- 
strahlung von dem magnetischen Stromring # = « aus in den freien Raum hinaus erfolgt. 
Die diesem Fall entsprechende allgemeine Lésung (18) setzt sich, was den Faktor kK, anbetrifft, 
aus Gl. (16a) und betreffs der Funktion K, aus Gl. (13a,) zusammen. Es mu nur darin 
noch der Grenziibergang 6—> oo vollzogen werden. Unter der voriibergehend gemachten 
Voraussetzung Im (k) > 0 strebt aber fiir b> oo die Funktion h® (bk) +0, dagegen die 
Funktion 7, (6 k) > oo, und daher ist im Grenzfall: 

Hk) \J,¢' k)- BY (r' k)| 


eee ee (22) 


Teh Tt aa, =: F, 
=| Ls Bw oo) h(a k) BAC k) AY (a k) | 
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Es versteht sich, daB in Gl. (22) hinterher wieder & rein reell sein darf. Aber auch bei reellem 


k = ky strebt K, in Gl. (22) fiir 7 > oo gegen Null. 
Ist vorliufig noch B £ « und B £2, so ergibt sich damit auf Grund von Gl. (18) die folgende 


SchluBgleichung: 


Fx) Fase) 
ee 2 
u(r.) = "ee [| HP | es) ee 
2yrlr' J} |4,(ak) AY (ak) | |P,_a (cosa) P,_ 1 (— cosa) | 
¢ ; 
P,_1(cosB)  P,_1 (— cos §)| 
2 2 : 
BEND ES Gabe ph ae eee (23) 


a (ak) #®#—1/4 
Fiir 0 <7 <7 < ow hat man nur nétig, in Gl. (23) unter dem Integral, aber nur dort, 7 mit 7’ 
zu vertauschen. 

Denkt man sich in Gl. (23) unter dem Integralzeichen den Quotienten aus den Funktionen 
H® und h') in Zahler und Nenner mit dem Faktor exp (¢- 7/2) erweitert, so wird deutlich, 
daB der Integrand eine in ¢ ungerade Funktion ist. Ersetzt man also im Integranden ¢ durch 
—1?’, so lautet der Integrand in ¢# mit Einschlu8 des Vorzeichens genau ebenso wie in Gl. (23) 
der in ¢, nur verlauft jetzt der Weg © wie der in Bild 4 dargestellte Weg —. Die beiden 
Wege ©, —€ lassen sich dann wieder wie friiher zusammenfassen. Die teils reellen, teils 
rein imaginaren Polstellen ¢ = 9, und ¢ = +70, der Polkette 3 in dem Text zu Gl. (18) sind 
jetzt verschwunden. Ihnen entsprechen neuerdings die unendlich vielen komplexen Null- 
stellen der Funktion h{! (a k) in bezug auf t. Sie haben etwa die in Bild 4 angedeutete Lage 
in der ¢-Ebene. 

Der Vollstandigkeit halber und in Riicksicht auf die spatere Verwendung werde die Dar- 
stellung (23) noch erganzt durch die Gleichung 


tee (— cos #) HS (v k) fh 


oe fel (7) an ie Ae 
P,_1 (— cosa) WO (ak) 2—a1/4 


= | 23a 
2} 5, (ak) WY (ab) ae 
& 


u(r, 2) 


(6 OS CO 0 ot) 


fiir den besonderen Fall 8 = a und durch die Gleichung 


web) ee of EO aie rags (23b) 
, ark bPiak) Po we cosa P77, 3 


2 


C 
(Gis OS) a OOO i eee ere) 


fiir die weitere Spezialisierung von 7’ in (23a) auf den Wert a. Im Falle der Gl. (23b) liegt 
also der magnetische Stromring in der Nahtlinie zwischen dem Kegel & = a und der Kugel 
ry =a. Auf diesen Sonderfall weist auch SCHELKUNOFF [13] hin. Auch in (23a, b) kann € 
durch den Weg ©” in Bild 3 ersetzt werden. 
Riickt andererseits der Stromring in die Kugeloberflache, so bleibt #” 4a, wahrend 
r’ = awird. Mithin ist dann 
21 


Kp =". HO(r Ryn (ah) 


mak 


Wird noch auBerdem 6 = z gesetzt, so ist jetzt nach Gl. (18): 


—] , : 
Oe: Uy-sin 8’ [eee Prat 08 8) P,_1 (cos x) fae Cos a) frat 
a a) : © 
2k. a ¥ ¢ h; (a4 hk) Prt (— cos a) Py (cos D) Pie (— cos 0) cos (xt) 
Ss 2 
(234) 


(a<0<W <n; rf =a; pHa). 
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Fir ein a = # < 9<a steht unter dem Integralzeichen vor der Determinante im Zahler 
des Bruches P,_1 (—cos #) und in der zweiten Zeile der Determinante in derselben Reihen- 
folge ; 

8 Be (coe a eR sts (= cos) . 

2 2 

Macht man hierin wieder # = a, so erhalt man wie vordem die Gl. (23b). 

Schrumpft die Kugel mit dem Radius a auf Null zusammen, so gilt nach Gl. (23a) fiir das 
fiir den magnetischen Stromring in der Kegeloberflache 9 = a erzeugte Strahlungsfeld die 
Darstellung: 


P,__1 (—cos #) 
U t t+ dt 
) ee 4 (1) ») . u 25 fn ee = 
u(r, 9 Bre i Te kh) HOP B) Pi1(—008a) P= /4 (238) 
& 

(Ons Fee Faso Oa aC yp LO, oe OO ah) : 

Wir gehen in (23a) vom Integral zur unendlichen Reihe iiber. Sind dann ¢ = 7, wie friiher 
mit g = 1, 2,3... die unendlich vielen reellen Nullstellen der Gleichung 

P;,(— cosa) = 0, (24a) 


in bezug auf t, so ergibt sich durch Zusammenziehen des Weges auf die positiv relle Achse die 
Reihenentwicklung : 


ee) Srp atte ee Bee eles 
Vrir ay (t+ 1) Weg + + (a k) (@P,(— cos a)/0%)x = 19 
ENA (r’ k) test (r’ k) (a <0 <z) tay 
jg (2B) W@W] (@< 7 <1 <o0) 


Fiir ein 7 <<’ sind in der letzten Gleichung tiberall y und 7’ miteinander zu vertauschen. 
Bei den Winkelbereichen fiir a, wie sie im Falle des Kegels als Antennenleiter in Frage kommen, 
kann in Gl. (24a) in erster Naherung die Berech- 
nung der Werte t,nach der einfachen Formel 


1 4 
Tq = q -+- a . arctg (n=) 


Spree ae 
een cg are te( canis) 


mwq+f(2-In(2/a)) (¢=1,2,3...) (24a) 


erfolgen. Eine zahlenmaBige Ubersicht iiber die 
Werte von 7, fiir g=1...13 in Abhangigkeit 
von a im Intervall 0° = « < 180° gibt in anderer 
Bezeichnung die Abb. 5. 

Statt des Weges © kann in Gl. (23a) aber auch 
der Weg ©’ in Bild 4 gewahlt werden. Es ist dann 
ohne Wertiinderung moglich, diesen zur reellen 
Achse parallelen Weg durch einen sehr groBen 
Halbkreis in der oberen ¢-Halbebene zu einem 
geschlossenen Weg zu erganzen. Dieser Weg um- 
schlieBt jetzt nicht mehr die Polkette 7, von 
Gl. (24a), sondern die Polkette v 


die durch die 0 


n? 


o° 20° 40° ~=—60° go? 100° 120° 140° 160° 180° 


komplexenNullstellen’y, mit w= 1,2, 3°. .< der f pea 
Funktion Bild 5. Dieses Bild veranschaulicht die Verteilung der Null- 
d stellen s=z, der Funktion Ps (cos 6) =o in bezug auf s 
(1) eS 1/2 1) a bei konstanten Werten von coso innerhalb des Winkel- 
h;, (a Rp) — (- (x 2. HN (x)) ews = @ (25a) bereiches 0° S 6 S 180°. 
: an 0 


in bezug auf den Parameter ¢ in Gl. (23b) verursacht wird. Auf die Berechnung dieser Null- 
stellen v,, die bereits in einer neueren Arbeit von W. FRANZ und R. GALLE [14] angenahert 
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durchgefiihrt worden ist, wird weiter unten im Zusammenhang mit der Entwicklung anderer 
hierher gehériger Formeln eingegangen. Werden die Werte », vorlaufig als bekannt angesehen, 
so berechnet sich fiir die Funktion u(r, 9) noch ohne jede Vernachlassigung die folgende 
unendliche Reihe: 


2504) Uaieg Nese il Were ty (0) Re hy) HY (=. hy} 
u(r, 3) = area Dear? Hp (=. a tm (2 @ Ro 


<r Pup toen a 
(dh? (a ko)/@)1—v,, Va es (— cos a) 


denn wegen 2+ 47,(a hy) = h(a ky) + h(a kj) reduziert sich in Riicksicht auf Gl. (25a) 
jy, (@ ko) auf die Funktion hi (a ho). 

Soll die Aufgabe auch numerisch weiterverfolgt werden kénnen, so braucht man nach 
Gl. (25) Formeln sowohl fiir die Nullstellen selbst als auch fiir die Funktionen h!) und H!, 
wenn der Parameter ¢ = »; ist. Mathematisch ist dies wesentlich die gleiche Aufgabe wie die, 
der sich G. N. WATSON in seiner unten erwahnten Arbeit gegentibersah. Die hier vorliegende ' 
Aufgabe ist nur leichter insofern, als das Argument ak, im jetzigen Falle rein reell ist. Sie ist 
aber insofern wieder schwieriger, als der Betrag von ak hier in dem wesentlich niedrigeren 
Zahlenbereich von etwa 50---500 liegt. Die Berechnung der Nullstellen ¢ = »,, von Gl. (25a) 
muB8 daher jetzt mit gréBerer Annaherung erfolgen, als dies in der Arbeit [10] notig war. Der 
hier zu stellenden Forderung nach gréBerer Genauigkeit bei der Berechnung der Nullstellen 
von h\) (a ky) entsprechen am ehesten die Rechnungen in der schon oben unter [14] genannten 
Arbeit. GemaB den Gl. (ga) und (gb) dieser Arbeit auf Seite 375/376, die sich auf die Funk- 
tionen H'(x) beziehen, berechnet sich fiir h')(x) die Naherungsdarstellung: 


WN og) A er eae 


pe. * ag (28 + 2B) 4] 


12 27 +280 3780 243 - 800 
ss 38-9? cy RR g ; 
a 12 te 81-700 a 2835 21600 on 243 - 108 =] ay 
1 148 - 61-¢ op site 
=p (= - 700 me 27:27 +700 = Dif oo ecos) A (2) + | ‘ Gee 
In Gl. (26a) stehen die GréBen », x, y, g in dem Zusammenhang: 
Ga amet) (26a) 
y = els. (%/6)— 8 ; (268) 
und es ist A(g) die Arrysche Funktion, die entweder durch das Atrysche Integral 
+ co 2. 
ee a cee —g—3i a Baa 
A(q) : { exp [—7 (*8— q-1)]- drt = 3-3. [ es (+ t ae 1). di (26A) 


oder durch die Differentialgleichung 


d?A (q) q abe 2% 2 y a 1 

ae ty 4@=0, 40)=2/[r(Z), a =Z/r() 268) 
definiert werden kann. Sie ist mitsamt ihrer ersten Ableitung nach g gema8 den beiden 
Gleichungen 


Aq) = (a13)"* (Fars (2+ (913)°*) + Jus (2- (13M = 4, - 4i(—F), (e6c,) 
Ag) = — 5-9 Van: (a13)"") — Jal? (/3)*")] = — 5 av(—5) (ese, 


durch die BEssELschen Funktionen J £13 und J, ,;, darstellbar. 


Res H. Bucuuorz: Elektromagnetische Felder mit axialer Symmetrie 43 


In der der Arbeit [14] entstammenden Gl. (26a) ist A(x) in eine nach Potenzen von y? 
oder x—*/3 fortschreitende Reihe entwickelt, deren Glieder umso kleiner ausfallen, je gréBer x? 
ist. Die entsprechenden Reihen fir H®(x) und H®’(x), auf die oben Bezug genommen 
wurde und aus denen sie sofort hervorgeht, werden in der oben erwadhnten Arbeit aus den 
SOMMERFELDschen Integraldarstellungen mittels der Sattelpunktsmethode abgeleitet. Nach 
der Auflésung des Beifaktors im Sattelpunktintegral nach Potenzen der Integrationsvariablen 
lassen sich die einzelnen, neu entstehenden Integrale alle durch das Airysche Integral (26 A) 
fiir die Funktion A(q) selbst oder deren Ableitungen ausdriicken. Wegen der Gl. (26B) kénnen 
aber sdmtliche héheren Ableitungen von A(g) von der zweiten an durch die nullte und erste 
Ableitung dargestellt werden. Nach diesen Angaben diirfte der Aufbau des Ausdruckes (26a) 
im wesentlichen verstandlich sein. 

Aus der Definitionsgleichung (26C,) fiir A(g) ist nun aber sofort zu ersehen, daB die Funk- 
tionen A(g) und A’(q) hinsichtlich g unendlich viele reelle Nullstellen haben. Sie mégen 
mit &, bzw. &, bezeichnet werden, so daB z. B. die Relation gilt: 

Ale Nes Cages Dap ete ele (27a) 
Diese Nullstellen sind, was die €, anbetrifft, in dem bekannten Buch [15] von G. N. WATSON 
auf S. 751 zu finden. Von den Nullstellen der Funktion A’(q), auf die es hier ankommt, sind 
die ersten zehn, da sie in der Literatur gar nicht so leicht zu finden sind, in Zahlentafel 1 
zugleich mit den Werten der Funktion A(é,) auf 6 Dezimalstellen genau angefiihrt. Sie 
wurden nach den Angaben im Tafelwerk [17] berechnet, das u. a. die auf 8 Dezimalstellen 
genau berechneten Werte der ersten 50 Nullstellen der beiden AtRyschen Funktionen A? 
und Az der Gl. (26C,, C,) enthalt. 


Zahlentafel 1. Die ersten 10 Nullstellen g = &, 
der Funktion A’(q), Gl. (26A). 


q e A(é;) 


1 + 1,469 353.7 +1,166 798.8 
2 + 4,684 711.6 —0,912 724.1 
3 + 6,951 786.0 | +0,828 623.7 
4 + 8,889 027.4 | —0,779 616.1 
5 10,632 519.5 | +0,745 620.7 
6 12,242 516.3 —o,719 862.7 
Ff 13,752 497-3 +0,699 269.6 
8 15,183 513.7 —o,682 198.8 
9 16,549 895.5. | +0,667 674.0 
10 17,861 955.7 | 0,655 069.5 


Ist x sehr groB, so liegt nach Gl. (26a, f) bereits in g = &, selbst, und somit ausgedriickt 
in vy, in dem Wert 
y, =e + Ee et #3. (x/6)13 
eine erste Naherung fiir die Nullstellen »y = v, der Funktion h')(x) in bezug auf ¢ bei festen 
Werten von x vor. Man kann diese Formel fiir die Nullstellen von A‘"(x) aber offenbar 
noch verbessern, indem man fiir y eine Potenzreihe in y? oder x—3/? mit unbekannten Koeffi- 
zienten ansetzt. Diese zundchst noch unbekannten Koeffizienten treten dann auch in den 
einzelnen Potenzen von g innerhalb der Gl. (26a) sowie in den Funktionen A(q) und A’(q) 
auf. Geht man mit den daraus herstellbaren TAYLoRschen Entwicklungen um die Stellen é, 
herum in Gl. (26a) ein, so lassen sich dann die Koeffizienten in der angesetzten Reihenent- 
wicklung bestimmen, und man findet schlieBlich die Formel: 
3 Set Sot ene 
800 + (x/6)3/3 i 9 ai ae (x/6)5/3. 16000 


exp (+ 77/3) 1 ie ian, oO te TRE gt er 5O2 a 
E15 £2 567 n F n ]° 
b 


= 
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Auch sie ist bereits in der oben erwahnten Arbeit [14] als Gl. (14) aufgefiihrt. Da % reell ist 
so liegen mithin die Nullstellen y, in der oberen ¢-Halbebene. 

Auf welche Feldkomponenten es nun auch immer spaterhin ankommen mége, nach Gl. (25) 
sind auBer den Werten von y, auch noch Formeln fiir die Funktionen h (ako) und éh\” (a kp) /6t 
an den Stellen ¢ = y, erforderlich. 

Die Berechnung der zuletzt erwahnten GréBe erfordert zunachst eine Differentiation des 
Ausdrucks (26a) fiir h(x) nach » bei festem x und damit auch bei festem y. Nach Gl. (26a) 
ist aber offenbar 0/év = y+ é/ég, so daB es mithin nicht schwierig ist, die Differentiation aus- 
zufiihren. Da hinterher der Wert der Ableitung nur an der Stelle y = »,, d.h. fiir ein g = &, 
benotigt wird, verschwinden wegen (27a) in dem differenzierten Ausdruck alle Glieder mit 
A’(é) als Faktor, und man bekommt zwischendurch zunachst die Formel: 


r y ‘3 
(1) Borin Way goal bor eae) Pe A Pe Ay) 
(oh, (*)/0”), <1, aa 1 61/2 fn i 180 y 4 e y 11340 48 600 


Bey a ee a Oe ee LS BO. se 8 eral et (eaves 
4 y ( Sips 168 - 8100 Sn Cp ECaoon é, | (S,) - (28a) 
Beachtet man hierin noch den Zusammenhang zwischen y und x auf Grund der Gl. (268), so 
ergibt sich endgiiltig fiir die GréBe dh /ét im Nenner von Gl. (25) die Naherungsdarstellung: 


(ah (x) /6r), —, = a - A(é,)- [ —€é. 


g'2 _ (6/#) 413 e738 ‘ (s =) (6/x) 9/3 (— 1d ay, ae ae 8954 ES) ae | ee (28) 


3240 i Go 3600 Men miee 102060 =” 


(6/2)713 | 
oe 


e2% 1/3 


Die andere, noch fehlende Formel fiir h(a ko) 1aBt sich etwa auf Grund der Identitat 
WE (aig) = Ii. Cake) Im(k)) = 0 (29 a) 


berechnen. Darin bezeichnet der iiber y, und der iiber dem ganzen rechtsstehenden Ausdruck 
gezogene horizontale Strich, daB einmal », durch v, zu ersetzen und danach fiir die ganze, 
darunterstehende GroBe der konjugiert komplexe Wert zu nehmen ist. 


Um hierin 45) (a ko) zu berechnen, kann die Gl. (26a) herangezogen werden, in der dann 


CE BN Pare (30a) 
We ETS Hao Ya 2, (30b) 
X= ake (30c) 


zu setzen ist. Entwickelt man daraufhin mit der angegebenen neuen Bedeutung von q in 
Gl. (26a) die Funktionen A(g) und A’(g) nach Potenzen von y und faBt alle Glieder mit gleich 
hohen Potenzen von y? zusammen, so ergibt sich nach langerer und ziemlich miihseliger 
Rechnung, die aber wegen ihrer zweifelfreien Durchfiihrbarkeit hier nicht wiedergegeben zu 
werden braucht, unter Achtgabe auf Gl. (29a) die folgende SchluBformel!: 


h® (a by) | (2 esi) = AE. et2nils) ya. Sn ptanils , A'(é, et 2xil3) 


a 
ia ag 61/2 45 


2400 Sh 63 
pea et 6xi/3 3} ! Gr 122 Comme 89 6 
12.000 2 Ei 2Gl 27) 63g oie eee 
rae HG, eee) (30) 


Hierin hat y wiederum die Bedeutung von Gl. (30b). > 


1 Fiir die Vornahme der sehr zeitraubenden und miihseligen Kontrollrechnungen an den Gln. (28) 
und (30) méchte ich gleich an dieser Stelle Herrn Dipl.-Ing. ConsTANTIN meinen besten Dank sagen. 
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3.1. Die Beziehungen fiir die Feldkomponenten E,, E,, und H,, bei Einspeisung durch 
einen zonalen Spalt im Kegelmantel. Das Integraltheorem von KonrorowicH-LEBEDEW 


Erfolgt die Erregung durch einen zonalen Stromring in der Oberflache des Kegels wie in 
Bild 6, so hat als maBgebend die Gl. 1 (23a) zu gelten. Nach den Gl. (2) ist dann im Hin- 
blick auf Gl. (7a,) und Gl. (21) 


P,__1 (—cos #) () 
PN (720) Oe (A eg hye es ee ea (1a) 
24-Vrlr' Peas (— cos a) hy? (a k) 
& 
ie (—cos 8) pcp (r k) Pedi 

ie 2 ile ; 

alt, 0 tk ke [a rk, ak) Poin Pe nee El, (1b) 

2 


P}_ 1 (— cos 8) 7p (v R) t- dt 
, 2 t ek , 
: [4 a) [= (—cosa) hiP(ak) #—1/4 (29) 


@ 


P(r, 3) 


Yu 


Bild 6. Dieses Bild gibt die Anordnung von Kegel und Kugel an, wie sie den Rechnungen im Abschnitt 3 
entspricht. Es ist zu beachten, daB in Bild5 s = 180°—a ist mit der Bedeutung von « wie im Text. 


und hierin ist weiter 
J, (7 R) H\(r' k) 


Airy ky a k) = 
AS CaP OO 


(eee St a CY, OO 1) = aa dE 


Fiir 7 <7’ sind unter dem Integralzeichen die Veranderlichen 7 und 7’ zu vertauschen. Der 
Weg © umschlingt gemaB Bild 2 die reelle Achse der /-Ebene mit den Punkten + 1/2, + 3/2 
AES el 

Fir @ =a und a—o reduziert sich der Ausdruck fiir E,(r, «) zufolge Gl. (1a) auf die 
Form: 


{E,(r, 6) }a=0 ae ; | t+ J, (rk) » H® (rk) «dd (<n). (2) 


€ 


Nun gelten bekanntlich, falls bei endlich bleibendem z = |z|-¢'* in ¢ = |t|- e'* die GréBe |t| +00 
strebt, fiir die Zylinderfunktionen J,(z) und H(z) die beiden folgenden asymptotischen 
Darstellungen: 

It| 


[J.(e)| ~ (2: |e) - exp 


—|¢|- cos T- In (7 ao ek sin T+ «—o} 
(—a<1,0< +2) (2a) 
[2209() | ~ (2la| ey"? - exp {+ |#]-c08 © In (21) [2] sine (@—9) 


(--2< ae +2). (2b) 
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Die erste folgt unmittelbar aus der fiir ¢-> oo geltenden Naherungsgleichung 
T(z) ~ (2/2)'/D (t + 3) (jarc(t)| <2), (2a) 
die nichts anderes als das erste Glied der Reihenentwicklung fiir J,(z) darstellt, die zweite 
aus der bekannten SOMMERFELDschen Integraldarstellung [18] 
y=t-Cofa 


ice me aes y.(a—Tqa) + 3 rt Re See 
(2) (== Sin ) ‘ olny (2 ae Wey ) a 


Zz 


die hier nur fiir den einfachen Fall |y|S>j|z| in Anspruch genommen zu werden braucht. 


Geht man mit den Gln. (2a, b) in Gl. (2) ein, so zeigt sich, daB das Integral iiber den Weg © 
fiir alle y 47 auch jetzt noch absolut konvergiert, und da in (2) der Integrationsweg, ohne 
singuldre Stellen zu iiberqueren, auf einen beliebig kleinen Kreis zusammengezogen werden 
kann, so ist fiir a = o die Vektorkomponente E,(r, «), solange r £7’ ist, identisch Null. 


Um zu dem Wert des Integrals in Gl. (2) fiir y = 7’ zu kommen, l4Bt sich das Integral- 
theorem von LEBEDEW/KONTOROWICH [19] heranziehen. Hiernach ist 


(kn) = +— festen([ 10 HP) (kr) 802 (3) 


—140 0 


(k= 1h | eh 25" Oar, CIC) = x0 Oe 


Ist nun fiir einen sehr kleinen, positiven Wert ¢<7’ 


ene (OC. Set oP OT a) 
C (’ —eZ2rsrte) 
so gilt unter Berticksichtigung von Gl. (3) der folgende Grenzwert: 
+100 r' +e 
lim | Baul | H\) ea ea = 4, (3a) 


Fir y = 7’ kann aber in der Tat der Weg © im Integral (2) zu dem Weg —ico +--+ +100 
von Gl. (3a) aufgebogen werden. Das Verschwinden des Integrals von Gl. (3a) fir r 47’ 
erfolgt jetzt durch standigen Vorzeichenwechsel. Fiir ¢— 0 ergibt sich an der Stelle 7’ der 
oben angegebene Wert. 


Es wird dann damit im Hinblick auf Gl. (2): 


5 / 
{E,(r, Cs ao ae He (7 =f, B= a) ; (4) 


Drape ks 2€ 


Dieser Wert entspricht den Erwartungen, denn hiernach erscheint die Feldstarke E, am Spalt 
des Kegelmantels als die durch die Spaltbreite 2¢ dividierte Spannung U, zwischen den 
kugelférmigen Begrenzungsflachen 7’ + ¢ des Spaltes. 


3.2. Die durch den Spalt abgestrahlte Energie 


Da die gesamte abgestrahlte Energie aus dem Inneren des Spaltes heraustritt, so ist es fiir 
die Berechnung dieser Energie das Einfachste, den EnergiefluB am Ort des Spaltes selbst zu 
berechnen. Offenbar ist dann: 


Sp =a27r-U,- H,(r’, a) -sing. (5) 
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Unter der ausdriicklichen Voraussetzung eines nunmehr reellen Wertes von , ergibt sich 
dann einmal fiir S, die Reihenentwicklung: 
Ss = 7? w &-sina+ U2-/7'3/a . » 


q=1 


Taq(% Ro) HY (r’ Ro) 
Iq (a Ro) hY(a Ro) 


t+ 1/2 Pl,(—cosa) — Hag + (0 he) 
T° (t+ 1) (AP, (—cosa)/dx)z = cq hey (aho) ’ 


i 


x 


denn fiir ein solches h, ist ic ioeer (7 hey\ te Ay) 4 Role. 


Ersetzt man in Gl. (1c) den Weg © durch den oberhalb der reellen Achse verlaufenden 
Weg ©’ in Bild 4, so kann man jetzt das Integral in Gl. (1c) dadurch auflésen, daB man den 
Integrationsweg in der oberen ¢-Halbebene durch einen Halbkreis zu einem geschlossenen 
Weg erganzt und dann auf die dort liegenden Pole von 1/h) (a ko) zusammenzieht. Man 
erhalt dann 


Cc 


=* ae Det t (L) (4/ 
Sp = — 7? w &-sinasU2-)r'3/a . oe Fig 0 Bo) Hog (7° Fe) 


se a ky) hy. (a Ry) 


n=1 
x Ho, (a Ro) . Vy : ye — 2) (6b) 
(ar3>, (a ho) /O%)x =»,  v2—1/4 1g eres (= os a) 


In der angegebenen Schreibweise gibt allerdings die linke Seite nicht den komplexen Wert 
von S, an wie in Gl. (6a), sondern den konjugierten komplexen Wert von S,. In der 
ersten Spalte der Determinante dieser Gleichung diirfte es noch in der Regel eine Verein- 
fachung bedeuten, wenn die Besselsche Funktion durch die Hankelschen Funktionen zweiter 
Art ersetzt wird. 

Die in Gl. (6b) auftretenden Zylinderfunktionen sind im Falle einer numerischen Be- 
rechnung durch die schon friiher hergeleiteten Naherungswerte zu ersetzen. Ebenso kénnen 
in (6b) die jetzt komplexwertigen Kugelfunktionen naherungsweise aus der Gleichung 


Py—t (cos y) =~ fener a4) 2 er HA) , (2 my-sin mae , (7) 


berechnet werden. Die rechte Seite stellt hierin das erste Glied einer asymptotischen Ent- 
wicklung fiir die Funktion P,(cos y) bei groBen Werten des Index dar. Diese Voraussetzung 
kann fiir die Werte », in den Gln. (6a, b) als erfiillt angesehen werden. 


SchluBbemerkungen 


In der vorliegenden Arbeit wurde in der Regel von den beiden Winkeln a und B angenom- 
men, daB « ein spitzer und f ein stumpfer Winkel ist. An der 4uBeren Form der Lésungs- 
gleichungen Andert sich aber nichts, wenn a und f beide in dem Winkelbereich o und z/2 
liegen. Sind dann keine Kugelflachen vorhanden und ist tiberdies a < & < B, so liegt der 
Fall des kegelf6rmigen Doppelhorns vor. Dieser Fall wurde bereits in den dlteren Arbeiten 
des Verfassers behandelt. Die Entwicklungen der Lésungen nach den Eigenwerten 1, der 
Kugelfunktionen sind im Fall des Horns viel besser brauchbar als im vorliegenden Fall der 
Antenne, weil die Nullstellen t, nach Bild 5 dann in viel groeren Spriingen gegen oo streben 
und dadurch die Konvergenz der entsprechenden Reihen wesentlich besser ist. 

In einer weiteren Arbeit werden die hier aufgestellten Gleichungen dazu Verwendung 
finden, um den EinfluB der dem Kegel aufgesetzten Kugel auf die Intensitat der ausgesandten 
Strahlung unter verschiedenen Bedingungen zahlenmaBig zu bestimmen. 
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Der Warmewiderstand zwischen Zahn und Joch 


Von 
J. Hak, Briinn 


Mit 8 Textabbildungen 


(Eingegangen am 20. Mai 1959) 


Ubersicht. Es werden Unterlagen zur Bestimmung des Warmewiderstandes zwischen dem Zahn und dem 
Joch einer elektrischen Maschine gegeben und gezeigt, welchen EinfluB die Zerlegung des Eisens in zwei Warme- 
quellen auf die Genauigkeit der Berechnung von Warmequellen-Netzen ausiiben kann. 


1. Einleitung 


Wenn man im Warmequellen-Netz einer elektrischen Maschine das Stator- oder Rotor- 
Eisen durch eine einzige Warmequelle abbildet, so ist dies eine Naherung, bei der man die ver- 
schiedenen Eisenverluste sowie die verschiedenen Querschnitte im Zahn und Joch vernach- 
lassigt. Die zur gréBeren Genauigkeit fiihrende Zerlegung des Eisens in zwei Quellen erfordert 
die Berechnung oder wenigstens die geniigende Ab- 
schatzung des Warmewiderstandes zwischen Zahn 
und Joch. 

Wird eine auf dem Umfang der Maschine symme- 
trische Kiihlung und Warmeerzeugung angenommen, 
so geniigt es, die Halfte einer Nut und eines Zahnes 
nach Bildia zu betrachten. Die Warmeabfuhr er- 
folet mit den Warmeiibergangszahlen 0,, 0 in den 
Luftspalt und mit # in den Statorkérper bzw. in 
die Kiithlkandle in demselben. Die Wa&arimefliisse 
sind dementsprechend 4,9, G39, Joo und auBerdem 4g, 
zwischen Zahn und Joch, bzw. q,3 und q33. Die 
Richtung der Warmefliisse qo, 913, 932 ist im all- 
gemeinen unbekannt, sie kbnnten auch inentgegenge- 
setzter Richtung, alsim Bild1a eingezeichnet, flieBen. —_Bild1. a Warmestrémung im System Zahn-Joch-Nutkupfer. 

‘ 2, = Zahn, z, = Joch, z; = Nutkupfer. b entsprechender 

Der entsprechende Ausschnitt aus dem ganzen Ausschnitt aus dem Warmequellen-Netz der Maschine. 
Warmequellen-Netz der Maschine ist im Bild 1b 
angedeutet mit den Quellen z, = Zahn, z= Joch, z; = Kupfer und mit den Warme- 
widerstanden 7,. usw. zwischen ihnen. Die Str6mungen und Warmewiderstande in der Axial- 
richtung der Maschine kénnen hier unberiicksichtigt bleiben. Sie gliedern sich einfach zu den 
vorhergehenden parallel, wie es im Prinzip der Warmequellen-Methode liegt: 2; ist also mit 
z, = Wickelkopf iiber einen inneren axialen Widerstand 7g, verbunden, zu 79 liegt parallel der 
‘der Abkithlung der Stirnflachen des Statoreisens entsprechende Widerstand 72, usw. 


2. Vereinfachung des Problems 


Bei Vernachlassigung der Axialstr6mung liegt ein zweidimensionales Problem vor. Fiir die 
Temperaturverteilung im Kupfer und im Eisen gelten zwei Differentialgleichungen (PoIsson- 
sche Gleichungen) mit zwei Veranderlichen, deren gemeinsame Lésung bei den komplizierten 
Grenzbedingungen und bei der ungleichen Warmeerzeugung im Zahn und Joch undurchfithr- 
bar ist. Die Methode der konformen Abbildung, auch bei Vernachlassigung der Neigung zwi- 
schen den Begrenzungslinien (Winkel y im Bild 1a) fiihrt auf praktisch unlosbare Integration, 
Eine Erganzung dieser Methode durch graphische Integration ist nur bei maBiger Anzahl von 
Grenzbedingungen anwendbar und hier unbrauchbar, ungeachtet dessen, da andere Schwierig- 
keiten infolge der Warmeerzeugung im abgebildeten Polygon entstehen. 
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Es wird sich zeigen, daB die Widerstandserhéhung, die durch die Querschnittserweiterung 
des Warmeflusses g,. verursacht wird, nur klein ist bei den Abmessungsverhaltnissen, die uns 
interessieren, und wir konnen, neben der Vernachlassigung des Neigungswinkels, noch weitere 
griindliche Vereinfachungen zulassen. Wir untersuchen nur den Warmewiderstand zwischen 
zwei Staben verschiedener Querschnitte ohne innere Warmeerzeugung und verlangern sie 
beiderseitig ins unendliche (Bild 2a, b, c). Die fiir c gefundene Widerstandserhohung, die durch 
Lésung einer Potentialstromung leicht zu erhalten ist, wollen wir dann naherungsweise auch 
in allen tatsachlich vorkommenden Fallen verwenden. 


ie —_— 
Cc 
Bild 2, Vereinfachung des Problems. a Neigungswinkel y zwi- Bild 3. Lésung der Potentialstr6mung. 


schen den Symmetrieachsen der Nut und des Zahnes. b Nei- 
gungswinkel nicht beriicksichtigt, c Verlangerung beider Stabe. 


3. Losung der Potentialstromung 


Die Lésung der einfachen Potentialstr6mung nach Bild 2c ist bekannt und in mehreren 
Lehrbiichern zu finden (z. B. [1] S. 358, [2] S. 125, [3] S. 262), so daB die Ableitung der End- 
formel hier nicht notwendig erscheint. Fiir den Zusammenhang zwischen der Grundebene 
z= x +72y und der Abbildungsebene v + 7 w gelten die Beziehungen 


ae ple h BH 
Lis bt+i1 4H p+ h/H i 


(1) 


allie a= (Hib, f= eer | 


G 
Die Lage des Koordinatensystems z und die Abmessungen h/, H sind im Bild 3 eingetragen. 
Beim Verhaltnis H/h = 2 verlaufen bei v = + 2 und v = — 2 die Potentiallinien schon prak- 
tisch gentigend genau geradlinig, so daB bei |v| > 2 nur mit reellen Variabeln gerechnet werden 
kann. 


Fir groBe positive v = v, ist € = e*™ sehr groB und man kann 


ah ae Se 
p=): Pia aoa) 


setzen. Der entsprechende Wert von z = x, ist, wie eine einfache Berechnung nach zulassiger 
Vereinfachung ergibt, 


tal a—t1 h 1+ h/H 
ea (in ; -f. ee He MT Ferae (6) (2) 
Fiir groBe negative v = v, ist € sehr klein, so daB . 


al ral tL 


XLY. Band 3 5 Z : 
Hert Bees J. Hax: Der Warmewiderstand zwischen Zahn und Joch 51 


Pind sures — 4, erhalt man 


| 
ECE iki. ede. a 


gh mee ( 1— h/H h 1 4 


= )—hy>o. (3) 
Der Warmewiderstand zwischen den beiden Staben von der Lange || und x, ist, da v, — 2 
den Temperaturunterschied darstellt und der Warmeflu8 = 1 (W) angenommen wurde, 


ae °C/W 
Pe gee Ooi) (4) FL 
wobei 2 die Warmeleitfahigkeit des Materials beider 
Stabe (W/m °C) und u(m) die zur Bildebene senkrechte 40P 
Abmessung der Querschnitte bedeutet. 
Setzt man in (4) die aus (2) und (3) berechneten ast 


Werte v, und v, ein, so ergibt sich 


_ 1 (lal iL 2 ae 
é Au ( ely ale ak R, (5) SS 
mit Of 
a See eo 1+4+h/H , H/h— h{H 
as cE ray} 1—h/H ' as 4 BAS 02 


Es ist also 7, = R,/Au die VergroBerung des Warme- 
widerstandes oder der Zusatzwiderstand zwischen den 0 02 04 06 0,8 4,0 
beiden Staben, der durch die Querschnittsanderung her- hli—— 

vorgerufen wird. Werte von R, nach der Gl. (6) sind fiir er gs ek tee ia a 
verschiedene Verhaltnisse /H im Bild 4 dargestellt. : 

Eine Nachrechnung ergibt, da8 bei H/h = 2 auch in der Nahe von |v| = 1, wo die Po- 
tentiallinien schon leicht gekriimmt sind (s. Bild 3), der Wert von R, nur unbedeutend von dem 
fiir groBe |v| berechneten Wert nach (6) abweicht. Die kleinsten fiir elektrische Maschinen in 
Betracht kommenden Verhaltnisse x/H liegen bei 1,2, wahrend H/h in der Nahe von 2 liegt. 
Nimmt man fiir den ungiinstigsten Fall x,/h = x,/H = 1,2 an, so betragt die perzentuale Ver- 
eroBerung durch den Zusatzwiderstand 


100 FR H 100k, 

1 5 mote ie Hh’ (%o) (7) 
x fz ar 2 

166 h !} 


also bei H/h = 2... (1/1,2) (100- 0,25/3) = 7%. Nur in Spezialfallen, z. B. fiir einen Rotor 
mit sehr schmalen Zahnen, kann das Verhaltnis H/h bis auf 4 ansteigen. Auch dann betragt 
die perzentuale VergréBerung nach (7) nur rund 12%. 

Dieser verhaltnismaBig kleine Einflu8 des Zusatzwiderstandes rechtfertigt die Annahme, 
da8 man praktisch mit R, nach Gl. (6) auch in allen tatsachlich vorkommenden Fallen rechnen 
kann. Bei innerer Warmeerzeugung in beiden Staben ergibt R, etwas zu hohe Werte von 7,, 
wenn der Zahn warmer als das Joch ist, da der breitere Querschnitt in Wirklichkeit besser aus- 
genutzt wird. Fiir gréBere Maschinen, wo die Verhaltnisse x,/h und x,/H gréBer als bei dem 
oben berechneten ungiinstigsten Fall sein kénnen, betragt die VergréBerung nur einige v. H. 
der ohne R, berechneten Werte. 


4. Gang der Berechnung 


In einer fritheren Arbeit [4] wurde gezeigt, daB die inneren Warmewiderstande, die in 
einem Warmequellen-Netz vorkommen kénnen, nur durch eine fortschreitende Naherung genau 
berechenbar sind, da sie namlich von den Temperaturen aller benachbarten Quellen abhangen. 
Zur praktischen Anwendung der Methode ist es deshalb erforderlich, sie durch Naherungsaus- 
driicke zu berechnen, welche man bei Vernachlassigung einiger Quellen des Netzes oder bei Ab- 

4* 
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schatzung einiger Temperaturen benachbarter Quellen erhalt. Auch der gesuchte Warmewider- 
stand zwischen Zahn und Joch ist ein innerer Widerstand, der ahnlich wie die inneren axialen 
Warmewiderstande im Kupfer zu berechnen ist. 

Wir betrachten nur den Fall einer rechteckigen Nut. Bei anderen Nutformen mu8 man 
sich mit einer Abschatzung der fiir die Berechnung von 7,3 und 7,3; maBgebenden Abmessungen 
begniigen. Fiir den Zahn wird mit einer mittleren konstanten spez. Warmeerzeugung 8, (W/m?) 
gerechnet. Die Beriicksichtigung der in Wirklichkeit 6rtlich veranderlichen Warmeerzeugung 
41 = f(x) ist umstandlich und scheint praktisch nicht notwendig zu sein. Wichtiger und sehr 
wiinschenswert ware dagegen die Bestimmung von 4, durch Messung. Bekanntlich miissen die 
berechneten Eisenverluste der Maschine, welche den an Proben ermittelten MeBwerten ent- 
sprechen, durch einen Korrektionsfaktor erhéht werden, der den Einflu8 der Bearbeitung dar- 
stellt. Diese Bearbeitung beeinfluBt in viel héherem MaBe die Zahne als das Joch, so daB 
zwei verschiedene Korrektionsfaktoren vorteilhaft waren, um den Einflu8 der Bearbeitung 
besser auszudriicken. 

Wird y = o angenommen, so liegen zwei Stabe nach Bild 5 vor, mit den Langen J, /, (m) 
mit den Querschnitten s, = hw, s, = H u (m?), mit der Warmeleitfahigkeit 2 (W/m ° C) und 
mit der spez. Warmeerzeugung 4, 3. (W/m'). Die Abszissen x sind separat in derselben Rich- 
tung genommen. Der Verlauf der beiden Temperaturen ¢,, f, wird mit Hilfe der Porssonschen 
Gleichung fiir eindimensionale Warmestrémung gefunden. Fiir den Ubergang zwischen den 


beiden Staben setzen wir 
By (a Se eee (8) 
1\ dy x=l, 2 \dx ‘=o ewe 


wobei g,. der vom Zahn ins Joch gehende WarmefluB (W) ist. Den Einflu8 des Zusatzwider- 
standes 7, konnen wir naherungsweise so darstellen, daB wir 

ty —leo = 7, Ue (9) 
setzen (Bild 5). Der tatséchliche Temperaturverlauf erfolgt nach einer Kurve k, so daB die 
Tangente zu ¢, bei x = 0 in Wirklichkeit etwas kleiner ist, als angenommen wurde. Man be- 


r=0 r=l,0=0 
ae -—e 
va 
Bild 5. Zwei Stabe, die naherungsweise den Zahn Bild 6. Zwei Stabe wie im Bild 5. Warmeabfuhr nur auf dem Endquer- 
und das Joch darstellen. schnitt des zweiten Stabes (Joch). Dieses Bild gilt auch fiir die Berech- 


nung bei seitlicher Warmeabfuhr vom ersten Stab (Zahn), s. Absatz 6, 


rechnet die beiden mittleren Temperaturen ¢,,, fg, der Stabe und der gesuchte Warmewiderstand 
ist dann 


Y. ieee OY OY 11) 
12 (10) 


Der Temperaturkoeffizient der Eisenverluste ist bekanntlich negativ [5] und kleiner als 
derjenige fiir Kupfer. Wir kénnen also mit konstanter Warmeerzeugung im Eisen rechnen. 
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5. Erste Naherung 


Wenn angenommen werden kann, daB die Abkiihlung iiberwiegend nur auf der Statoreisen- 
Oberflache stattfindet, daB also 0,, 0,3, 0 im Bild 5 vernachlassigt werden kénnen, so hat man 
zwei Stabe nach Bild 6 zu behandeln. Die Potssonsche Gleichung fiir die Zahntemperatur ¢, 
vereinfacht sich zu 


d?t 
und ist durch die Parabel 2° 
4 = a — GX" 
gelést. Fiir ¢, gilt dieselbe Differentialgleichung, es ist jedoch 
ty = Ay — by % — Cy X? (12) 


zu setzen. Die beiden Differentialgleichungen liefern dann c, = 3,/2/, cp = 39 /2A. Die Be- 
ziehung (8) ergibt 


ues 
oe sy A 
und aus der Gl. (9) erhalt man 
Ag = a, — 0, 1? —7, Oye (13) 


mit qd). = 3; 5,4. Der bei x = /, iiber s, abgehende FluB ist der gesamten Warmeerzeugung 
gleich, so daB 


81 Sip + 82 Spl = Dy Sa bi 
ist, wenn man die Umgebungstemperatur, in welche #, miindet, zur Bezugstemperatur der Be- 
rechnung wahlt (f; = 0). Aus dieser letzten Gleichung erhalt man einen Ausdruck fiir a,, der 
zusammen mit (13) den Koeffizienten a, oder die Maximaltemperatur am Zahnkopf 
1 S$ A if A 
a = hm = 81) 4 (ie )|4 ie, (2 bs] bY, 1 Sy 


Ss Ds Bs 


ergibt. ; 
Die mittleren Temperaturen beider Stabe sind 


i, = 4; ¢, 02/3 und %, = 4,—b,h/2— 0,223 


und der gesuchte Warmewiderstand zwischen Zahn (J,) und Joch (/,) nach Gl. (10) ist nach ent- 
sprechender Umstellung, wenn noch s,/s, = 4/H und s, = 4 uw geschrieben wird 


ee ed ee ee (14) 
BAS 2A Spee OVA Sy Bn HCN le BSE 34,24 


6. EinfluB der Wicklung 


Wenn wir den Einflu8 des Warmeiibergangs zwischen dem Zahn und der Wicklung beriick- 
sichtigen wollen, so muB die mittlere Nutkupfertemperatur ¢; (Bild 6) abgeschatzt und die 
Warmeiibergangszahl #,, in die Rechnung eingefiihrt werden. Der Wert von #3 ist aus der 
Isolationsdicke, der mittleren Warmeleitfahigkeit der Isolation und der Ubergangsflache be- 
rechenbar. 

Fiir die Zahntemperatur ¢, gilt statt (11) die Differentialgleichung 

dt 


AS; ae = Sj — & ys (4 — 43) 
und ihre Lésung (vgl. [4], [6]) ist 
t, = a, — b, Cos (@ x) 
mit 
Wars 


a oS * / : 
sane em ae @ Ne 
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Fiir das Joch gilt die Differentialgleichung (11) mit Index 2 statt 1 und ihre Lésung ist die 
Gl. (12) mit cy = 3,/2A. Die Beziehung (8) liefert 


bac Oy OD SiS (15) 
wenn zur Abkiirzung S$, = Sin (9 J,) geschrieben wird. Aus der Gl. (8) erhalt man gleichzeitig 
Gua =A 5 (dt, dx), 2, = 15; OVO Sie (16) 
Die Gl. (9) fiir den Ubergang liefert dann 
a, — 6, C, — a, = 7,15,6, 058, (17) 
mit C, = Cos (oJ,). Fiir den Ubergang iiber %, gilt 
— A sy (dt,/dx) 1, = Be So toy 
und daraus, mit Benutzung von (15), 
a, — b, 0 Sy (lg + A]Ds) Sy/sq = Coletta FH 2A/0e) = Ts (18) 


wobei J zur Abkiirzung. Die beiden Koeffizienten a, und 0, sind aus den Gl. (17) und (18) zu 
berechnen. 


Die Differenz der beiden mittleren Temperaturen ist 
te—ty, = a —b, Sifoly— a, + by l,/2 + ce 15/3 - 
Wird a, — a, aus der Gl. (17) berechnet, 6, nach (15) und cy = 3,/2 A eingesetzt, so folgt nach 
einfacher Umformung 
i Ne eay ee (c eed Ten yee ae Fk oS) +*4 
1s oe Sed Rar ~~ os 351 @ 1 note 


Der nach (10) berechnete Warmewiderstand zwischen Zahn und Joch ist dann 


ey : & | 


MONS Olbil  PLRSS 


ih Lah Coilt /)Sa@) sary as GHGs LANs. 
2 eps oF Cy/( 145) a 2 ofA) [Sz (19) 
1 


Mit der Reihenentwicklung C,/S, = 1/0 l, +. 1,/3 — 0° 2/45 + 2 0° 2/945 — 07 4/4725 +--- 
im ersten Glied hat man 


wobei 


aa Osis) z Ada ae 
R= Ol —2 oh + eh) + 


derselbe Korrektionsfaktor ist, den wir bei der Berechnung von inneren axialen Warmewider- 
standen gefunden haben [4], [6]. Das letzte Glied in (19) entspricht dem EinfluB der seitlichen 
Warmestré6mung zwischen Zahn und Nutkupfer. Wird dieselbe Reihenentwicklung auch im 
letzten Glied verwendet, so kann der Ausdruck fiir 7,. auf die Form 


n: Cae ae ly erie a 
f= GeO RGt2 +5. 422) +74042)4 fens b= )2 (20) 
uberfiihrt werden, wobei 
ye 


6 A (a,— T) 

bedeutet. 

Wenn ¢; > ¢,, ist, wenn also vom Nutkupfer ein Warmestrom in den Zahn abgeht, so ist 
a, — T gréBer als bei entgegengesetzter Stromrichtung. Der Warmewiderstand 7,, ist also bei 
tz > ts kleiner als bei f; << ,,. Eine ahnliche Auswirkung ist von dem vom Luftspalt in den 
Zahn oder umgekehrt flie8enden Warmestrom zu erwarteh. 

Wenn der Warmewiderstand 7, nach (19) oder (20) mit einer geschatzten Nutkupfer- 
temperatur fs, berechnet wurde, und wenn die Lésung des Warmequellen-Netzes ts > ts, ergibt, 
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so wurde mit einem zu hohen Werte von 7,2 gerechnet. Ergibt sich jedoch t; < ¢3,, So wurde mit 
zu kleinem 7,, gerechnet und die Lésung kann dann mit dem korrigierten Wert 7,. wiederholt 
werden, wobei es geniigt, in a, die durch die erste Berechnung gewonnene Nutkupfertemperatur 
einzusetzen. 

Der vom Nutkupfer direkt in das Joch abgehende WarmefluB q,, (Bild 1a) ist gew6hnlich 
viel kleiner als g,, und verkleinert nur maBig den nach (14) oder (19) berechneten Wert von 749. 
Fir den Warmewiderstand 7,5 zwischen Nutkupfer und Joch kann man also geniigend genau 

0 IE 


— = 21 
"28 (ss— Sy) ra 21s, Ca 


setzen, wobei 6 die Isolationsdicke (m) und /; die mittlere Warmeleitfahigkeit der Isolation 
sind. 

~Bemerkung. Der Einfachheit halber wurden alle Warmewiderstiande, insbesondere 745, 793 
usw. fiir eine halbe Nut und einen halben Zahn berechnet. Da bekanntlich die Warmequellen- 
Netzmethode entweder fiir die ganze Maschine oder fiir eine Nut verwendet werden kann, ist 
eine entsprechende Umrechnung notwendig. Wird die ganze Maschine behandelt, so sind die 
hier angegebenen Werte mit (2 x Nutenzahl) zu dividieren. 


7. EinfluB des Neigungswinkels 


Der Neigungswinkel y zwischen der Nut- und Zahn-Achse wird in erster Naherung dadurch 
beriicksichtigt, daB man fiir s, und s, die Mittelwerte der in Wirklichkeit veranderlichen Quer- 
schnitte einsetzt, also Werte die sich fiir x = 1,/2 bzw. x = 1,/2 ergeben. Dabei wird R, fiir das 
tatsachliche Querschnittsverhaltnis am Ubergang bestimmt. 

Fiir den einfachen im Absatz 5 behandelten Fall kann der Einflu8 der Neigung naher unter- 
sucht werden. Die Querschnitte beider Stabe seien s,9, S99 bei x = 0 und sj), sy, bei x = 1, bzw. 
4 = Igibild 7). bir beliébige-w ist also s,, =-s,9 +. € %} Sox = Soo 1-€ %, Welln 


€ = (S37 — Syo)/4, = (Sar — S20) /ls 
bezeichnet. Die mittleren Querschnitte sind 
Sis = (Sto + S12)/2 und sy, = (Sgq + Sq))/2. 


Wie eine einfache Ableitung zeigt, gilt fiir die Tem- 
peratur ¢, statt der Gl. (11) die Differentialgleichung 


d dt 
ioe (S152) a (22) 


-also 
ty \ E dt, fj 
dx | Sig + EX Ax 7 


und ihre Loésung ist 


= 0 (23) 


t = — By (Su + € x)? + oy In [Q, (Sto +e %*)]. (24) 
Die Einsetzung in (23) ergibt 
by = Ml4re (25) 
: Bild 7. EinfluB des Neigungswinkels zwischen der Nut- 
und die Grenzbedingung (d?,/dx),_, =o liefert c und Zahnachse. 


=—y Styl a AE 
Fiir die Temperatur ¢, im Joch gelten die Gl. (22) bis (25) mit dem Index 2 statt 1. Die 
Gl. (8) liefert nach einfacher Berechnung 


_ 32530 14 Sis 
2A& dé 


Co 


Die gesamte Warme muB tiber den Endquerschnitt s,, abgefiihrt werden, was die Gleichung 


81 Sts ly + bo Sos ly = Dy Sei tor = Dy Si [— Dp $31 + Cy In (Qp S21) ] (26) 
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ergibt. Fiir den Ubergang bei x = J, gilt die Gl. (9), also 
— by st) + cy Im (01 S11) + 2: S39 — CeIn (2 S29) = 7, br S154 + (27) 


Da die Koeffizienten 8,, b5, c,, cg bekannt sind, kann aus den Gl. (26) und (27) auch In (@, s,;) und 
In (0, Sy) bestimmt werden. Die durch Integration berechnete mittlere Temperatur im Stab 


L, ist 


eens a 
=: 2) eee = In (sy1/Sy9) — 4 + 4, In (0; 51) - 


t 
1s 31,8 ‘é 


Gleicher Ausdruck mit dem Index 2 gilt fiir die mittlere Temperatur ¢,, im Stab /, und nach 
Einsetzen aller Koeffizienten erhalt man fiir die Differenz 


ths — tes = 1%, Ne + Py gild + Po bold 
mit 
Py = —P2/3 — 3 by 8y9/4 € — Sig/2 €® + (Sio/2 4, €%) In (511/810) + (A Sis Sau/ts €*) 1m (Sgi/S20) | 
Py = [3/12 + Ip 89/4 € + S3q/2 €? — (S39 Se1/2 Fe €°) In (521/550) - 
Wie eine einfache Berechnung zeigt, sind die Verhaltnisse J, ¢/sy) und J, ¢/sy9 von der Nuten- 
zahl unabhangig und zwar ist 
Li C189 ==. 2 1] Loy PUN ele s3n = 2 a ae 


wenn D den inneren Stator-Durchmesser und D’ = D + 21, bezeichnet. Es ist also $,)/sy) = 
I + 04, S51/S99 = I + o2 und in den Polynomen P,, P, kann die Reihenentwicklung 


In (11/539) = In (4 + 0,) = 0, (1 — 04/2 + 09/3 +--+) 


und eine ahnliche fiir In(s,/sy9) verwendet werden. Der gesuchte Warmewiderstand 7,, ergibt 
sich dann nach Gl. (10), mit ¢,5 = 3, 4, s,, und nach entsprechender Umformung zu 
ly 


eo l, — 2 32 ly | 8 
We (4 RUST Gee (1 R,) 4 Cree (are Coates, (28) 


hi = 


wobei die den EinfluB der Neigung y darstellenden Korrektionsfaktoren sind 


3 1 1 1 
Ky = 6 o o2 + ---], 
it D Ht 4 5 Tages ia 


1 1 
Ry = @,(S S46 0, —-- 0 4 ; 
2 2 6 wees 2 & 20 
R 1 Iasicoy bale 7 
02 On Op eG 


Bei y = 0 geht der Ausdruck (28) in die Formel (14) iiber. Praktisch vorkommende héchste 
Werte von 2 1,/D bzw. 2 1,/D’ liegen bei 0,4 bis 0,5, so daB die Vernachlassigung des Neigungs- 
winkels y nur ausnahmsweise zu gréBeren Fehlern des fiir die mittleren Querschnitte s,,, Sos 
berechneten Widerstands 7,, fiihren kann. 


8. Auswirkung der Zerlegung von Eisenverlusten in Zahn- und Jochverluste 


Betrachtet man die drei im Bild 1b eingezeichneten Warmequellen z, = Zahn, ie —w Ochs 
2; = Nutkupfer, so sieht man, daB ohne Zerlegung nur Zwei Quellen (z, + 2) = Eisen und 
2; = Nutkupfer (Bild 8a) vorliegen, die durch einen Warmewiderstand 75, verbunden sind. 
Die Zerlegung der Eisenverluste ergibt dagegen das Schema Bild 8b. Da alle Verluste tiber den 
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Abkihlungswiderstand 7,, in beiden Fallen abgefiihrt werden miissen, wird die Eisentempera- 
tur /, im ersten Falle der Jochtemperatur é, im zweiten Falle gleich: 

ty = ty = Yao (2% + %2 + 2s) 5 
d. h. die Zerlegung hat keinen EnifluB8 auf diese Temperatur. 


Nach einem einfachen, allgemein giiltigen Satz zur direkten Reduktion in Warmequellen- 
Netzen kann das Netz Bild 8b in der Weise vereinfacht werden, da man die Verluste z, im 


Teal Wet M3 
t3 2,= Ta 
1 T2443 
Cc 


Bild 8. Auswirkung der Zerlegung von Eisenverlusten. a Warmequellen ohne Zerlegung, 
b Warmequellen bei Zerlegung, c Vereinfachtes Netz nach Elimination der Quelle Zahn. 


Verhiltnis 7,./7,, zu denQuellen z, und z, addiert und diese beiden verbleibenden Quellen noch 
mit dem Widerstande (7,, + 73) verbindet, wie es im Bild 8c angedentet ist. Die Nutkupfer- 
temperatur ohne Zerlegung ¢, ist nach Bild 8a 

Tere fac Ue 


wihrend bei Zerlegung sich nach Bild 8c 


Y. 1 i 
[pe eaaiey poy led bees Slay 12 )( 
: a (- “te + hs eae = 
ergibt. Da aber 7,, = (1/7,;3 + 1/re3)~4 ist, falls man wie tiblich nur die Warmewiderstande der 


Nutenisolation beriicksichtigt, und da ¢, = ¢, ist, betragt die Erhdhung der Nutkupfertempera- 
tur, welche durch die Zerlegung der Eisenverluste zum Vorschein kommt, 


—i 


1 45%. lig 
peat pooc ee IEE (: +2 s_). 2 
: M9 + %13 + 23 : : 131 V8 (29) 
Ohne Zerlegung ergibt die Berechnung immer zu niedrige Nutkupfertemperaturen. In man- 
chen Fallen kann die Zahntemperatur ¢, hoéher als die Nutkupfertemperatur #, liegen, wenn auch 
ohne Zerlegung eine Eisentemperatur ¢, < ¢, berechnet wurde. Aus den elementaren Glei- 
chungen, die fiir das einfache Netz Bild 8b gelten, findet man sofort, daB ¢, > ¢, bei 
es (30) 
2 ip 
ist. In diesem Falle ist also eine eingehendere Berechnung mit Zerlegung erforderlich. Aber 
auch bei ¢, < /, ist der Fehler (ts; — #3) nach Gl. (29) in vielen Fallen nicht unbedeutend. Diese 
Betrachtungen gelten nur zur Orientierung, da im vollstandigen Warmequellen-Netz der Ma- 
schine, wo sich zu den drei hier betrachteten Quellen andere Warmequellen und Widerstande 
reihen, die Endergebnisse nicht so einfach sind. 


Anwendungsbeispiel. In einer mittleren Maschine hat man im Stator die Verluste 
2 = 250, 2 = 380, 23 = 350 (W) und die Warmewiderstande 7,, = 0,06, 713 = 0,012, 723 = 0,04 
(°C/W). Es ist 2/23 = 250/350 = 0,715 und 793/715 = 0,04/0,06 = 0,667, so daB nach (30) die 
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Zahntemperatur ¢, etwas hdher als die Nutkupfertemperatur liegen wird. Der Fehler, der ohne 
Zerlegung des Eisens entstehen wiirde, betragt nach (29) 


tz — ty = 220-224 (250 + 350 ot) tieGz, 


0,052 


Zusammenfassung 


Bei zuldssiger Vereinfachung kann die Erhohung des Warmewiderstandes zwischen zwei 
Staiben mit ungleichen Querschnitten und ohne innere Warmeerzeugung mit Hilfe der be- 
kannten Lésung einer Potentialstromung gefunden werden. Diese Erhohung ist bei den in 
elektrischen Maschinen vorkommenden Verhiltnissen geniigend klein, so daB sie auch fiir Stabe 
mit innerer Warmeerzeugung und folglich fiir den Warmeiibergang zwischen Zahn und Joch 
naherungsweise verwendet werden kann. Die eigentliche Berechnung des inneren Warme- 
widerstandes zwischen Zahn und Joch wird mit Hilfe der fiir Warmequellen-Netzmethode all- 
gemein giiltigen Beziehungen fiir zwei Sonderfalle durchgefiihrt. Die Temperaturverteilung 
bei Beriicksichtigung des Neigungswinkels zwischen der Nut- und der Zahnachse wird fiir den 
einfachsten Fall ohne seitliche Warmeabfuhr aus dem Zahn berechnet, und es zeigt sich, daB 
dieser Winkel nur einen geringen EinfluB auf den fiir mittlere Querschnitte ermittelten Wert 
des Warmewiderstandes zwischen Zahn und Joch haben kann. Eine Orientierungsberechnung 
ergibt, daB man ohne Zerlegung der Eisenverluste immer eine zu niedrige Nutkupfertemperatur 
erhalt, und daB in vielen Fallen die Zahntemperatur hoher als die Nutkupfertemperatur liegen 
kann, ohne daB dies bei der Berechnung ohne Zerlegung zum Vorschein kommt. 
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Versuch einer genaueren Vorausberechnung 
des zeitlichen Erwarmungsverlaufes elektrischer Maschinen 
mittels Warmequellennetzen 
Von 
A. KESSLER, Bratislava 
Mit 9 Textabbildungen 
(Eingegangen am 5. Januar 1959) 

Ubersicht. Auf einer Untersuchung von MeBergebnissen (des Erwarmungsverlaufs) sowie der wichtigsten 
Berechnungsmethoden der Enderwarmung elektrischer Maschinen aufbauend, wird eine allgemeine angenaherte 
Loésung des Problems an Hand der Methode der Warmequellennetze aufgestellt. Es wird gezeigt, wie diese, 
praktischen Zwecken angemessen und unter einem MindestmaB an Rechenarbeit ausgewertet werden kann; 
die Ausfithrungen werden sowohl durch praktisch erreichte Ergebnisse, als auch durch eine vollkommen durch- 


gefiihrte Vorausberechnung des Erwarmungsverlaufs einer Maschine mit MeBergebnissen verglichen und einige 
Anmerkungen tiber die Anwendung von automatischen Rechengeraten erganzt. 


Ubersicht der wichtigsten Bezeichnungen 


t Zeit 

INC INCA) Temperatur: in (#) resp. im Punkte (¥, y, 2) 

T(%, t), T(x, y, z;t) desgleichen zur Zeit ¢ 

Wh, Ch ONC) mittlere Erwarmung: eines Kérpers, der Warmequelle (vy) zur Zeit ¢ 

5 Ge FAG) einspaltige Matrix der mittleren Erwarmung resp. Erwarmungsverlaufs (eines Systems 

von Warmequellen, Gl. (14)) 

p, O Zeitkonstante und Koeffizient der Exponentialreihe Gl. (1) 

w, W, W, W* Warmeentwicklung: spezifische, in einer Warmequelle, einspaltige Matrix Gl. (15) bzw. (5) 

Ais dae Warmeleitzahl, Warmeiibergangszahl 

Q, 4; q(t) Intensitat eines Kithlluftstromes, eines Warmestromes; zur Anwarmung notige Warme- 
leistung zur Zeit ¢ 

Taps Vip yye le Warmewiderstand: allgemein, tibergangs-, innerer, gesamt 

AN, aie Matrix (reziproke Widerstande) Gl. (13) resp. (5) 

(mc), (mc),, Cg Warmekapazitat: eines Kérpers, der Warmequelle (vy), Matrix Gl. (14) 

tig, Wy IBS, WV Koordinate, Lange, Umfang, Querschnitt, Oberflache, Volumen 

We Volumen der Warmequelle (v) 

a Koeffizient der Temperaturabhangigkeit der Warmeentwicklung 

E, Transformationsmatrix Gl. (10), Normalform von A* Gl, (16) 

Einleitung 


Der Frage nach der Bestimmung des zeitlichen Erwarmungsverlaufes + in elektrischen Ma- 
schinen wurde bereits eine nicht zu unterschatzende Anzahl von Arbeiten gewidmet, was jeden- 
falls zeigt, daB dieser Teilfrage des Erwarmungsproblems elektrischer Maschinen schon immer 
eine ziemliche Aufmerksamkeit geschenkt wurde. Heute ist ihre Aktualitat wenn méglich 
noch angewachsen, da die elektrische Maschine eine immer breitere Anwendung als Antrieb 
unter den verschiedensten Arbeitsbedingungen und Zyklen findet. Dabei geht es hier ebenso 
wie bei der Bestimmung der Enderwarmung im Grunde genommen um die richtige Bemessung 
und wirtschaftliche Ausnutzung der Maschine. 

Gehen wir die bisher veréffentlichten, einschlagigen Arbeiten kritisch durch, gewinnen wir 
notgedrungen den Eindruck, daB eine genauere Vorausberechnung des zeitlichen Tem- 
peraturverlaufes fast véllig auBerhalb jeder praktischen Méglichkeit liegt. Dies ist aber keines- 
wegs tiberraschend, da ja schon die genauere Vorausberechnung der Enderwarmung auf groBte 
Schwierigkeiten st6Bt [1]. Dabei ziehen wir schon in Betracht, daB einige Autoren [2], [3], [4] 
u. a.m. den meist einfachheitshalber vertretenen Standpunkt, die Maschine als einen einzigen 
K6rper zu betrachten, bereits haben fallen lassen. Die Frage der méglichen Auswertung dieser 
breiteren Ansatze wurde aber bisher noch nicht geniigend weitgehend verfolgt, trotzdem gerade 
die Beschrankung auf die sogenannte Einkérper-Lésung anscheinend das gré6Bte Hindernis fiir 
eine erfolgreiche Weiterentwicklung der Berechnungsmethoden darstellt. 


1 Weiterhin nur Erwarmungsverlauf im Gegensatz zur stationaren Enderwarmung — weiterhin nur End- 
erwarmung. 
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In den letzten Jahren wurde nun der Versuch unternommen, die sogenannte Methode der 
aquivalenten Warmequellennetze fiir die Berechnung der Enderwérmung elektrischer Maschi- 
nen weitgehend durchzuarbeiten [5], [6]2. Die verhaltnismaBig giinstigen Ergebnisse, die er- 
zielt wurden und vor allem die weitere Ausbaufahigkeit dieser Methode [1] gaben AnlaB, auf 
der gleichen Grundlage auch das Problem des Erwarmungsverlaufes in Angriff zu nehmen 
[7], [8]. 

In ee wollen wir nun versuchen einen kurzen Abri8 der an Hand dieser Methode 
aufgestellten allgemeinen approximativen Lésung des Erwarmungsverlaufes und deren An- 
passung an praktische Erfordernisse zu geben. Um ein méglichst vollstandiges Bild sowohl der 
erzielten Ergebnisse, als auch der praktischen Handhabung zu gewahren, wird auch der Er- 
warmungsverlauf eines kleinen Induktionsmotors bei Anwarmung, Abkiihlung sowie bei inter- 
mittierendem Betrieb praktisch berechnet und mit MeBergebnissen verglichen. Da die bei der 
Vorausberechnung zu leistende Rechenarbeit fiir die Praxis von ziemlichen Belang ist, sind im 
Nachtrag auch noch einige Hinweise betreffs der Anwendung von modernen Rechenautomaten 
resp. elektrischen Netzmodelen® beigefiigt. 


1. Einleitende Betrachtungen 


Bei der Untersuchung einer Reihe durch Messung ermittelter Anwaéarmungs- und Ab- 
kithlungskurven* [#(t)] hat es sich gezeigt, daB diese bei einer technisch annehmbaren Ge- 
nauigkeit durch Exponentialreihen mit konstanten Koeffizienten — wie dies bei Ausgleichs- 
vorgangen tiblich ist — 


6) = 0,,—F O,e-”", (1) 


von mindestens zwei bis drei Gliedern der Summe in Gl. (1) dargestellt werden miissen, sollen 
10—50% Abweichungen von dem tatsachlichen Verlauf vermieden werden. 


Fiir die verschiedenen Teile der Maschine ergeben sich dabei gleiche Zeitkonstanten (,) 
und verschiedene Konstanten bei den Exponentialen (@,). 


Dies zeigt, daB man von einer Losung des Erwarmungsverlaufes an Hand eines Einkorper- 
modeles nur im giinstigsten Falle mehr als eine grobe Annaherung erwarten kann; und es liegt 
nahe, in dem Systhem von Funktionen der Form (1), die den Erwarmungsverlauf zumin- 
dest angenahert darstellen, die Lésung eines Systems linearer, simultaner Differentialglei- 
chungen mit konstanten Koeffizienten zu suchen, wie sie sich bei Mehrkérpersysthemen er- 
geben [3]. 

In diesem Zusammenhang sei noch auf einen wichtigen Umstand hingewiesen. Die Lésung 
der Enderwarmung und die Lésung des Erwarmungsverlaufes haben einen festen mathe- 
matisch-physikalischen Zusammenhang, dem Rechnung getragen werden sollte. Die End- 
erwarmung miBte namlich, wie schon der Name besagt, unbeschadet einer bloBen Annahe- 
rung, fiir ¢-»oo mit dem Grenziibergang lim *#(t) der beziiglichen zeitabhangigen Lésung 


t> co 


identisch sein. Nun hat sich bei der Berechnung der Enderwarmung eine Reihe von ,,Mehr- 
kérpermethoden™ bewdahrt, sei es an Hand eines eindimensionalen Warmestrémungsbildes 
oder eines Warmequellennetzes; siehe z. B. [11], wo auch bei eindimensionalem Warme- 
strémungsbild schon bloBe Teile von Maschinen als zwei resp. drei Stabe, die thermisch gekop- 
pelt sind, betrachtet werden. Es liegt also auch aus diesem Grunde theoretisch nahe, die Lé- 
sung des Problems unbeschadet mathematischer Schwierigkeiten an Hand einer Mehrk6rper- 
methode zu suchen. 


2 AuBer [5], [6] siehe unter [9] und [10] einiges Schrifttum zur Anwendung der Methode der Warmequellen- 
netze zur Berechnung der Enderwarmung. ‘ 


5 Modele des Warmequellennetzes. 


* KessLER [8], oder LicutBanp [4] und unverdffentlichte Messungen des Verfassers an Gleichstrom- 
maschinen. 
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2. Untersuchung und Bewertung der zur Berechnung der Enderwarmung 
hauptsdchlich benutzten Berechnungsmethoden 


Wenn wir die stationire Warmeableitung aus der elektrischen Maschine, d.h. aus einem 
System mehrerer thermisch gekoppelter Kérper physikalisch exakt beschreiben wollen, miissen 
wir, da es sich um isotrope Warmeleiter handelt, ein System von Differentialgleichungen 
Poissonschen (w > 0) und Laplaceschen Typs (w = 0), 


RATA 2) == 2", (2) 


lésen. T(x, y, z) bezeichnet die Temperatur, A die Warmeleitzahl und w die spezifische Warme- 
entwicklung. Das System besteht aus je einer Gleichung fiir jeden wesentlichen und unter- 
schiedlichen Teil der Maschine; Statorkupfer, Isolation, Statoreisen usf. mit Randbedingungen, 
die jeweils fiir zwei oder mehrere von ihnen gemeinsam sind. Um das Problem der Warme- 
ableitung tiberhaupt praktisch lésenzu kénnen, ist es aber unumganglich nétig von 
einer einfacheren Formulation auszugehen. 

Dies kann z. B. bewerkstelligt werden, indem wir die Warmeleitungsgleichungen wie folgt 
vereinfachen; da in Folge der Rotationssymmetrie bei Maschinen ohne ausgepragte Pole die 
Temperaturverhaltnisse um jede Nut die gleichen sein sollten, betrachten wir bloB einen einzi- 
gen Statorausschnitt mit einer Nut. Die Kupferwicklung in der Nut, deren Temperatur schlieB- 
lich maBgebend ist, sehen wir als einen homogenen Kupferstab an, dessen Temperaturverlauf 
in derQuerrichtung wir in Anbetracht der Lange der Wicklung vernachlissigen. Desgleichen 
betrachten wir das Statorpaket als einen Eisenstab, in den der Kupferstab eingelegt ist, usw., 
usf. 


4 Richtung der Warmeabgabe @ Kiihlluftstrom @ Kupfer @ Eisen 


Bild 1. 


Wir gelangen auf diese Art zum sogenannten eindimensionalen Warmeleitungs- 
bild der Maschine Bild1, das grundsatzlich auch eventuell um Rotoreisen und Kupfer 
erweitert werden kénnte. Die exakte Lésung des Erwarmungsverlaufes ist aber bei eindimen- 
sionaler Warmestrémung schon fiir den Fall eines einfachen Stabes [12] auBerst kompliziert 
und das Ergebnis ware bereits bei diesem einfachen Ansatz technisch viel zu unpraktisch. Der 
Erwarmungsverlauf fiir zwei parallel liegende gekoppelte Stabe (Bild 1) wurde in analytischer 
Form tiberhaupt noch nicht gelést und abgesehen davon hat die eindimensionale Lésung noch 
gewisse andere Nachteile, wie sie z. B. J. Hax [1] S. 307, 308 aufzeigt. 

Eine weitaus giinstigere Vereinfachung des mathematischen Problems ergibt sich, wenn die 
Warmeleitzahl in den metallischen Warmequellen im Hinblick auf die der Isolation z. B. in 
allen Richtungen als praktisch iiber alle MaBe gro8 angesehen wird. Daraus folgt dann, daB 
die Temperatur der Warmequellen als raumlich konstant zu betrachten ist (dT/dx = 0, 
0T/dy = 0) und wir erhalten mit Hilfe der Randbedingungen von Gl.(2) ein System alge- 
braischer Gleichungen oder bildlich gesprochen ein Warmequellennetz?. 

Die mathematischen Schwierigkeiten waren damit zwar behoben (wie wir sehen werden, 
auch in bezug auf die Lésung des Erwarmungsverlaufs), aber es bleibt natiirlich dahingestellt, 
ob in Anbetracht einer so weitgehenden Vereinfachung und der damit verbundenen Fehler 
durch allzu ungenaue Formulierung die gewiinschten Ergebnisse iiberhaupt zu erzielen sind. 
Schon ein einfaches Beispiel, ein Stab, dessen zwei Abschnitte verschieden intensiv gekthlt 


5 Wir wollen hier aufmerksam machen, da auch eine modelmaBige Lésung auf dem Wege tiber das Warme- 
quellennetz bedeutend einfacher als ein direktes Modelieren ist. Vergleiche z. B. KESSLER [10] und [13]. 
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werden, zeigt ([7], S. 290), daB die véllige Vernachlassigung des inneren Warmegefalles der 
Warmequellen zu betrachtlichen Fehlern fithren kann. Um so entstehende Schwierigkeiten 
zu beheben, gehen wir von einer etwas modifizierten Festlegung des Warmequellennetzes aus. 
Aus der formalen Analogie des Fourierschen Warmeleitungsgesetzes und des Gesetzes der 
elektrischen Strémung g = —Agrad T, J = —o grad @ bestimmen wir ahnlich wie dies bei 
elektrischen Stromleitern geschieht, fiir die verschiedenen Warmeleiter in der elektrischen Ma- 
schine deren ,,Warmeleit- bzw. Warmeiibergangswiderstand“ sowie den ,,inneren Warme- 
widerstand‘‘ der Warmequellen. Wir gelangen so zu einem quasi Ohmschen Gesetz der Warme- 
leitung, das die gegenseitige Abhangigkeit von Warmewiderstand 7, Warmegefalle 49(= 3,—0, 
oder = # — %, das ist die mittlere Erwarmung der Warmequelle) und Warmestrom q festlegt ; 


(a ae ; (3) 


Dabei sind die Warmewiderstinde wenigstens in der ersten Annaherung nur von den geo- 
metrischen Abmessungen und von den physikalischen Konstanten der beziiglichen Materialien 
resp. Kithlluftstr6mungen abhangig. 


Dank der verhaltnismaBig einfachen Geometrie der aktiven Teile der elektrischen Maschine 
ist es méglich, alle wesentlichen Warmewiderstande in der Maschine wenigstens annahernd 
zu bestimmen. Auf diese Art kann dann die elektrische Maschine thermisch durch ein 4qui- 
valentes Warmequellennetz, das den Warmewegen der Maschine mehr oder minder gut 
entspricht, als Mehrkérpersystem dargestellt werden. Den aktiven Materialien entsprechen 
Warmequellen, die durch die einzelnen Warmewiderstande miteinander und mit den Warme- 
senken verbunden sind; als Warmesenken betrachten wir die Kiihlluftstrémung oder bei Kiih- 
lung durch natiirliche Konvektion, die umgebende Luft. Wenn wir den Wert der Warme- 
widerstande angeben kénnen, kénnen wir, sind die Warmeverluste = Ergiebigkeit der Warme- 
quellen bekannt, auch deren mittlere Temperaturen bestimmen. Dabei gelten fiir das Netz 
analoge Gesetze wie die Krrcuuorrschen bei elektrischen Netzen. Einzelheiten siehe z. B. bei 
Hak [5], [6] und KEsseEr [7], [8]. 


Die Hauptaufgabe und maBgebend fiir das Ergebnis dieser Methode ist vor Allem die Be- 
stimmung oder Abschatzung der Warmewiderstande und der Entwurf des ge- 
eigneten Warmequellennetzes; damit wurde die mathematische Schwierigkeit der L6- 
sung komplizierter und mehrdimensionaler Falle auf die Bestimmung der Widerstande tiber- 
fiihrt. In ungiinstigen Fallen kénnte dies zwar auf das Gleiche hinauslaufen, aber meist ist es 
moglich, durch theoretische Analyse, Analogie, durch Erfahrung und auf MeBergebnisse ge- 
stiitzt, den Wert der Widerstande mit einer, den Verhaltnissen angemessenen Genauigkeit, an- 
zugeben. Dafiir sprechen verschiedene praktische, mittels dieser Methode zufriedenstellend 
geléste Aufgaben [9], [10]; andererseits harrt aber natiirlich noch eine Reihe von Schwierig- 
keiten der Lésung [1]. Von Vorteil ist auch, daB die Anzahl der zu lésenden Gleichungen, wie 
schon im Zusammenhang mit der Méglichkeit des Modelierens betont wurde, im Verhaltnis zu 
anderen numerischen Lésungen gering ist. Unter Umstianden ist es sogar méglich, die Lésung 
in analytischer Form anzugeben. 


3. Ein Ansatz zur approximativen Lésung des nichtstationdren Problems 


Auf Grund der Erwagungen des letzten Abschnittes wollen wir voraussetzen, daB die End- 
erwdrmung, wenigstens theoretisch, im Allgemeinen mit zufriedenstellender Genauigkeit be- 
rechnet werden kann und daB sich also fiir jede Maschine ein angemessenes Warmequellennetz 
angeben 1aBt. 


Fassen wir die Enderwarmung der einzelnen Warmequellen @,, wo y die einzelnen Warme- 
‘ 


quellen bezeichnet, zu einer einspaltigen Matrix #, die reziproken Warmewiderstaénde gemaB 
ihrer Bezifferung zu einer quadratischen Matrix A und die Verluste, d. h. die Ergiebigkeit der 
einzelnen Warmequellen sowie die Temperaturen der einzelnen Warmesenken zu einer ge- 
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eigneten einspaltigen Matrix W zusammen (siehe Abschnitt 7 Gl. (11) oder [7]); dann kénnen 
wir das Problem der Bestimmung der Enderwarmung der Netzquellen folgend ausdriicken 

ASS W.. (3a) 
Jede einzelne Zeile der Matrizengleichung (3a) stellt die Warmebilanz der mit der Zeile gleich- 
lautend bezifferten Warmequelle dar, so daB die linke Seite von Gl. (3a) physikalisch als die 
aus den Warmequellen abgefiihrte Warmemenge bezeichnet werden kann. 

Wollen wir nun den Erwarmungsverlauf der Netzquellen ermitteln, so geniigt es, auBer der 
in jedem Augenblick aus den Warmequellen abgeleiteten Warme A @, noch die im zeitlich ver- 
anderlichen Falle zur Anwarmung der Warmequellen nétige Warmeleistung, gemaB den ihnen 
eigenen Warmekapazitaten (mc), zu beriicksichtigen. Bei einer einzelnen Warmequelle er- 
halten wir fiir die fiir die Anwarmung bendtigte Warmeleistung q,(¢), wenn wir mit T,(x, y, z; 2) 
die Erwaérmung im Punkte (x, y, z) zur Zeit ¢ bezeichnen und weiter @ die Dichte und c die 
spezifische Warme ist 


dq, (¢) =eef{ [tz T(x, 9,251) a dx dydz=o-c He LGN et Gied 2% 
V, i 


Fihrt man den mittleren Erwarmungsverlauf #,(t) im Kérper durch Integration tiber sein 
ganzes Volumen V, ein #,(t) = — i T,(%, y,2;t) dV ist, da wir die Reihenfolge der Inte- 
a 


v 
gration und Derivation vertauschen kénnen weiters 


dg,() =9-¢ = 1] T,(%, y, 23) ap df= ne Veal 
dt at 
Wes 
und es kann allgemein fiir den Erwarmungsverlauf der Netzquellen geschrieben werden 
dat) fa = 
Co + Ady) = W, 


oder nach Multiplikation von links mit der zu Cy reziproken Matrix Co*, wobei Cg* A = A* und 
Col. W = W* ist; ie 
d9(t) 
cay 
C, ist die Matrix der Warmekapazitaten des Netzes, deren Determinant von Null verschieden 
ist. Gl. (4) ist, wie erwartet, ein System linearer Differentialgleichungen mit konstanten 
Koeffizienten und die Gleichung geniigt auch der Bedingung 


= — A* ot) + W*. (4) 


lim 0(t) = 0 

t— oo 
Fiir ¢ > oo ist, sofern Gl.(4) und Gl. (3) eine endliche Lésung haben, lim d&(é)/dt = 0, so daB 
Gl. (4) in Gl. (3) tibergeht. eae 


Es muB allerdings betont werden, daB die Lésung der Gl. (4), d. h. der Erwarmungsverlauf 
des aquivalenten Netzes die tatsachlichen Verhialtnisse nicht voll widerzugeben vermag. 


H 
Tes" 


U 


S53 


Bt) % 
(mec)( W 


Bild 2. 
Zeigen wir dies statt an einem vdllig allgemeinen Fall z. B. an dem konkreten Beispiel eines an 


der Oberflache gektihlten Stabes Bild 2. Die zeitliche Warmeleitungsgleichung lautet hier, bei 
eindimensionaler Warmestrémung (KEssLER [12}) 


aroma 6,209 tH. Tx) =— mF fi ta Teed}. (5) 
x 
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wenn die Warmeentwicklung im Stab mit « von T(x ¢) abhangt; die Anfangs- und Randbe- 


dingungen seien T(x 0) = 0; ay ° =o; T(xl) =o. Fihren wir die Integration der 


O* |x= 
Gl. (5) nach x von o bis / durch, erhalten wir mit Hilfe der Randbedingungen 
Se _ Fim) P+ SHA) = Wf +0 A}. (6) 


Dabei ist (mc) =o cV; Wy =m-V; S=U-l, und die Anfangsbedingung #(0) = 0. Nun 
bedeutet der erste Ausdruck Gl. (6) physikalisch den in der X-Richtung abgeleiteten Warme- 


strom (q), der laut Definition Gl. (3) im stationaren Falle g = —A Fe pepe 


dx x=] Vin 
zeichnet die mittlere Temperatur des Stabes, 7,,, dessen (inneren) Warmewiderstand in der Langs- 


richtung und es gilt sichtlich rae = ligt ae alee Daraus ergibt sich ohne weiteres, 
£=1 {> ~ \e=1 

daB Gl. (6) durch Einsetzen von &(é)/7;,, und 7, = 1/S H(siehe Bild 2) der Gl. (4) fiir das ein- 

fache Netz Abb. 2 entsprache, doch bliebe festzustellen, ob der Gesamtwiderstand R—1/R 


= 1/r;, +1/r, — wie in Gl. (4) von allen Widerstanden schweigend vorausgesetzt wurde, 
auch wirklich konstant und fiir jedes ¢ dem, fiir den Limes ¢ > oo definierten Widerstand, 
gleich ist. 


Es zeigt sich, daB im Allgemeinen zwar der Wert der Warmewiderstande nicht konstant sein 
muB (siehe z. B. Nachtrag, Abschn. 10), aber doch in erster Annaherung als konstant betrach- 
tet werden kann. Wir kénnen daher sagen, daB Gl. (4), d. h. der zeitliche Erwarmungsverlauf 
des aquivalenten Netzes den tatsdchlichen Erwarmungsverlauf angendhert darstellt und 
zwar,wie leicht einzusehen ist, um so besser, je weniger die einzelnen Warmewiderstande 
y = y(t) von limr(t) = 7 = konst. abweichen. Als den tatsachlichen Erwaérmungsverlauf be- 


t— oo 


trachten wir dabei die Lésung eines Systems zeitanhangiger Warmeleitungsgleichungen der 
Form 


; F) 
AT(x,9, 231) —2= - 5 T(x, y, 3) =— , (7) 


das die thermischen Verhaltnisse in der Maschine physikalisch exakt beschreibt. DaB eventuell 
schon die Werte der Warmewiderstande in bezug auf die Enderwarmung nur angenidhert be- 
kannt sind, wollen wir hier auBer acht lassen. 

In bezug auf den zeitlichen Verlauf ist die Annaherung um so besser, einen je gréBeren Wert 
¢ wir ins Auge fassen und dann auch noch fiir ganz kleine Werte von ¢ (siehe Nachtrag). 


4, Die allSemeine Lésung des Erwarmungsverlaufes 
des 4quivalenten Warmequellennetzes 


Um die Lésung von Gl. (4) in geeigneter Form durchfiihren zu kénnen, bestimmen wir eine 
geeignete Matrix c durch die die Matrix — A* auf ihre Normalform n iiberfiihrt wird 


c1.— A*.c=n. (8) 


Setzen wir nun in Gl. (4) ot) =c¢ T(t) ein und multiplizieren wir von links mit c7, ergibt sich, 
wenn wir noch c-! W* mit M bezeichnen, wie leicht nachzurechnen ist 


=n-Ti)+M. 


Bei dem Ansatz T(t) = Ty e*! erhalten wir dann mittels Variation der Konstanten als Lésung 
von Gl. (4) wofern nur |A*| ¥ 0 ist 


Hl) = (e fi et. Mat +. Tol. (0) 


DaB Gl. (9) wirklich eine weitgehende allgemeine Lésung darstellt, ist leicht zu zeigen. Erstens 
ergibt Gl. (g) eine Lésung sowohl fiir Anwérmung, als auch fiir Abkithlung, denn sie gilt fiir 
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jegliches W*, also auch fiir W* = o und fiir beliebige Anfangsbedingungen #(0). Auch bleibt 
die Form der Gl. (9) unverandert, wenn W* = W*(¢) eine Funktion der Zeit ist. Es ist dann M 


durch M(t) = c¢1W(t) zu ersetzen. Anderungen der Betriebsverhiltnisse z. B. bei inter- 
mittierendem Betrieb, kénnen durch zwei Losungen, fiir die Anwaérmung und fiir die Abkiithlung 
durch Einsetzen der beziiglichen Anfangswerte in jeder Periode fiir beide Vorginge dargestellt 
werden usw. Nur ist zu fordern, daB die Elemente von A* wenigstens annahernd konstant 
sind, was wir ja vorausgesetzt haben. 

Bedenken wir weiter, daB das Integralzeichen vor der Matrix nichts anderes bedeutet, 
als daB jedes Element von e—"* von o bis ¢ zu integrieren ist, sehen wir, daB nach Durchfithrung 
der Integration und Multiplikation mit e"’ Gl. (9) die Form von Gl. (1) annimmt. Dabei wird 
die Anzahl der Glieder der Reihe durch die Anzahl der Gleichungen des Systems Gl. (4), d.h. 
durch die Anzahl der Warmequellen bestimmt. Untersuchen wir diese Form-Gleichheit etwas 
ndher. 


Q Korper ohne Warmeentwicklung 


o\ 70 Warmewiderstand 


@ Netzknotenpunkt ohne Warmekapazitat 


Q  Kihlluftstrom 


Bild 3. Bild 4. 
1 GuBeisengehause, 2 Statorbleche, 3 Wicklung in der Nut, 4, 5 Wicklungskopfe ,6 Rotoreisen, 7 Kupfer- 1 Statoreisen und Gehause, 2 Stator- 
wicklung oder KurzschluBkafig, 8, 9 Wicklungskopfe oder KurzschluBringe mit Kihlschaufeln, 10, 11 Luft wicklung, 3 Rotoreisen und Kupfer, 
zwischen Lagerschild und Wicklungskopfen. 4 Luft zwischen Wicklungsképfen 


und Lagerschild. 


Gingen wir bei der praktischen Bestimmung des Erwarmungsverlaufes z. B. bei Induktions- 
motoren von dem der Berechnung der Enderwarmung zugrundeliegendem Netzbild (Bild 3) 
aus, ergabe sich als Lésung eine aus neun Gliedern bestehende Reihe. Vergleichen wir dieses 
Ergebnis mit MeBergebnissen, die auBer O,, hochstens drei Glieder der Reihe Gl. (1) ergeben®, 
ist, da beide Ergebnisse den gleichen Erwarmungsverlauf mit nicht allzu unterschiedlicher 
Genauigkeit approximieren, klar, daB in diesem Falle in Gl. (g) die Glieder vom dritten resp. 
vierten Glied an nur einen technisch belanglosen Beitrag zum Losungsverlauf darstellen kénnen. 
Es sollte darum auch ein Gleichungssystem — wir wollen es ein ,Ersatzgleichungssystem‘ 
nennen — von zwei bis drei Differentialgleichungen geben, dessen Lésung gerade die einfachere, 


durch Messung ermittelte Reihe lefert. Denn ist I(t) in der Form Gl.(1) gegeben, kann man 


zweifelsohne aus O(t) und d@/dé eine Differentialgleichung der Form Gl. (4) bilden, sofern die 
Determinainte der Koeffizienten 9, nicht Null ist. 

Da zw schen den beiden Lésungen, im technischen Ma8stab gesehen, kein allzu wesent- 
licher Unterschied besteht, kénnen wir desgleichen von den beziiglichen Differentialgleichungen 
erwarten. Demgema8 kénnen auch die den Gleichungen zugrundeliegenden Warmequellen- 
netze in ihren Eigenschaften nicht wesentlich von einander abweichen, doch miiBte das Ersatz- 
netz aus entsprechend weniger Warmequellen bestehen. 


® Siehe Abschn. 1 oder [8]. 


on 
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Haben wir also auch in Gl. (9) eine Lésung des Erwarmungsverlaufes gefunden, die eine 
annehmbare Annidherung der exakten Lésung liefert, so miissen wir doch noch, um unser eigent- 
liches Ziel zu erreichen, vornehmlich aus praktischen Griinden, einen Weg zur Vereinfachung’ 
des der Enderwarmung zugrundeliegenden Netzbildes anstreben, um ein passend vereinfachtes 
Warmequellennetz, daB dem Ersatznetzbild entsprache, zu gewinnen. 


5. Der Entwurf eines vereinfachten Warmequellennetzes 
zur praktischen Bestimmung des Erwaérmungsverlaufes 


Bei der praktischen Suche nach einem geeigneten vereinfachten Warmequellennetz wollen 
wir wieder von dem fiir die Enderwarmung maBgebenden Netzbild ausgehen, da dieses ja mit 
dem gesuchten ,,verwandt sein muB und als bereits festliegend und bestimmt angesehen 
werden kann. Beziiglich seiner Bestimmung siehe insbesondere Hak [1] (Auf Probleme, die 
mit der Berechnung der Enderwirmung zusammenhangen, wollen wir hier nicht weiter ein- 
gehen). Dabei wollen wir festhalten, daB sich die ,, Verwandtschaft* nattirlich vor allem darin 
manifestieren muB, daB die mittleren Enderwarmungen der Warmequellen des vereinfachten 
Warmequellennetzes sich gegeniiber den urspriinglichen Werten nicht wesentlich andern 
diirfen. 

Ein Anhaltspunkt fiir den praktisch einzuschlagenden Weg bietet sich nun an Hand fol- 
gender Erwagungen: Ist in einem Warmequellennetz der Wert eines Widerstandes, der zwei 
Warmequellen verbindet, verhaltnismaBig gering und eventuell die Warmekapazitat einer der 
Warmequellen im Verhaltnis zur anderen ebenfalls klein, dann andern sich die Parameter des 
Netzes sicher nur wenig, wenn wir die beiden Warmequellen zu einer einzigen zusammen- 
fassen®. Ist weiter der Wert eines Widerstandes wiederum im Verhdltnis sehr groB, andert 
sich praktisch nichts, wenn der Widerstand weggelassen wird. Treffen die besagten Umstande 
bei einem Warmequellennetz zu, bietet sich die Méglichkeit, das volle Netz zu vereinfachen. 
Da in Wirklichkeit natiirlich das Verhaltnis zu allen Widerstanden eine Rolle spielt usw., lohnt 
es aber nicht, auf mathematischem Wege quantitative Regeln oder MaBstabe zu suchen, auBer 
der Kontrolle der Endtemperaturen bei dem versuchsweise vereinfachten Netz. Auf diese 
Art kénnten wir unter giinstigen Umstanden offensichtlich zu einem vereinfachten Ersatznetz 
gelangen, dessen Erwarmungsverlauf nicht allzusehr von dem des urspriinglichen Netzes ab- 
weicht und das, wie gefordert, in dem betrachteten Fall aus zwei bis drei oder héchstens vier 
Quellen bestiinde. Praktische Untersuchungen in dieser Richtung ergaben [8], daB bei einer 
Vereinfachung wirklich die gréBeren Zeitkonstanten verschwinden (Zahlentafel 1), daB aber 
insbesondere die letzte und auch eventuell die vorletzte Zeitkonstante bei dem Ersatznetz 
unter Umstanden betrachtlich von dem urspriinglichen Werte abweicht. Dabei wurde der 
Forderung nach Einhaltung entsprechender Enderwarmungen eventuell auch durch eine 
Korrektion der urspriinglichen Werte der Widerstande Rechnung getragen. Es 
scheint, daB zur Berechnung von zwei bis drei Messungen entsprechenden Zeitkonstanten, ein 
Ersatznetz von mindestens drei oder vier Warmequellen notig sein wird. 


Zahlentafel 1. Ps 
Anzahl der Quellen | 
im Ersatznetzbild : a ie WW 
Zeitkonstanten 0,0482 0,0410 0,0395 0,0414 
0,231 0,245 0,251 
1,19 
2,61 2,60 


? Gehen wir den entgegengesetzten Weg und versuchen wir das Ergebnis Gl. (9) zu verbessern, indem wir 
das Netz verfeinern, so gelangen wir schlieBlich zu einem Punkt, wo Gl. (4) in ein aus Gl. (7) abgeleitetes Dif- 
ferenzengleichungssystem tbergeht. 

8 Die Anderung der Zeitkonstanten in einem Netz mit zwei Warmequellen verglichen mit dem Wert der 
Zeitkonstanten bei zusammengelegten Quellen ist in [8] ganz allgemein berechnet worden. Es ergab sich, daB 
von einem gewissen Verhaltnis des Zahlenwertes des die beiden Quellen verbindenden Widerstandes zu den 
anderen usf. an, sich wirklich praktisch nichts andert, wenn man die beiden Quellen zusammenfaBt. 
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6. Praktische Ergebnisse 


Entsprechend den Ausfiihrungen in Abschnitt 5 wurde z. B. fiir die Berechnung der Zeit- 
konstanten einer Reihe kleinerer, véllig gekapselter _Induktionsmotoren folgendes, in Bild 4° 
dargestelltes Ersatznetzbild benutzt. Das volle, fiir 


die Berechnung der Enderwarmung benutzte Netz- 40 + gemessen 


bild” zeigt Bild 3. Seas Sheree 
Die an Hand des angefiihrten Ersatznetzbildes 4 aS a rer 

mit nur drei Warmequellen durchgefiihrten Berech- 5 oe 

nungen ergaben die in Bild5 mit den gemessenen 0+ 8. NSN 

Werten verglichenen Werte der Zeitkonstanten. Als ogL® “Se 

Kontrolle fiir die Giite der Abschatzung der Warme- 03 

widerstande wurde die berechnete Enderwarmung . PET era ft Sa 


des Ersatznetzes genommen. Es wurde festgestellt, 
daB Abweichungen von mehreren Graden, also Ab- 
weichungen der Werte der Widerstainde, denn die 
Verluste, Warmekapazitaten usw. sind ja relativ gut 
bekannt, schon zu betrachtlichen Abweichungen von 
den durch Messung ermittelten Werten der Zeitkon- 
stanten fiihrten (siehe Zahlentafel 2). Bei der Be- 
urteilung der angefiihrten Zahlen ist zu bedenken, 
daB die letzte und gréBte Zeitkonstante (hier 0,230) Bild s. 

aus MeBergebnissen bestimmt wurde und zwar an Zeitkonstanten bei Anwarmung. 

Hand einer graphischen Mittelung in semilogarith- 

mischem MaBstab. Die Bestimmung weiterer Koeffizienten wurde bloB mit Riicksicht auf 
die Genauigkeit der MeBergebnisse abgebrochen. Fiir den in Zahlentafel 2 angefiihrten Fall 
z. B. ware der dritte Koeffizient 0, bereits bedeutend kleiner als der vorausgesetzte MeBfehler. 
Die weitere Aufspaltung von @, und #, wiirde ein kleineres f, und ein einigemal gréBeres pz je 
nach den Umstanden ergeben, was mit den berechneten Ergebnissen gut tibereinstimmen wiirde. 


Zahlentafel 2. 


cu stat PFe stat 9Cu rot pi Po Ps 

gemessen 65,0 80,0" 83,0 — 0,0269 | — 0,230 — 
60,5 | 45,7 | 77,8** | —0,0276 =O; 102 — 0,745 
berechnet 68,0 | 58,3 | 80,7 | -—0,0193 — 0,0975 — 0,740 
86,0 69,0 96,0 — 0,0153 = 0,09 CONN 103.7150 


* Erwdrmung der Oberflache. Mittlere Erwarmung ca. 5—6° C héher zu nehmen. 
** Mittlere Erwarmung des Rotors. 


Es sei noch hinzugefiigt, da8 ein auf nur vier Warmequellen vereinfachtes Ersatznetz, be- 
stehend aus Stator und Rotor, Eisen und Kupfer auch richtig drei annehmbare Werte der Zeit- 
konstanten lieferte (f,, #2. und 3), wohingegen eine etwas gewaltsame Vereinfachung auf zwei 
Warmequellen keinen einzigen annehmbaren Wert ergab. An Hand des Ersatznetzbildes von 
drei Warmequellen konnten auch je zwei Zeitkonstanten bei Abktihlung, Uber- oder Unter- 
belastung der Maschine usf. mit relativ gutem Ergebnis berechnet werden. 


7. Vorausberechnung des zeitlichen Erwarmungsverlaufs eines kleinen 
Induktionsmotors bei Anwaéarmung und Abkihlung 


Es soll der Erwarmungsverlauf des vollig gekapselten 4,5 kW Induktionsmotors aus Zahlen- 
tafel 1, [8] berechnet werden. Als Ausgangspunkt soll uns das Warmequellennetz Bild 3 sowie 
das Ersatznetz Bild 4 dienen. 


® Bild 3, 4, 5 und Zahlentafel 1 sowie 2 sind [8] entnommen. 
10 Vergleiche auch Hax [1] Bild 6, 7. 
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Werden nimlich die Zahlenwerte der Warmewiderstande, die Warmekapazitaten usf. fiir 
das Netz Bild 3 bestimmt, zeigt sich, da8 der innere Warmewiderstand der Statorwicklung 
und des KurzschluBkafigs des Laufers einen bedeutend kleineren Wert als die anderen Warme- 
widerstiinde haben. Auch der Warmewiderstand zwischen dem Statorblechpaket und dem 
GuBgehduse, sowie zwischen KurzschluBkafig und Rotorblechen ist, wenn auch nicht aus- 
gesprochen, so doch merklich kleiner, als die tibrigen Widerstande. Wenn wir den Wert der 
genannten Warmewiderstande vernachlassigen und weiters die Warmeiibergangswiderstande 
von den Wickelképfen resp. von den mit Kiihlschaufeln versehenen KurzschluBringen zu- 
sammenfassen, erhalten wir das Ersatznetzbild Bild 4. 

Da die an Hand dieses Netzbildes berechnete Enderwarmung bereits merklich sowohl von 
den Werten, die sich fiir das Netzbild Bild 3 ergeben als auch von der durch Messung ermittel- 
ten Wirklichkeit abweicht, wurde der Wert des Widerstandes 7,, und 7,3 (Bild 4) korrigiert, 
um die Vernachlassigung der Widerstande 7,. und 7,7, (Bild 3) aufzuwiegen. Fiir das Ersatz- 
bild Bild 4 wurden so folgende Werte bestimmt: 


Tap == 050203 W,(Fe) = 160 | Schatzungswerte: 
fan = 070405 Wo2(Custat) = 2807 W Oe 160. 
77052200 °C/W W,,(Curot) = 150 Oy = ALG ‘ee 
i100 250 (HC) 2224 10% ies 
194 = 0,3120 (mc), = 2,47 + 10° a ; x = 3,93 
14, = 0,6700 (HC) eee Oe "Og Or 
Berechnen wir nun die Enderwaérmung konkret aus A § = W, dh. 
(-+-+2) cae Ee at 0 \ 
Uses Uap} (as 12 Tals 
a (ea Moa 0 es 
12 12 24 "94 
ae 0 ie 7 — 6 Wos] : : 
"13 13 "34 . a 
‘ pat $3: hip ( eee a 1 ) 
Yea ¥3 4 Vaba piisiay ania 
oq 
(%) (m+2s 
. dD, Woe 
zs 11 
d| | Wos 
0, 
O, o + — ) 
Yad 
Dazu setzen wir die Zahlenwerte (10) in (11) ein; wir erhalten 
O40 15221555 145 0 228,4 
— 21, 28, ) — 8,0 — 280, 
A= 2 ‘ ; W= 2 (12) 
— 4,5 O ies re 150,0 
) —8,0 —3,2 12577 6,7 


und wenn wir das Gleichungssytem auflésen, ergeben sich die folgenden, annehmbaren Werte 


; a, Dy D5 04 
Derechnet "47/0 =32.7. 1 6475 Gee oO 
gemessen 18,0 33,5 45—50 — 


11 Die Temperatur der Umgebung und der in die Maschine eintretenden Luft ist 20°C. Auf diese Tem- 
peratur bezieht sich auch W),, Wo, und a. 
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Um den Erwarmungsverlauf zu berechnen, stellen wir nun Gl. (4) zusammen. Wir erhalten 
nur drei Differentialgleichungen, da dem Knotenpunkt ,,4“ praktisch keine Warmekapazitat 
zugeschrieben wird (Luft), mit einer Nebenbedingung, die durch Umformung, wie folgt, 
explizit aus (12), vierte Zeile, ausgedriickt werden kann. 


H, = 0,531 + 0,630 3% + 0,252 Bs. (13) 
Nach Einsetzen von (13) fiir #, in (12) haben wir 
64,0 —21,5 —4,5 228,4 2 DA Ah Oe fo) 
Neen 215 9 28535 2,07 or We 28d Wee Cee O 2,47 +108 fo) ; 
See ye 0, O25, 151,7 f°) ) S717 10: 


was dann weiter 


2,857:107-% —9,598-10-* —2,009:- 1074 1,020+ 107 
A* =| — 8,704- 10-3 9,457:10-> —8 178-1074]; Wee yt sts .0>+ 
—8,704-10-4 —3,907-1074 NBA ROU e St 220 3A AO a 

ergibt. 

Die Eigenwerte der Matrix — A* ermitteln wir am einfachsten, indem wir die charak- 
teristische Determinante |— A* + # E| ausmultiplizieren und die Wurzeln des Polynoms 
z. B. an Hand der GRAEFFEschen Methode bestimmen 

D(p) = |— A* + pE| = 63 —1,3531- 10-2 £2 + 3,3156- 10->  — 1,8782- 10% =o 
Die = 1, 0AZ Oe Pp, = 2, 100 0, ps == © B00 1 Ons 


oder in Minuten ausgedriickt 
py = — 0,623, f,=0,130 und £, = —0,0489. 


Nun bestimmen wir noch die Matrix c, die die Matrix — A* auf die kanonische Normalform 
transformiert, indem wir (8) umformen 


Cra yo Crs (i “2 “Ss sal, O42 110m 0 ; 0 
AP C5 00 aCe Ni aay . O — 2,169 + 1073 fo) 
O 0) — 8,306 - 1074; 


Die sich durch Ausmultiplizieren ergebenden drei mal drei homogenen Gleichungen fiir c,, 
lésen wir, wobei wir willkiirlich c,, = 1, cg = 1 undc,3 = 1 bestimmen. Wir erhalten 


1 1 1 /0,1200 —0,1103 0,0080 
c = | —7,9310 0,9923 1,3515 | und hernach c7} = [ossis 0,0755 —0,2069 
0,2421 —1,3203 3,6298 0,2280 0,0348 0,1988 


Fihren wir nun M =c71. W* in Gl. (9) ein und fiihren wir die Matrizenoperationen, die 
die Gleichung vorschreibt, aus, erhalten wir vorerst 
; = 1,07 85 (crt 2) 
Tent (eee! aide = |e a aiee (tee) 


ie) 


14,6390 (1 ee 8,306 - 10™*t) 


und schlieBlich aus 0(t) = c- T(t) als Lésung des zeitlichen Erwarmungsverlaufs 
8, (t) i577 ae 1,08 bas 1,048- 10°72 4,21 +6 2160 JOne ts ae 14,64 | eae 8,306 +1074 ; 
D,(t) = 32,52 —e 8,55 AV 1,048 - 1057 % = 4,18 Cae 2, L6G 10°F 19,78 p e 8:306 - 10 4t ; 
0,(t) == 4731 =e 0,26 Sem 1,048 + 1072 ¢ . 5,56 ens 2,169-107°% 53,14: am 8,306 +10 4¢ - 
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Fassen wir nun den Abkiihlungsvorgang ins Auge, so ist zu bedenken, daB die Kihlluft- 
stromung entfallt, so daB sich die meisten Warmewiderstande bedeutend vergr6Bern 


ya = U3333% 24 = 0,357 5 
114 = 0,046, Fo5 = 1,000, 


Wy = Wo Wee 0; 
b= 4, = 0. 
73°04 70s Ty, = 2,500, 


Bestimmen wir bei gleichem Rechnungsvorgang wie bei der Berechnung des Erwarmungs- 
vorganges die Eigenwerte der Matrix — A*, erhalten wir 
pd, = —1,002-10%, f,=—6,782-10-* und , = — 1,0863- 104 
sowle 
1 1 Don 
c =| —9,2021 0,8840 1,0059 
—0,0766 —3,6182 1,2094 


In Gl. (9) ist M =, wohingegen die Anfangstemperaturen gegeben sind 


a 17577 a 
0(0)} = 132,525] = ¢c- 7(o). 


47,31 


Lésen wir diese Gleichung nach den unbekannten Elementen der einspaltigen Matrix T(o) 


auf, erhalten wir nach Einsetzen von T(o) in c- e"!- T(o) die gesuchten Beziehungen fiir die 
Abkihlungskurven ‘ 


a " r cass * ee — . aie 
Oat) = 21537 tee 0 5,04 + e~ 07821088 4 24,19 + Ee 2089+ 107% 
a 2+107 — + 1074 
0,(t) so 12,64 + @ 1:002 + 10 ee 4,46 5 e 9782 10-4¢t =f 24,32°€ 1,089 + 10 fy: 
04(2) —— —— (61a c Ca 1,002+1077¢ ar 18,24 Sea 6,782: 10°4¢ ae 29,25 2 e7 1,089 1074 t 


In Bild 6 sind sowohl die durch Messung bestimmten, als auch die berechneten Erwar- 
mungs- und Abkiihlungskurven eingetragen. Bei dem Vergleich der beiden Ergebnisse ist zu 


= — — — 
30 gat \ 


s gemessen 
——— berechnet 


SS 


Statorbleche 


SS 


~Abkithlung[min } 
; 0 50 00 
Laufzeitt-Anwiirmung| min] 
0 50 100 50 200 
Bild 6. 


bedenken, daB nicht nur die berechneten Kurven mit einem Fehler behaftet sind, sondern auch 
die durch Messung ermittelten. Mit einer naheren Analyse der MeBfehler wollen wir uns hier 


aber nicht beschaftigen, sondern uns mit einer Schatzung zufrieden geben — fiir kleine Zeiten 
rund 7%, weiteres dann abnehmend auf 34%; 
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Da natiirlich die Abstimmung der Widerstaénde entsprechend der Enderwarmung nicht 
vollig eindeutig ist — bei verschiedenen Werten der Widerstande kénnen sich die gleichen 
Enderwarmungen ergeben — sind die Korrektionen in gewissem Grade willkiirlich. Es wur- 
den darum versuchsweise noch einige andere Alternativen in Erwagung gezogen, indem die 
Werte der Widerstande innerhalb gewisser wahrscheinlicher Schranken gedndert wurden und 


deren Lésung dann je nach dem, bessere oder weniger gute Annaherungen ergaben (siehe z. B. 
Bild 7). 


=) \ 

Say GNU | ee \ Anwarmung 

= Za \ (0,026) 

= emessen Ta 10,032 \ 0,02 
30 £ * g : 12 = 0,075 (0,046) 

S NOP aes berechnet og = 0,243 (0,224) 
25 ira 134 = 0,096 (0,125) 


Yb = 0,114 (0,667) 
W, = 160 (160) 
Wo2 = 280 (280) 
W 3 = 165 (150) 


Abkthlung 
111 = 0,400 (0,333 


= ) 
Typ == 0,075 (0,046) 
10 143 = 0,667 (0,476) 
To4 = 0,356 (0,356) 
3g = 1,175 (1,060) 
5 43 = 0,715 (2,500) 
8, = 18,1 (17,9) 
i i 8, = 35,5 (32,7) — 
Laufzeit t{min] ie 
0 50 100 150 200 


Bild 7. 


~ 


8. Berechnung des Erwarmungsverlaufes bei intermittierendem Betrieb 


Soll noch der Temperaturverlauf bei intermittierendem Betrieb berechnet werden, ist zu- 
nachst notig, den Erwarmungsverlauf fiir beliebig gegebene Anfangstemperaturen zu bestim- 
men. Zu diesem Zwecke gehen wir von dem zweiten Ausdruck Gl. (9) aus. Es ergibt sich 


{c: e'T(o OV = Ho 0) wobei H(0) die einspaltige Matrix der Anfangstemperaturen ist, die 
wir nun ganz allgemein auffassen miissen. Da T(o) = c7! (0) ist, erhalten wir ohne Weiteres 


T,(0) = 0,1200 #,(0) — 0,1103 9,(0) + 0,0080 #5(0) , 
T,(0) = 0,6519 0,(0) + 0,0755 8(0) — 0,2069 ¥,(0) , 
T3(0) = 0,2280 #,(0) + 0,0348 3,(0) + 0,1988 ee 
Im gleichen Sinne setzen wir fiir den Abkiihlungsverlauf anstatt 8 (0) = (772 3o58 Ze 47) 
allgemein # (0) = [0,(0) 8,(0) @,(0)] und erhalten aus c T(o) = Ho ) (¢c fir die Abkihlung!) 
z. B. an Hand der CRAMERschen Regel 
T,(0) = 0,0960 8,(0) — 0,0982 #,(0) + 0,0027 8,(0) , 
T,(0) = 0,2280 #,(0) + 0,0230 #,(0) — 0,2077 8,(0) , 
T3(0) = 0,6761 %(0) + 0,0752 (0) + 0,2050 #,(0) . 
Setzen wir fiir T (0) jewels die jetzt berechneten Werte fiir Anwarmungs- und Abkiihlungs- 
vorgang ein, erhalten wir die recht umfangreichen Ausdriicke fiir den berechneten Erwar- 
mungsverlauf bei intermittierendem Betrieb 
0,(0) = 17,77 + [1,08 + 0,1200 8,(0) — 0,1103 8,(0) + 0,0080 (0) ] e~ 290? * 10"! 
= 4,25 -.0,0519 05(0) — 0,0755 0(0) 10,2069 0,(0)] e977" 4°" 
—[14,64 — 0,2280 8,(0) — 0,0348 8,(0) — 0,1988 8,(0)] e— 1°86 20°" 
#,(0) = 32,52— [--- usf., 
wobei wir jeweils bei der Anwirmung und Abkiihlung fiir (0) die bei vorhergehender Ab- 
kithlung oder Anwarmung gerade erreichten Werte einzusetzen haben. 
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Wollen wir uns mit der Berechnung der Endtemperaturen der einzelnen Zyklen begniigen, 
dann setzen wir die Laufzeit (z. B. 12 min) und Stillstandszeit (z. B. auch 12 min) in die oben 
angedeuteten Ausdriicke fiir Anwarmung resp. Abkiihlung ein, so daB wir Ausdrticke folgender 
Art erhalten 


One = 8,385 a> 0,239 Ode mile 0,033 ears a 0,083 een ) 


Oink =-0,772 Ding 0,084 Daye 1 :0;,0025 Usya  1USWe 


Aus diesen sich derart ergebenden Gleichungen, hier sechs, ist es dann auch moglich die 
Enderwarmung bei intermittierendem Betrieb zu berechnen; hier 


eg = 13,4 pe = 14,4 
Cine aa 25,8 ae hs oe = 14,6 a Gee 
93 mx = 33,7 B3mn = 27,0 


In Bild 8 sind die gemessenen und berechneten Erwarmungskurven fiir intermittierendem 
Betrieb verglichen. 


Stator - 
kupfer 


=. Erwirmung(°C] 


Stator- 
eisen 


gemessen 
——— berechnet 


96 —_Loufzeit ¢ [min] 
0 12 2 36 50 60 72 84 100 180 


Bild 8, 


9. Bewertung der Berechnungsergebnisse und SchluBfolgerungen 


Die in den Bildern 6 bis 8 graphisch dargestellten Rechenergebnisse zeigen, daB der, fiir ein 
Ersatznetzbild, das in bezug auf die Enderwarmung angeniéhert richtige Ergebnisse liefert, be- 
rechnete Anwarmungsverlauf in verhaltnismaBig guter Ubereinstimmung mit den Messungen 
ermittelt werden konnte (Bild 6). Dabei ist zu bedenken, daB die Temperatur des Rotors nur 
beilaufig und die des Netzknotenpunktes ,,4‘‘ iiberhaupt nicht gegeben war. DaB dieses Resultat 
nicht auf einem bloBen Zufall beruht, zeigt die bis zu einem gewissen Grade in den Werten der 
Widerstande fiir das gleiche Netzbild unterschiedliche Alternative Bild 7, sowie weitere, nicht 
angefiihrte Ergebnisse. Die Abkithlungskurven, bei denen keinerlei Anlehnung an die End- 
erwarmung moglich ist, stimmen schon weniger gut mit den Messungen iiberein (Bild 6 und 7), 
wenn auch prinzipiell eine gute Ubereinstimmung der Form des Verlaufes nach zweifelsohne 
besteht (Bild 7). 

Desgleichen kann in Anbetracht der Schwierigkeit des Problems das Resultat fiir den 
intermittierenden Betrieb (Bild 8) als recht gut angesehen werden. Aus Bild 8 ist deutlich 
abzulesen, daB die Abkiihlungskurven weniger genau sind und daB, waren sie zutreffender be- 
rechnet, die Ubereinstimmung mit den MeBergebnissen wohl hier noch weitaus besser ware. 
Demgegeniiber weist der aus den durch Messung ermittelten Anwarmungs- und Abkiihlungs- 
kurven auf die iibliche Art konstruierte Erwarmungsverlauf fiir den intermittierenden Betrieb 
3ild gy, zumal beim Eisen, einen im Detail den Tatsachen widersprechenden Verlauf auf, 
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Wie schon in Abschnitt 6 angefiihrt wurde, konnten auch die Werte der Zeitkonstanten 
ftir Maschinen einer fortschreitenden Leistungsreihe gleicher Bauart (Bild 5) recht gut be- 
stimmt werden. Dabei sei nochmals auf die Tatsache verwiesen, daB z. B. bei einem auf nur 
vier Warmequellen vereinfachten Ersatznetzbild fiir das Netz Bild 3 alle drei durch Messung er- 
mittelten Zeitkonstanten mit annehmbarer Genauigkeit rechnerisch bestimmt werden konnten. 


= Erwarmungsverlauf gemessen: 

30S Statorkupfer bei Anwiirmung und Abkiihlung 
2 4 ——-— bel intermittierendem Betrieb 

x s \ — -— Erwarmungsverlauf bei intermit- 
ie tierendem Betrieb konstruiert 


Laufzeit t-Abkithlung [min] 
0 50 100 

Loufzeit t-Anwirmung| min | 

0 50 100 150 


Bild 9. 


Ohne die Schwierigkeiten, die noch zu titberwinden bleiben, bevor von einer richtiggehenden 
Vorausberechnung die Rede sein kann, unterschatzen zu wollen, kénnen die angefiihrten prak- 
tischen Ergebnisse doch dahingehend bewertet werden, daB eine zunachst beschrankte Voraus- 
berechnung an Hand der angegebenen allgemeinen Loésung immerhin im Bereich des Méglichen 
liegt ; beschrankt in dem Sinne, daB die Parameter des Warmequellennetzes, hauptsachlich die 
Werte der Warmewiderstande, die den Ausgangspunkt der Vorausberechnung bilden, derzeit 
noch nicht allein an Hand von Konstruktions- und Rechnungsunterlagen geniigend verlaBlich 
bestimmt werden koénnen. Es kénnen zwar einige einzelne giinstige Ergebnisse nicht als ein 
erbrachter, allgemein giiltiger Beweis hierfiir angesprochen werden, die Beispiele zeigen aber 
immerhin, daB bei der angewandten Methode das thermische Verhalten der Maschine auf dem 
Rechnungswege in den wichtigsten Ziigen erfaBt werden konnte, wie dies an Hand der Erwa- 
gungen Abschnitt 2—5 auch erwartet werden konnte. 

Hinsichtlich des etwas anspruchsvolleren Rechnungsganges, der sich bei der praktischen 
Anwendung ergab, kénnte sicher in Zukunft bei der Erlangung einiger Erfahrung und Ubung 
und bei Benutzung moderner Rechenmaschinen und Modelen [14] soweit mit Verbesserungen 
und Vereinfachungen gerechnet werden, daB die zu leistende Rechenarbeit kein Hindernis fiir 
die praktische Vorausberechnung bilden miiBte. Denn auch wenn die praktische Seite der 
Lésung von Mehrkérperproblemen — und das Warmequellennetz stellt ein Mehrkérpersystem 
dar — sicherlich einer der Griinde war, warum man dieser Lésung auswich, ist dieser Schritt 
doch im Hinblick auf die steigenden technischen Anforderungen sicher gerechtfertigt. 


Nachtrag 


10. Zum Vergleich der exakten und der approximativen LOsung des 
Erwarmungsverlaufs fiir den Stab Bild 2 


Ermitteln wir die exakte mittlere Enderwarmung des Stabes Bild2 (siehe KESSLER [12]), 
erhalten wir 


1 w 1? UH 
= & i Ys EB Mer fee ae a fh 
9 [ 7) dv = 25 [Ay — teh Al; 49= 55 
Fiir A, < a/2 entwickeln wir tgh A, nach den Potenzen von A, und erhalten nach der Addition 
von A, und Division durch A3 


1 6 l Of Tali 1 l U H1\2 
§=wlF-——|1—*(__. \4 | ee [ae 
3AF 5 \3 AF 1 105 \3 ) 
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Setzen wir fiir 7;, = 1/327 F und als Warmeiibergangswiderstand 7, = 1/UH / (in Abschnitt 3 
angefiihrte Werte) als Schatzungswerte ein und wenn wir die gesamten Warmeverluste im Stab 
mit W = wlF bezeichnen, erhalten wir, wenn wir noch den Ausdruck 7;/(7;,, + 7,) vor die 
Klammer setzen 


eae aes f 0,2 (2) | 0,26 (2) — ae dE 
Y. 


Vin + %q% ‘i ii / 
wohingegen wir nach den Regeln der Methode der aquivalenten Warmequellennetze 
W 5a Sen ipa 


es Yen t Va 


“i 


erhalten. 

Die Summe vom zweiten Glied in der Klammer des exakten Ausdruckes an gibt also den 
Fehler an, den wir begehen, wenn wir die angefiihrten Werte von 7;, und 7, anstatt der genauen 
Werte von 7,, und 7, zur Berechnung benutzen. Man sieht, daB z. B. fiir (7;,,/7;) < 0,5 der abso- 
lute Wert des Fehlers schon sehr gering ist. 

Fiir den Erwarmungsverlauf lautet in dem gleichen Falle die Lésung von Gl. (5) (siehe 
KESSLER [12]) nach durchgefiihrter Mittelung 


1 
t pe I Re “a 
B(i) = - | T(x, t) dx = Shien [A, —tgh Ag] 
—— es ae n+1 sit Ore 2 oe. 2 2 
4 oS Age ae ae 4}, 
— UH 2. Sie, j ~ (2n—1)a 
An ae Lancia und O52 


Fiir den zweiten Ausdruck entwickeln wir den Bruch x (eal in eine Reihe. Wir erhalten 
n n 2 

so nach Einfiihrung von W,7,,,7, und (mc) =o-F-+1-c auBer dem schon angefiihrten 

Ausdruck 


—wIF 3. B(x sind, 6 [5 x ( 1 +E —|\ = 


Wi ced oa dis nO! geen ya 
aed ERE Ae ete ee 1 SAD Oe dute Oe ae 
Oe Wes i Ie) 1 olFc I 3 1 
Sal ee ep ga ae (=) 54 i) 162 Sy ee are a er a 
oat los On \ a On \ "a On \%e (mc) \tin 3 Vi 
oder wenn wir noch den Ausdruck 7,/(7,,, ++ 7,) vor die Klammer setzen und = 0,1, 2,...n 
einsetzen 
Vine ‘"y. 7,,\2 a [oss : + | t 
W —* * 410,985 —0,21 (“) + 0,27 (24) i | Toe "in" 
Yin + Ve VG VG 
1 1 
V. v.,\2 Sac [7 peel 
= 11,24 - Lon? — 2) Oye Ome (“*) = 8,91 Ome (2) a Se | 52 in ti 
VG VG 
ay al Al 
| - =o —3 [Vin 5 (Vin 2 i ~ (me) [287 a zal é 
Sia Le Ole LOmes aly Ou lO = | —— FAO LO ea tercone one in ee eer 
VG Ve 
usf. 


Es ist augenscheinlich, daB es nicht mehr, wie im stationaéren Fall, méglich ware, durch 
Einsetzen von richtigen Werten fiir 7,;, und 7, die Reihen in den eckigen Klammern zu liqui- 
dieren und daB die Reihe nach ¢ ebenfalls bestehen bliebe. Es ist jedoch der absolute Wert des 
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zweiten und der weiteren Glieder der Reihe gering, und es schwindet auch der EinfluB der 
zweiten und der folgenden Zeitkonstanten praktisch schon nach kurzer Zeit nach Beginn der 
Anwarmung. 

Die Lésung mittels der Methode der aquivalenten Warmequellennetze ergibt 


1 1 1 
OG) = W Min "ti f et (mc) Ga cs Al | 


Yin VG 


Der Verlauf der exakten Lésung fir 7,, = 0,07° C/W, 7, = 0,20° C/W, W = 800 W und " 
(mc) = 3,104 kWsek/° C ist in Zahlentafel 3 enthalten. Gleichzeitig ist der jeweilige Fehler der 
angenaherten Losung angegeben, sowie der Anteil, der durch den 4°%igen Fehler der Approxi- 
mation der Widerstande in bezug auf die Enderwarmung entsteht. 


Zahlentafel 3. 


¢ (min) ) 2 5 | 10 | 20m AO | 80 | 120 oa) 
Erwarmungsverlauf exakt (° C) fo) 1,86 | 5,94 | 11,15 | 19,1 | 28,7 | 36,3 38,4 | 39,1 
Gesamtfehler (° C) der aquiva- | | | 
lenten Lésung (0) 0,96 0,96 | 1,45 | 2,2 | 2,89 Soules 1,6 
Der von der Approximation von | ~4% 
v;, und 7, herrihrende Fehler | 0,17 | 0,39 0,80 1-2 A eel 1,5 1,6 


11. L6sung des Erwarmungsproblems mit Hilfe von Modelen 
und automatischen Rechengerdten 


Um erstens bei der Wahl des Netzbildes nicht aus Griinden iibermaBiger Rechenarbeit ein- 
geschrankt zu sein oder auch um einfach die Rechenarbeit auf ein Minimum herabzudriicken, 
ist es natiirlich wichtig, sich bei praktischer Anwendung angemessener Hilfsmittel zu bedienen. 
Wir wollen darum hier einige Anmerkungen zu dieser Frage, der natiirlich eine selbstandige 
Arbeit gewidmet werden sollte, beifiigen. Numerische Rechenmethoden, die ebenfalls als 
Hilfsmittel dienen kénnen und die auch beim Matrizenkalkiil in Frage kommen, sind allgemein 
bekannt und wir wollen uns daher mit ihnen nicht naher befassen. 

Zur Berechnung der Enderwarmung kann mit Vorteil ein Gleichstromnetz benutzt werden, 
das ohne besondere Schwierigkeiten realisierbar ist. Dabeiist beim Entwurf eines solchen nicht 
einmal notig, z. B. die Temperaturabhangigkeit der Warmeentwicklung in den, das Kupfer 
darstellenden Warmequellen oder die Abhangigkeit der Kithllufttemperatur #, von der in die 
Kiihlluft abgegebenen Warmemenge usf. durch kostspielige elektronische Hilfsmittel zu 
modelieren. Wie leicht zu zeigen ist, konnen die Werte von W, und #, an Hand der ersten An- 
naherungswerte fiir #,, » = 1, 2,... durch das Model, nach und nach festgelegt werden. Die 
erste Annaherung von #, erhalten wir, wenn wir statt den Unbekannten W, zunachst Wy, und 
schatzungsweise fiir #, einsetzen usw. Dabei geniigt meist schon die erste Korrektion. Von 
der Giite der Ergebnisse kann man sich durch Berechnung der Warmebilanz ftir die einzelnen 
Netzknotenpunkte iiberzeugen. 

Man kann natiirlich auch fiir die Berechnung eine beliebige, zur Losung simultaner Glei- 
chungssysteme entworfene Maschine benutzen, was aber gewisse Nachteile mit sich bringt. 
Der Vorteil des Modelierens des Netzes besteht darin, daB jedes Netzelement direkt einem 
Widerstand, den Verlusten in der Warmequelle usw. entspricht. Es sind darum die am Netz 
durchgefiihrten Operationen anschaulich. Man kann z. B. ohne groBe Miihe die Parameter des 
Netzes so 4ndern — was technisch sehr wichtig ist — daB man fiir gewisse gegebene Elemente 
und Temperaturen leicht die ihnen entsprechenden Werte der restlichen Parameter quasi ver- 
suchsweise findet u.a.m. Einen gewissen Nachteil bilden allerdings die oben erwahnten 
Korrektionen, wenn man mit einem ganz einfachen Netzmodel ohne selbstandige Korrektionen 
fiir « und #, arbeitet. Bei einer Maschine, die fiir die Lésung von Simultangleichungen ent- 
worfen ist, entfallt zwar die Berechnung von Korrektionen, bei der praktischen Arbeit wird 
sich aber der Mangel an Anschaulichkeit unliebsam bemerkbar machen. Auch ist es, sollen 
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beispielsweise mehrere Alternativen gerechnet werden, nétig, immer vom Neuen die Werte der _ 


Koeffizienten der Gleichungen zu berechnen, da z. B. die Hauptdiagonalkoeffizienten des 
Systems von dem Wert aller Widerstande der betreffenden Zeile abhangen usf. Sollen z. B. 
zu einer gegebenen Lésung die Parameter des Problems gesucht werden, st6Bt man so auf 
eroBe Schwierigkeiten. 

Bei der Berechnung des Erwarmungsverlaufs sind ebenfalls zwei Wege gangbar, bei den 
gleichen Vor- und Nachteilen, wie sie eben erwahnt wurden. So wurde die Lésung Abb. 7 mittels 
eines Analog-Rechengerates ermittelt. Es sei bemerkt, daB man, hat die zur Verfiigung stehende 
Maschine nicht geniigend Elemente fiir die Berechnung der Form des simultanen Gleichungs- 
systems Gl. (4), diese Gleichung in eine Gleichung héheren Grades umformen kann, wodurch 
man evtl. mit weniger Verstarkern usw. auskommt. Will man z. B. den Temperaturverlauf bei 
intermittierendem Betrieb ermitteln, so mu8 man zwei Gleichungssysteme — Anwarmung und 
Abkiihlung — in der Maschine haben, so daB man durch Einsetzen der im vorhergehenden 
Halbzyklus ermittelten Enderwarmung den nachsten Halbzyklus ermitteln kann. 

Andererseits zeigte es sich aber praktisch als unméglich, durch Variation der in der Maschine 
eingesetzten Koeffizienten die Werte der Parameter des Problems, einem gemessenen Verlauf 
entsprechend, zu bestimmen. Es scheint, daB dies bei einem eigentlichen Model eher méglich 
ware. 
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Multiflex - Galvanometer 
Empfindlich wie ein Spiegelgalvanometer — ab 
5-10-11A—. Einfach im Gebrauch wie ein 
Zeigerinstrument. Verschiedene Systeme von 
20-6000 Ohm. Austauschbare Lichtzeigerskalen 
mit parallaxenfreier Ablesung. Skalenbeleuch- 
tung umschaltbar fir Netz- und BatterieanschluB. 
Zusatzgerdte: Photographische Registrierein- 
richtung. Photoelektrischer Nachlaufschreiber. 
Photoelektrische Schaltkontakte. 
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Electric Contacts Handbook 


by RAGNAR HoLM, Ph. D., aided by 
Etsr Ho.m, Ph. D., both of St. Marys (Pa.) 
Third completely rewritten edition of ,,Die 
technische Physik der elektrischen Kontakte“ 
by R. Holm 


Mit 194 Abb., XVIII, 522 Seiten Gr.-8°. 1958 
In englischer Sprache. Ganzleinen DM 52,50 


AUS DEN BESPRECHUNGEN 


,,..-Aufbauend auf eigenem Erfahrungsmaterial 
aus mehr als vier Jahrzehnten hat Rk. Holm — mit 
Frau E£. Holm — das heutige Wissen wohl er- | 
schépfend dargestellt, und zwar in einer Weise, 
wie es eben nur der kann, der sich das Gebiet 
selbst erarbeiten muBte. Sorgfaltige Behandlung | 
der bestimmten Werkstoffeigenschaften und | 
vieler physikalischer und chemischer Vorgange 
sind erforderlich, das Verstandnis fiirdie Vorgange | 
an den Kontakten zu schaffen; anscheinend | 
langst verstandene Begriffe, wie Harte, Reibung 
und VerschleiB, erscheinen dabei unter einem 
neuen Licht. . .‘ 
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Bavelemente 


Wir fertigen in unserem Freiburger 
Werk fir die verschiedensten Anwen- | 
dungen der modernen Elektronik: 


Silizium-Transistoren 
Silizium-Dioden 
Silizium-Zenerdioden 

Silizium-Leistungs-Gleichrichter 


Germanium-Transistoren 
Germanium-Flachendioden 
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Gesellschaft fir Metallurgie 
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Freiburg/Brsg., Hans-Bunte-Str. 19 
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C. CONRADTY NURNBERG 


Elektroden, elektrische und galvanische Kohlen 


Der lichtelektrische Effekt und seine Anwendungen 


Unter Mitwirkung von Dr. K. W. BOrr, Berlin, Dr. F. ECKART, Berlin, Reg.-Rat Dr. W. LEo, 
Braunschweig, 

bearbeitet und herausgegeben von Dr. H. SIMON, Professor an der Humboldt-Universitat und stellv. 

Direktor des Instituts fiir Festkérperforschung der Deutschen Akademie der Wissenschaften zu Berlin, 

und Dr. R. SUHRMANN, o. Professor und Direktor des Instituts fiir physikalische Chemie und Elektro- 


chemie der Technischen Hochschule Hannover 


Zweite, vollkommen neubearbeitete Auflage des Buches ,, Lichtelektrische Zellen und ihre Anwendung“ 


Mit 599 Abbildungen. XII, 747 Seiten Gr.-8°. 1958. Ganzleinen DM 97,50 


AUS DEN BESPRECHUNGEN 


,,Die vor 25 Jahren erschienene erste Auflage trug den Titel ,,Lichtelektrische Zellen und ihre Anwen- 
dungen‘‘. Inzwischen hat die Forschung eine Fiille neuer Erkenntnisse gebracht, so daB die nun vor- 
liegende Ubersicht iiber den heutigen Stand einem dringenden Bediirfnis entgegenkommt. Die An- 
derung des Titels soll zum Ausdruck bringen, da das Buch in seiner zweiten, vollig neubearbeiteten 
Auflage weit iiber Photozellen im engeren Sinn hinausgreift. Auf dem Gebiet. des 4uBeren Photoeffekts 
beruhen die in den vergangenen Jahren erzielten Fortschritte auf der Anwendung ,,zusammengesetzter 
Photokathoden und der Legierungskathoden sowie auf der Verstarkung photoelektrischer Str6me durch 
hochemittierende Sekundarelektronen-Kathoden. In Belichtungsmessern hat das Selen-Sperrschicht- 
Photoelement, das sich trotz vieler Neuschépfungen auf dem Halbleitergebiet als unschlagbar erwiesen 
hat, weiteste Verbreitung gefunden. W4ahrend also hier ein gewisser Abschlu8 erreicht zu sein scheint, 
ist auf dem Gebiet des inneren Photoeffekts die Entwicklung noch in vollem Flu8. In den entsprechen- 
den Abschnitten wird deshalb weniger Vollstandigkeit als eine Ubersicht tiber die Grundlagen ange- 
strebt. Dazu werden die erforderlichen Erkenntnisse der Festk6rperphysik in erfreulich klarer Weise 


dargeboten. 


Bei der Fiille des aufzunehmenden Stoffes war es ein gliicklicher Gedanke der Herausgeber, zur Abfassung 
einiger spezieller Abschnitte weitere Spezialisten (K. W. Borr, F. Ecxart, W. Lo) heranzuziehen. 
Das Buch bringt in wohlabgewogener Weise Grundlagen und Anwendungen. Wieweit das von den 
Anwendungen umfaBte Feld sich ausdehnt, kann nur durch einige wahllos herausgegriffene Stich- 
worte angedeutet werden: Sekundarelektronenvervielfacher und Szintillationszahler, Photowiderstande, 
Halbleiterbildwandler, Rontgenbildverstarker, Dosimeter, Ikonoskop und Vidicon. Natiirlich wird 
auch auf die Verwendung von Photoelementen zur Umwandlung von Sonnenstrahlung in elektrische 


Energie eingegangen. Der praktisch erreichbare Wirkungsgrad betragt bei Siliziumzellen 11%. 


Das ganz dem neuesten Stand entsprechende Buch stellt eine wahre Fundgrube dar und kann daher 
warmstens empfohlen werden. Aber nicht nur der reiche Inhalt ist hervorzuheben, sondern in gleicher 


Weise auch die Prazision der Darstellung. Zeitschrift fiv angewandte Physik 
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